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ABSTRACT 


The use of VLBI observations for the estimation of geo- 
detic and astrometric parameters is discussed. Analytic models 
for the dependence of delay and delay rate on these parameters 
are developed and used for parameter estimation by the method 
of weighted least squares. 

Results are presented from approximately 15,000 delay and 
delay-rate observations, obtained in a series of nineteen VLBI 
experiments involving a total of five stations on two continents. 
The rms scatter about the weighted mean in the estimates of 
baseline lengths can be summarized as follows (Haystack (Mas- 
sachusetts) is the reference station] : 
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The station in Alaska had by far the poorest system sensitivity 
and that in Sweden the next poorest. The corresponding scat- 
ter in the source coordinate results is less than 0?1 for all 
coordinates except the declination of the low-declination 
sources, which have a scatter of the order of 0',*2. The closure 
of baseline triangles is investigated and found to be consis- 
tent with the scatter of the various baseline-component results. 
Estimates are made of the wobble of the Earth * s pole and of the 
irregularities in the Earth's rotation rate. The differences 
between our estimates of the pole position and the BIH estimates 
have a scatter of about 1 m in both the X and Y components , and 
the corresponding (Al-UTl) differences have a scatter of under 
3ms. Estimates are also made of the precession constant and of 
the vertical Love number, h, for which a value of 0.55±0.05 was 
obtained. 


Thesis Supervisor: Irwin I. Shapiro 

Title: Professor of Geophysics and Physics 
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CHAPTER 1 


INTRODUCTION 


One of the main themes of astronomical research through- 
out history has been the desire to construct instruments 
capable of increasingly finer angular resolution. In the 
field of radio astronomy this desire has culminated with the 
development of very-long-baseline interfercmetry (VLBI) . A brief 
sketch of the history of radio astronomy should help to put 
this development into perspective. 

The science of radio astronomy got its start in the 
1930 's with the discovery by Karl Jansky of extra-terrestrial 
radio signals. Jansky, \;ho was then working at the Bell 
Telephone Laboratories, built a directional radio antenna 
for the purpose of locating sources of radio interference. 

This antenna had an angular resolution of about 30®. Even 
at this level of resolution Jansky was able to recognize 
that some of the sources of the signals he received were extra- 
terrestrial, and he was able to identify the newly discovered 
galactic center as the source of some of the signals. 

It was obvious that if much progress were to be made 
in studying celestial radio sources, it would be necessary 
to improve the resolving power of the receiver systems. In 
particular, in order to compare radio observations with opti- 
cal observations, it would be desirable to have a resolving 
power in the radio range comparable to optical telescope reso- 
lution, Increasing the resolving i>ower of 


1 




! 


radio antennas to this level was to prove a dif- 
ficult task, for the angular resolution in radians of any re- 
ceiver system is roughly equal to the ratio cf the wavelength 
of the signals received to the diameter of the antenna system. 
The wavelengths of the radio signals being employed ranged 
from a few centimeters to several meters, whereas optical wave- 
lengths are typically 5 x 10 m. Thus, in order to make a 
parabolic antenna to receive 3 cm radiation with the same re- 
solving power achieved by a 15 cm diameter optical telescope, 
the radio antenna would have to be 10 km in diameter. build- 
ing such an antenna would be a herculean task, even if the an- 
tenna were to be fixed and not steerable. These difticulties 
were sufficiently obvious to make early researchers in radio 
astronomy turn to interferometry as the only practical way to 
achieve resolution porfonrance con^^rable to tiiat of optical telesco^xis. 

With interferometry, angular resolution is achieved by 
comparing signals from widely separated antennas. The limiting 
factor on the resolution is then the distance between the an- 
tennas rather than simply the size of any single antenna. 

These early interferometers needed a real-time communication 
link to transmit signals from one antenna to the other. Techni- 
cal problems with the communication link limited the antenna 
separations to tons of kilometers. The caj^abilities of these 
interferometers still provided a dramatic improvement over 
the capabilities of a single antenna, but it was clear that 

still better icsolution could be attained if the communication 

1 2 


link restriction could be overcome. 

The problem of eliminating the communication link was 
solved with the development of very-long-baseline interferometry 
(VLBI) . With VLBI, the need to transmit signals from one 
station to the other in real time is bypassed by the use of 
highly stable atomic clocks at each station. The clocks are 
used to control the local oscillators which heterodyne the 
signals down to a frequency which can be sampled and recorded 
on a standard digital tape recorder system. (The clocks also 
control the sampling equipment and the tape recorders, but it 
is in the control of the local oscillators that the frequency 
stability of the clock is most important.) Thus instead of 
transmitting the radio signals from one antenna to the other 
in real time, it is only necessary to transport the tape re- 
cordings to a common location at leisure. With this develop- 
ment, the deployment of radio interferometer stations was 
limited only by the size and shape of the Earth. 

There was an added, though not unexpected, bonus 
which resulted from the development of VLBI. The interfero- 
meter observations could be used to estimate the antenna separ- 
ations as well as the positions of the radio sources observed. 
Antenna separation can be estimated from short-baseline inter- 
ferometer data as well, but in the case of short-baseline in- 
terferometers, the accuracy of the antenna position determina- 
tion is not substantially better than that obtainable with 


conventional survey techniques. 
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Over the longer baselines used in VLBI work, the conven- 
tional survey techniques are much less accurate than they are 
over shorter baselines. Thus the development of VLBI promised 
an improvement in the ability to determine both radio source 
positions and radio antenna locations. The purpose of this 
thesis is to report on a series of such determinations. 

The work that led up to this thesis began in 1968 with a 
series of experiments involving stations in Tyngsboro, Mas- 
sachusetts, and Green Bank, West Virginia. These earlier ex- 
periments are discussed in Ffeferences 5, 6, and 21 • The 
experiments reported on in thi s thesis were carried out between 
April 1972 and January 1975 under the sponsorship of the 
National Science Foundation, the Advanced Research Projects 
Agency of the Department of Defense, and the United States 
Geological Survey. The stations employed were located in 
Tyngsboro, Massachusetts; Goldstone, California; Green Bank, 
V7est Virginia; Gilmore Creek, Alaska; and Onsala, Sweden. 

In the following Chapter, I will discuss the nature of 
the VLBI observables and the limitations of using them to estimate 
the astrometric and geodetic quantities of interest Chapter 3 
will deal with the analytic models used to represent the VLBI 
observables in the computer programs used to analyze the VLBI 
data. Chapter 4 concerned with tlie results of an analysis of 

a large body of VLEJ observations spanning the period from 
April 1972 to January 1975. Finally, in Chapter 5, I will discuss briefly 
the proutur*'. 1 ;; for future' {nogress in the development of VLBI. 
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VLBI OBSERVABLES AND PARAMETER ESTIMATION 

2.1 INTRODUCTION 

As considered in this thesis, the fundamental VLBI ob- 
servable is the measurement of the difference in arrival times 
at two stations of a radio signal from a distant radio source. 
This quantity is referred to as delay, and its time derivative 
is referred to as delay rate. Suitable combinations of such 
measurements can be used to estimate the locations of the 
stations and the radio sources. The nature and limitations of 
such estimates are the subjects of this chapter. 

2.2 THE DELAY OBSERVABLE 

In order to understand the delay observable, we construct 
the following idealized picture. Consider a point source of 
radio signals sufficiently distant that the signals arrive es- 
sentially as plane wave fronts, and consider two stationary 
receiving stations, capable of detecting and recording these 
signals. The stations are also equipped with accurate clocks 
capable of recording the arrival times of these signals (Fig- 
ure 1) . In the sinplest case, with a non-rotating Earth, perfect clocks and 
receivers, and signals of arbitrary bandwidth, the difference 
in arrival times of the signals could in principle be determined 
by simply recording each signal on a strip chart recorder (also 
of arbitrary bandwidth) and then comparing the recordings to 
determine what shift of time origin will make the signals match. 
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Figure 1 

Basic VLBI Geometry 
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Real VLBI observations differ from this simple picture prin- 
cipally by having a limited receiver and recorder bandwidth, 
by having Doppler shifts in the signals due to the relative 
motions of the receiving sites, and by having imperfect clocks 
at each of the sites. Nevertheless, this picture is adequate 
for understanding much of the discussion which follows, and is 
a useful preliminary for understanding those situations where 
a mere detailed model must be employed. 

In the presence of a dispersive propagation medium ( i .e , 
a medium in which the velocity of the incoming signal is a 
function of the frequency of the signal) there will be two 
characteristic velocities associated with the radio signals — 
the group velocity and the phase velocity. There will there- 
fore be two characteristic delays, the group delay and the 
phase delay. Our data-reduction programs are designed to es- 
timate both the phase delay and the group delay, as well as the 
phase -delay rate. In theory, the difference between the phase 
delay and the group delay could be measured and used to es- 
timate the magnitude of the effect of the dispersive medium. 

It has not been possible with the data being discussed in this 
thesis to resolve the closely spaced (-0.127 ns) ambiguities 
in the phase- delay data, and therefore the phase-delay in- 
formation has not been used (see Reference 21/ pps. 21-25). 

There are also ambiguities in the group-delay data, but 
the group-delay ambiguities are much more widely spaced 
(typically j usee to 1 psec -- see Table 2 ) , and can be 
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eliminated from the data (see Chapter 4). 

In the remainder of this chapter it will be assumed, 
unless specified otherwise, that the propagation medium is 
non-dispersive and that therefore no distinction need be made 
between phase delay and group delay. 

2.3 THREE COMPONENTS OF DELAY 

It is useful to think of VLBI delay observations as com- 
posed of three components. The first component is due to the 
geometry of the station locations and source location, as 
illustrated in Figure 1. This component contains all of the 
information about the geodetic and astrometric parameters 
that are of interest.* The second component is due to instru- 
mental effects, principally clock errors. The third component 
is due to effects of the propagation medium ( e . g . the atmos- 
phere and the ionosphere) . The detailed theory and computer 
algorithms required to deal with each of these components will 
be discussed in Chapter 3. 

2.3.1 Geometric Delay 

In order to understand the use of delay observations 
for the estimation of geodetic and astrometric quantities, we 
need to consider the geometric component of the delay observable. 


ft 

Many of the radio sources that we observe are known to have 
structure on the order of 10“3 arc-sec. Structure on this 
scale could affect the delays on transcontinental baselines 
by about 0.1 ns. The effects of source structure have not been 
considered in the discussions in this thesis since other un- 
certainties dominate our results. 
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From Figure 1, vre see that tJie arrival time of a signal at Station 2 
minus the corresponding arrival time at Station 1, i,e. the geonetric 
delay, can be vnritten as 



where § is the unit vector in the direction of the radio source, 
B IS the vector separation of the receiving stations, and c 
is the speed of light. (In the remainder of this thesis, the 
unit of time will be taken to be the second, and the unit of 
length will be the light-second, with these assumptions, 
c = 1, and it will not appear explicitly in equations such as 
Equation (2-1) . For « discussion of the specific coordinate 
system in which 2 and ^ are defined for this thesis, see 
Section 3.1.) Now for an inertial coordinate system we will 
assume that d is fixed, and that B will, to a very good approx- 
imation, rotate with a constant angular velocity 2, the angu- 
lar velocity of the Earth, is therefore a diurnal sinusoid: 

Tg K^sin(nt + (J)) + (2-2) 

Because at present there is no available method for setting 
the clock epoch and rate (see Section 2.4) at each station 
to the accuracy required for VLBI observation, two additional 
terms must be added to Equation (2-2) to represent the dif- 
ferences in the clock epochs and rates at the two sites: 

Tg = Kj^sin(f^t + 4>) + K 2 + Cq + c^t (2-3) 

Thus delay observations on a single source can be used to 
determine only four independent parameters ( e .g. 

(K 2 + Cq) , and c^; ft is assumed to be known). The parameters 
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that we wish to determine from these observations are; three 
components of the baseline, two source coordinates (components 
of 6), and the two clock parameters, Cq and Of these 

seven parameters, only six can be estimated, as the right as- 
cension of the source cannot be separated from the hour angle 
of the baseline. There are still six parameters to be estima- 
ted, and only four independent observations. This is clearly 
an under-determined case. In other words, VLBI observation of 
a single radio source cannot serve to determine all of the 
baseline and source coordinates that we wish to estimate. How- 
ever, if we observe more than one source, each source adds 
three new independent observations (K^^, K 2 and 0 for that 
source; and c^ are common to all sources) and only two new 
unknown parameters. It is thus easily shown (see Reference 17) 
that observations on a minimum of three sources are necessary 
to determine all of the parameters of interest (3 baseline 
parameters, 5 source parameters and 2 clock parameters, in 
that case). Furthermore, since there are ten unknown parame- 
ters, there must be at least ten observations in all, and each 
source must be observed at least three times. In actual 
practice these minimum numbers of sources and observations 
are generally greatly exceeded in order to provide redundancy 
and internal consistency checks on the data. 

If the observations made are delay rates rather than 
delays, a slightly different situation emerges. If we resolve 
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B into equatorial and polar components in Equation (2-1) , 

'g “*®EQ • ® ® ‘2-4) 

and differentiate this equation, noting that S was assumed 
to be rotating with an angular velocity n parallel to Bp , we 
see that 

% = 0 (2-5) 

and 



Thus delay- rate observations have no sensitivity to the 
polar component of the baseline. 

If we differentiate Equation (2-3) , we get 

T = flK,cos(()t + <|)) + c, (2-7) 

Following an analysis similar to that for delay observa- 
bles, we note that there are four unknown quantities which 
we wish to estimate from observation of a single source (two 
baseline parameters, one source parameter and one clock para- 
meter) , and there are only three independent parameters in 
Equation (2-7). If we observe more than one source, each 
new source adds two unknowns, but only two independent ob- 
servations ( e.g. and c<^ is common to all sources). Thus 
it is not possible with delay rates alone to estimate all of 
the parameters of geodetic and astrometric interest. The polar 
component of the baseline cannot be estimated, and, in ad- 
dition, at least one of the baseline coordinates or source 


21 


I 

f 


T 


coordinates must be determined by some other means. 

If the observations made are a combination of delays and 
delay rates, one might expect that the minimum number of sources 
to be observed would be less than the number required for de- 
lay data alone. This is not the case, however. If there are 
enough delay observations to determine the delay vs . time 
sinusoid, then the de’ay rates are also determined, and measure- 
ments of delay rates do not then add independent information. 
(The situation is the same as for observations of delay beyond 
the minimum number required to determine the delay vs . time 
sinusoid.) 

In all of the VLBI data that will be discussed in Chap- 
ter 4 the minimum criteria for the estimation of all the base- 
line coordinates and source coordinates will be more than met. 

In general, the data sets will include data which span a period 
of one to four days, involving two to four stations and ten to 
twenty radio sources. Also, both delay and delay-rate ob- 
servations will generally be included simultaneously in the 
solutions. 

There are other parameters which affect the geometric de- 
lay, besides the coordinates of the baselines and of the 
sources. These other parameters concern effects of the non- 
rigidity of the Earth and irregularities in the Earth's 
angular velocity vector, ti. These parameters are of interest 
in geophysics. Because these effects modify the curve defined 
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the arguments of the previous paragraphs 


by Equation (2-2), 
concerning the information content of a set of delay observa- 
tions on a single source need to be modified. That is to say, 
the delay vs . time curve for a given source contains informa- 
tion about the flexure of the Earth and the irregularities in 
S in addition to infonnation about the three geometric para- 
meters discussed. The analytic expressions for these irregu- 
larities and their effects on delay and delay rate will be 
discussed in Chapter 3. 

2.3.2 Instrumental Delay 

The instrumental delay can be thou^t of as consisting 
of two components. The first component is the time required , 
for the signals to travel from the antenna to the device 
which samples and records the signals. We will refer to this 
component of the instrumental delay as the "cable" delay, al- 
though the signal passes through amplifiers, mixers and other 
active components in addition to cables. The second component 
of the instrumental delay is caused by the error in the system 
clock at each station (VLBI is sensitive, of course, only to 
the differences between the clock errors at the various 
stations) . 

In an ideal VLBI station the cable delay should be care- 
fully measured and eliminated from the data. However, since 
the equipment necessary to make such measurements was not 
available during nearly all of the experiments reported on in 
this thesis, the cable delay must be dealt with. The cable 
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delay could be either dispersive (frequency dependent) 
or non-dispersive , and it could be constant or time- 
varying, If it is non-dispersive, then its effects will 
be indistinguishable from the effects of the clock errors. 
If it is dispersive, then it could well affect the group 
delay and the phase-delay rate differently. In the 
absence of direcc cable-delay measurements as a function 
of frequency, dispersive cable delays will have to be re- 
garded as a possible source of systematic error. 

We have attempted to model the errors in the system 
clock at each station in terms of a polynomial in time 
whose coefficients will be estimated from the VLBI data. 

If we refer to the clock error at Station i as then* 


"i 


+ B^t + a^t' 


( 2 - 8 ) 


where the first term will 


In the remainder of this thesis it will be recognized 
that, in the absence of measurements of the cable delay, an 
unknown portion of the apparent clock error is caused by 
(non-dispersive) changes in cable delays. 
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be referred to as the clock-offset error, and the second term will 
be referred to as the clock-rate error. Higher order terras 
will be referred to as "second-order terra", "third-order terra", 
etc., as needed. (The reader should be aware that an n order 
polynomial generally contains n + 1 coefficients, and that, for 
exeunple, the coefficient of the second-order terra is actually 
the third coefficient in the polynomial.) 

In the hydrogen- maser clocks generally employed in VLBI 

work, the clock error consists predominantly of the clock-of f i'.et 

and rate-error terras. If the hydrogen maser is operating close 

14 

to Its purported stability limit of one part in 10 , then 

the effects of error terras of second-order and higher should be 
less than one nanosecond in a day. Therefore the total in- 
strumental delay in a well designed and operated VLBI system 
should be able to be modeled to the level of a nanosecond per 
day by determining the coefficients of the offset and rate terms 
of the clock-error polynomial. 

In the data to be discussed in this thesis there arc many 

instances of clocks whose stability was no better than a part 

12 . • 
in 10 . The common types of clock errors include discontin- 

uities in both offset and rate, and long-term drifts in the 
clock error due to the presence of significant high-order deriva- 
tives. The problems caused by the high-order derivatives of the 
clock error could sometimes be handled by fitting high-order co- 
efficients of the clock-error polynomial. The discontinuities 
were handled by fitting a separate clock-error polynomial to 
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the data on either side of the discontinuity. The epochs at 
which the VLBI theoretical model changes from one polynomial 
to the next will be referred to as "breaks" in the c.lock-error 
polynomials, ’niere is nc^v no provision for constraining clock 
error polynomials to be continuous in any derivative across 
a "break", although this feature could easily be added. 

Tiie selection of the clock parameters for a set of VLBI 
observations with very badly behaved clocks is a problem whose 
solution is entirely up to the discretion of the person pro- 
cessing the data. The judicious resolution of this question 
and the determination of clock models which accurately leflect 
the behavior of the real clocks is probably the thorniest 
problem in the entire process of analyzing the VLBI data re- 
ported here. This problem will be discussed at some length 
in Chapter 4. 

There is a technique which can be used to eliminate 
clock problems altogether. If there is more than one 
antenna at each station, observations can be taken on more 
than one radio source simultaneously. If such observations 
are differenced the clock errors are eliminated from the ob- 
servations. This type of experiment has been performed and 
described elsewhere (References 3 and 7}. If there is only one 
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antenna at each station, then observations cannot be made 
simultaneously on different sources. However, if the antennas 
can be switched with sufficient rapidity from one source to 
another, then it is possible to form differences of pairs of 
observations which are nearly simultaneous. However, sinc',‘ 
the observations being differenced are not exactly simultaneous, 
the effects of short-term clock drifts (that is, clock drifts 
vM.ch build to a significant level during tlie time interval between the two 
observations vhich have been differenced) will still be present in the data, 
but the effects of longer term clock drifts will be largely eliminated. This 
type of experiment has also been performed and is discussed elsevdiere (Ref. 23) 
Some theoretical work we have done has indicated that in- 
clusion of observations on a large number of sources may allow 
the recovery of all of the di jnal geodetic and astrometric 
parameters from a set of pure differenced observations (see 
Appendix A) . While none of the results discussed in Chapter 4 
involved exclusively differenced data, differenced data were 
used to handle several experiments involving more than two 
stations, where the clocks at one or more of the stations were 
very badly behaved. The data on all the baselines involving 
a station with a badly behaved clock were differenced, and the 
data on the remaining baselines were left undifferenced. One 
advantage to this scheme is that the unditferenced data are 
sufficient to determine all of the source coordinates. The 


differenced data are then more than sufficient to determine 
the coordinates of the stations with the faulty clocks. 

It is possible, of course, to use more general linear com- 
binations of observations, rather than restricting ourselves 
to differenced observations. If, for example, the observations 
which have been differenced are also summed, it is then possible 
to form a data set consisting of sums and differences of obser- 
vation pairs. Such a data set will contain all of the informa- 
tion of the original data set, but with the novel property that 
the long term clock errors have a much greater effect on one 
half the data set (the sum data) than on the other half (the 
difference data) . It should be possible to use this combined 
data set with the estimated standard errors for the summed 
data adjusted to reflect the presence of the systematic clock 
ei’rors, and thereby improve our estimates of geodetic and astro- 
metric parameters. The implementation of this technique is 
beyond the scope of this thesis. 

There are also techniques for dealiiig with clock errors 
which involve the use of Kalman filters to develop stochastic 
models for clock errors. These techniques are also beyond the 
scope of this thesis. 

2.3.3 Propagation- Medium Delay 

The remaining component of the delay observable which 
must be dealt with is the part due to the propagation medium 
( e.g. the atmosphere and ionosphere) . The daytime ionosphere 




28 


zenith correction at X-^band is typically of the order of one 
nanosecond, and the nighttime zenith correction is typically 
an order of magnitude smaller. We have experimented with a 
model for the ionosphere which assumed that the ionospheric 
electron content followed a rectified diurnal sine wave with 
a maximum at local noon. The correction to the delay was then 
calculated with a simple ray-trace scheme. Since this model 
did not affect our results significantly, either in terms of 
the rms of the residuals or of the scatter in the baseline 
results, it was not used in obtaining the results discussed 
in this thesis, and will not be discussed further. The 
ionosphere must be regarded as a source of possibly systematic 
noise at the level of a fraction of a nanosecond. It should 
be remarked that it is possible to exploit the frequency- 
dependence of the ionospheric correction to eliminate this 
source of noise completely. By comparing measurements made 
simultaneously at two or more frequencies it should be possible 
to deduce the ionospheric contribution to the delay and elim- 
inate it. While multi-frequency observations have been proposed 
for future experiments, the necessary multi-frequency receiver 
hardware was not available for any ot the experimeiits reported 
here. 

The neutral atmosphere has proven to be somewhat more 
tractable than the ionosphere. The zenith correction is typically 
about 7 nanoseconds and a fairly simple tncdel has been shewn to produce 
a noticeable improvement in the postfit residuals. 
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The primary reason we have had success with a simple model 
for the neutral atmosphere is that the neutral atmosphere is 
fairly stable. The main cause of both spatial and temporal 
variations in the neutral atmosphere is the highly variable 
water-vapor content of the atmosphere. But the water-vapor 
contribution to the atmospheric correction to the VLBI delay 
observable is believed to be no more than 10-15% of the total 
atmospheric correction, or, in other words, about the same as 
the ionosphere contribution in the daytime. The atmosphere 
model we have employed uses the VLBI data to estimate the 
magnitude of p^, where is the correction to the delay in 
the zenith direction at Station i. The theoretical value for 
each delay observation is corrected by p^^csc 9^ at each station, 
where is the elevation angle for the source at Station i. 

(The actual model employs a small correction to the csc9^ term 
which is not important for the discussion in this chapter; the 
details are discussed in Chapter 3.) There are three fairly 
obvious weaknesses to such a model: first, the modified cscG^ 
law may not adequately model the elevation-angle dependence 
of the real atmosphere; second, there is no provision for 
azimuthal variations iii the atmosphere; and third, there is 
no provision for time variations in the atmosphere. The third 
objection can be ameliorated to some 
degree by solving for separate zenith 
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thicknesses on separate days or even at shorter time intervals. 
The selection of the time intervals over which a given value 
of p^will be applied is the second major variable (after the 
determination of the clock model) which is at the discretion 
of the experimenter. 

2 . 4 PARAMETER ESTIMATION 

The station locations and source coordinates are estimated 
from the VLBI observations by using the method of weighted least- 
squares, if the noise in the data is gaussian with zero mean, is 

equivalent to the method of maximum liJcelihood. The procedure 
employed is to calculate the theoretical value of the 
observable using the best available values of the parameters 
to be estimated. The weighted-least-squares or maximum like- 
lihood estimates of the parameters are then found through 
an iterative procedure based on linearized 

equations. Since this technique is quite standard, it will not 
be treated further here. The observation "weight", or formal 
standard error, used in the weighted-least-squares analysis is 
based on the signal-to-noise ratio for the observation (Refer- 
ence 2 2 ) . The term "residual" will be used to refer to the 
difference between an observation and its corresponding theore- 
tical value. The formal errors of the parameter estimates are 
calculated from observation formal standard errors that have 
been scaled so as to make the rms value of the weighted per* t fit 
residuals equal to unity. 
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(CRT) display for read out and a light gun which makes possible 
real time interaction with the computer through the CRT dis- 
play. Thus it is possible to use the light gun to set up the 
clock and atmosphere models and select the parameters to be 
estimated and then have the computer solve the linearized 
least-squares equations, print out. the parameter estimates, and 
graph the postfit residuals on the CRT display. The data can 
then be edited with the light gun, and a different clock model 
or other change in the parameter set can be set up and the 
solution can be run again, all in a matter of a few minutes. 
While this program cannot do cinything which cannot be done in 
some fashion by the batch-processing program, the advantages of 
easy set-up and fast turn-around are enormous. It is pos- 

sible with this program to explore more completely the range 
of possible clock models and other variations in the parameter 
set, and therefore to explore possibilities that otherwise 
would be neglected for want of time. Another advantage of 
real-time processing is that the CRT display of the plot of 
postfit residuals can be manipulated. Whereas in a batch- 
processed program the plotting scales must be pre-selected 
( according to some criterion such as maximum residual size) , 
with real-time processing the plot can be manipulated and 
scales can be expanded or contracted to display all or part of 
the residual set to better advantage. (These plots can be 
photographed, and such photographs are displayed in Chapter 4.) 
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The Haystack processing program has some limitations 
which ought to be mentioned. At present, it cannot handle 
difference and sum data, discussed in Section 2.3.2. Also the 
procedures for sorting and selecting data to be processed are 
restrictive; for instance presently all data are required to be 
in strict time sequence. A more fundamental limitation on this 
program arises from the fact that the computer's memory is too 
small to accommodate easily the programs necessary to calculate 
the theoretical values and partial derivatives required to form 
the linearized least-squares equations. In order to by-pass this 
problem, the theoreticals and partials are calculated by the 
IBM computer program and stored on magnetic tape for use by the 
CDC program. The inability of the CDC program to calculate 
theoreticals restricts the program to a single iteration of the 
iterative solution to the least-squares equations. Fortunately 
a single iteration has generally been found to be sufficient 
to produce convergence of the parameter estimates to within 
an error very much smaller than the formal errors of these para- 
meter estimates. 

The analysis upon which our data processing programs are 
based is discussed in Chapter 3, and Chapter 4 is devoted to a 
discussion of the use of these programs and to the results of 
processing the VLBI data. 

2.6 OBSERVATION SCHEDULES 

Before concluding this discussion of baseline and source- 
coordinate estimation, it is necessary to discuss the scheduling 
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of observations. The observation schedule is a factor which 
affects the process of parameter estimation, yet the problem of 
optimizing a schedule for astrometric and geodetic parameter 
estimation has rot yet been solved. A complete discussion of 
the process of creating an optimal schedule is beyond the scope 
of this thesis. There is one scheduling consideration, however, 
which proved to be of some importance in the analysis of the 
data to be discussed in Chapter 4. It is very useful to have 
interleaved observations on a large number of sources, in order 
to separate the effects of the clock errors from the effects of 
errors in the source coordinates. Many of our schedules did 
not call for such interleaved observations because the schedules 
were usually designed for the purpose of studying the structure 
of the radio sources being observed, rather than for the pur- 
pose of estimating baseline coordinates and source coordinates. 

Creating a schedule can be a tedious and trying project 
even without the problem of having to optimize the schedule for baseline and 
source-coordinate estinvation. At the very least, it is neces- 
sary to calculate which sources are visible at all sites at any 
given time, and how much time is required to switch each antenna 
from one source to another. In order to automate the routine 
calcula*" ions involved in creating a schedule, an operator- 
interactive computer program was written for the CDC 3300 com- 
puter at Haystack. With this program, the operator first 
nrlr’iUn < ho r. 1 .H i<MV4 ami st'iirc'os tio wishrs 


to include in a 


4 


f 


schedule^ The computer then displays on the CRT the sources 
which are visible at all stations at any given time. For 
each source the computer also displays such pertinent informa- 
tion as the time required to move the antennas from one source 
to another, the minimum elevation angle of the source, the time 
elapsed since the source was observed and the total number of 
observations already scheduled for that source. The operator 
need only select with the light gun the sources he wishes to 
observe, and the computer will automatically increment the time 
by the appropriate minimum time required to move the antennas, 
and record the source to be observed at that time, and then dis- 
play a new list of sources which are visible at all stations at 
the new time. This process is repeated until the schedule is 
complete. The program also has options which allow the operator 
to modify a schedule already created or to print a final 
schedule in a form to be used by the antenna operator at each 
site. 
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AU»IYSIS OF VLBl OBSERVABLES 


3.1 INTRODUCTION 

We wish to derive theoretical expressions for the dependence 
of the VLSI delay and delay-rate observables on parameters which 
are of geodetic and astrometric interest. We will then be able 
to use these expressions to estimate the values of these para- 
meters from the VLBI data using conventional least- squares tech- 
niques . 

3.2 DEFINITION OF COORDINATE SYSTEMS 

Before continuing with a discussion of delay and delay- 
rate observables, we must discuss the coordinate systems which 
will be employed in most of the analysis which follows. There 
are two coordinate systems which must be considered, one which 
is fixed to the "rigid" Earth and one which is "nearly" xnertial. 
We will defer discussion of the Earth-fixed coordinate system 
until Section 3.5. 

The "nearly" inertial coordinate system we will use will be 
defined with its origin at the solar-system barycenter, its 
Z-axis will be defined as parallel to the Earth's mean pole at 
1950.0 (positive north), its X-axis will be defined rn the 
direction of the ascending node of the Earth's orbit at 1950.0 
(the "first point in Aries") , and its Y-axis will be defined to 
complete a right-hand triad. This coordinate system will be 
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employed for all of the analysis in Chapter 3 unless otherwise 
specified. (As in Chapter 2 , the unit of length will be taken 
to be the light-second, and the unit of time to be the second, 
so that c = 1.) This coordinate system is not quite inertial, 
since it is subject to accelerations arising from the gravitation- 
al attraction of the nearby stars and from the orbital motion 
of the solar system about the center of the galaxy. The motion 
relative to the galactic center introduces a 

change in the apparent position of the radio source with a mag- 
nitude of about 10 ^ seconds of arc per year. The effect of the 
gravitational attraction of the nearby stars is several 
orders of magnitude smaller than the effect of the galactic or- 
bital motion. Both of these effects are small enough to be neg- 
lected in our treatmer u of VLBI observables. 

The reason we have chosen to use a solar-system barycentric 
coordinate system for our analysis here is that we hope at some 
point to be able to combine VLBI observations with spacecraft 
tracking data and interplanetary radar data in order to estimate 
the motions of the Earth more completely than can be done with 
VLBI data alone. Combining various types of data in this manner 
is simplified if the analysis of all of the data types is done 
in the same coordinate system. 

In Appendix B we deal with the problem of transforming time 
as kept by atomic clocks on the Earth to coordinate time, t, of 
our solar-system barycentric coordinates. We can re-write 
Equation (B-30) as 
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hT . •- t - 32.15 sec *• K • r + LPT (3-1) 

1 y 

where AT^ is he n luing of a (peru^cL) atomic clock at Station i c.* the 
Earth, P is tiv. position of the center of the Earth relative 
to the solar system barycenter, r is the position of the clock 
relative to the ceiiter of the Earth and LPT represents the 
long-period (non-diurral) terms in Equation O-l).* The reason 
for grouping the long-period terms in this fashion wjll become 
apparent later. The 32.15-sec term is introduced in order to 
have coordinate time agree as closely as possible with the pre-* 
vious definitions of ephemeris time (Reference 1, p, 15). 

Equation (3-1) represents the reading of a "perfect” 
atomic clock located at Station i. Recognizing that ''real" 
clocks are necessarily imperfect for a variety of unpredictaole 
reasons, we will model the errors of each real clock by a poly- 
nomial whose coefficients we will attempt to estimate from the 
VLLI observations. If T.,, represents the reading of the clock 
at Station i fo ■ observation K, then we write 


iK 


M I 
1 


^ “iK + 


6^^(t 


. ) 
'iK^ 






(3-2) 


The superposed dot, as in "R" signifies differentiation with 
respect to coordinate time only. The sole exception to this 
convention will be in the expression for delay rate, t, where 
the superposed dot signifies differentiation with respect to 
time as kept by an atomic clock on the Earth, as specified in 
Equation { 3- 11 ) . 
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where t.„ is the defined coordinate-tiine origin for the clock poly- 

XK 

ncmial. Equations (3--) and (3-2) can be re-written to give t in the 
following useful form 

t = T.j^ + 32.15 + S . - LPT - - 6.j,(t - t.^^) - ... 

‘(3-3) 

3.3 THE DELAY OBSERVABLE 

If we ;,hink of the observed radio signal as a plane wave 
(Figure 1) which is received at Station 1 at coordinate time t, 
and received at Station 2 at coordinate time t + At, then we 
can define the observed value of the delay as: 

iK(t) - T2j^(t + At) - (3-4) 

At is simply the geometric (coordinate) time delay 

At(t) =-B(t)*d + At^ (3-5) 

wher^ 

-B(t) = l5(t) + r^it)] - lR(t+At) + r^(t + At)) (3-6) 

and 

§ 5= unit vector in the direction of the radio source 
At^ = atmospheric (including ionospheric) delay 

where the ionospheric correction is presumed to be appropriate 
for group or phase delay depending on which type of delay is to 
be modeled. 

A simple method for using Equation (3-5 ) to determine 

At is to perform a regula falsi type of numerical iteration. 

Equivalently, one can expand the expression (Reference 14) , neg- 

*The orr«)rs introduced hero and elsewhere by adding vectors in 
n :\i i' tuM]lit|ihlc foi 1 h»' purposes of this thesis. 
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lecting terms smaller than 10 sec in delay. We can 
express the site position vector at Site 2 as 

S(t + Lt) + ?2(t + At) = S(t) + 5(t)At + |i(t)At^ 

+ ?2 (t) + r2 (t) At 

+ ^r2(t)At + ... 


Then Equation (3-5) can be re-written as 

At(t) = (r3^-r2)-d - (R+r2) *^At - (R+r2) + A- 

At(t) • {l+(^+r2) •§} = {rj^-r2)*6 - j(5+r2) (?j^-r2 

+ AT^ 

At(t) = (rj^-r2)*6 - [ (^+^ 2 ) I (rj^-?2) -6) 

. d-(R+r^).d] - i[(R+r^).^] 

H 




where the first term on the right hand side of Equation (3-8) 

_2 

is of the order of 10 seconds or smaller, and the second term 

“4 

is smaller by a factor of the order of 10 , and the third 

term is smaller by a factor of the order of 10 The part 

of the second term inside the third set of brackets consists 

of two terms, the first of which is unity and the second is 

-4 

smaller than unity by a factor of the order of 10 

Now if we similarly expand the first term in Equation (3-4) 
as 

4 - -»> 

T^K^t+At) = t + At - 32.15 sec - R(t)*r 2 (t) 

-CS(t)-r2(t) + ;l(t) -rjCt) lAt + LPT(t) 

+ g^[LPT]At + 


(3-9) 
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we can then combine Equations (3-1)# (3-2), and (3-9) with 
Equation (3-4 ) , using At of Equation (3-8) to obtain 


+ (LPT) 'At - 




(3-10) 


v^ere the first term of the right-hand side of Equation (3-10) 

-2 

IS of the order of 10 sec or smaller, the second term is 

-4 

smaller by a factor of the order of 10 , the third term is 
smaller than the first by a factor of the order of 10 
and of the two terms inside the brackets of the second term, 
the first is of the order of 10 sec and the second is 

3 

larger by a factor of the order of 10 . The magnitude of 

the clock-error terms depends, of course, on the behavior of the 

clocks. 

It is obvious from Equation (3-10) that and a 2 j^ can- 
not be determined separately from observations of Tj^ alone. 

It might seem that the higher order terms, could be 

determined for both stations (i=l •■'nd i=2) because of the 
type of term. However, for most VLBI observations $ 21 ^ ^ 

and At ~ 10 sec, thus this term has a magnitude of less than 
-15 

10 sec, which is well below the noise level in the measure- 
ment of Tj^. The contribution of the At term for the higher-order 
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coefficients is even smaller. Therefore observations of can 
be used to determine only the differences between the clock 
error coefficients at the two stations [ e . g« (c*2k ~ °^1K^ ' etc] . 
The only reason for writing two polynomials explicitly in 
Equation (3-10) is for convenience in handling the case where 
more than two stations are involved. In the procedures we em- 
ploy in processing '?LBI data, the values of the coefficients of 
the clock error polynomial at one arbitrarily selected reference 
are fixed at zero, and the polynomial coefficients at the other 
stations are understood to represent the differences between 
those coefficients for that station and those for the reference 
station. 


3.4 THE DELAY-RATE OBSERVABLE 

In order to discuss the delay-rate observable, t, we must 
first define which of several possible time scales has been used 
to form the derivative of x. The choice is somewhat arbitrary. 
We have defined the delay rate to be the derivative of delay 
with respect to time as kept by the clock at Station 1. We can 


write this explicitly as 
; (t) = ^ dt 


(3-11) 


Now it is a simple matter to differentiate Equation (3-10) to 
obtain (Reference 14) 
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dt 


A*t “ - S- (r2”r^) - (5* r2+2R* r2+J * ^ 2 ) At 

2 


-(R*r2+RT2) At + ^->2-(LPT)At + ^(LPT) A*t 

dt 


"^IK " - ... + $2^(1+ A*t) (3-i: 


+2y^^(t-t2y^*lit) (1+ t^'t) + ... 

where A*t is obtained by differentiating Equation (3-8) : 

• • • • • • 

A*t = (rj^-r2> •© - I (R+r2> •§! I (rj^-r2) 'e) 

-( (S+r2> -S] [ (rj^-icj) 

+ [(rj^-r2)-e] [(R+rj) ’e]^ 

+ 2[ (rj^-r2) -S] ( (R+r2> -6] ( (K+rj) -§] 

-((S+r2)-i)((?i-r2)*^H(?l-?2>'®l + 

and from Equations (3-1) and (3-2) 


= [1 - R-?j - R-Jj^ + ^(LPT) + + 2Y,v(t-tiv) + 


IK 


IK"' ‘'IK' 


« 1 + R-rj^ + S-r^ - ^(LPT) + ... 


(3-1 


Then, combining Equations (3-14), (3-11), and (3-12), and 
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keeping terms larger than 10 


we write 



(LPT) 1 


dT 


T (t) = 


K 


Et 


“ at" 


or, writing 
10-16 


T(t) = 


dx .. * • 

g^[R.ri+ R.rj^- 


d 

dt 


A*t(l-J^+ R.?^) - A‘t I^U-PT) 

<3t„ 

— ^ explicitly, again deleting terms smaller than 

• • • • * t ** 

A*t ~ [R- (r2“r^) + R* {r2'r^) 1 • (1+ A*t ) -R- ?2 * At 


" \k ■ ■ ••• 

+ 2Y2K(t-t2j^+At) (1+ A't) (3-16) 

where the first term is of the order of 10 ^ or less; the second term is 
smller by a factor of about 10 and the third term is smaller by a factor 
of about 10**^^, As in the case of the delay observables, only the differences 
between the coefficients of the clock-error polynanials (e.g. , 
estimted usefully from the VLBI delay-rate observable. 

3.5 CALCULATION OF THE SITE POSITIONS 

In order to evaluate Equations (3-10) and (3-16) , the position R of the 
center of the Earth and its time derivatives are calculated fron a tabulation 
on negnetic tape of an ephemeuris of the positions of all of the bodies in the 
in the solar system. 

In order to evaluate r (t) , we must begin by defining an 
Earth-fixed coordinate system in which we will express our 
station coordinates and baseline components. The Z-axis of our 
Earth-fixed coordinate system will be defined in the direction 
of the mean pole for 1900-05 (see Reference 10, p. 93-95). 

The X-axis will be defined perpendicular to the Z-axis in 
the direction of the Greenwich meridian. The Y-axis will be 
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defined to complete a left-hand triad. (The coordinate 
system is defined to be left-handed in order to be consistent 
with the spherical system involving radius, latitude and west 
longitude, which is also left-handed.) 

The transformation from Earth-fixed coordinates to coordin- 
ates aligned with the 1950.0 system described in Section 3.2 is 
done as follows j If we let [u] represent the 1x3 matrix 
containing the components of r^^in Earth-fixed coordinates, 
and (r^l represent the 1x3 matrix containing the components 
of r^^ in 1950.0 coordinates 


(r^)’’ = [P] IN) [S] [W] (u)’' 

where [Pi is the precession matrix, [N] is the nutation matrix, 
[SJ is the rotation matrix for the Earth's spin, and [W] is the 
rotation matrix which expresses the polar motion of the Earth. 
^htrix [P] has the following form (see Reference 10, p. 93-98): 


[PI = 


cos CqCOS w cos z cos ^qCOS w sin z 


-sin^QSin z 


+sin CqCos z 


-sin CqCos w cos z -sin^^cos w sin z 


-cos ^Qsin z 


-sin w cos z 


+COS ^qCOS z 


-sin w sin z 


cos ^Qsin w 


sin ^Qsin w 


cos w 
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where 

90®-^q = the right ascension of the ascending node of the 
mean equator of date in the 1950.0 coordinate 
system 

90®+2 = the right ascension of the ascending node of the 

mean equator of 1950,0 in the coordinate system 
defined by the mean equinox and equator of date 
w = the inclination of the equator of date with res- 

pect to the 1950.0 coordinate system. 

The values of z, and w are obtained from the series; 

= 2304V948T + 0’.'302T^ + 0'.'0i79T^ 

2 = 2304V948T + 1V093T^ + 0V0192T^ (3-17) 

w = 2004V255T - 0';426T^ - 0'.'0416T^ 


where 


T = time from 1950.0 measured in tropical centuries 

(36524 . 21988 days of 86,400 coordinate-time seconds each) 


Matrix [N] has the following form: 

{ 

1 COSE A\|; sin e 


[N) = 




-Ai{i cos e 
-A\|; sin e 


1 

-Ae 


Ae 

1 


where 


A\{; = the nutation in longitude 
Ag = the nutation in obliquity 
G = the true obliquity of date of the ecliptic 
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Alp and Ae are calculated from the trigonometric series given in 
Reference 4, pp. 44-45. In order to save computer time, these 
series are calculated and tabulated on magnetic tape at half- 
day intervals, and interpolated by a fourth-difference Everett 
interpolation to the desired epoch (see Reference 1, p. 224). 
The true obliquity of the ecliptic e is given by 

e = Cq + Ae (3-18) 

where, by Reference 4 , p. 98, the mean obliquity of the eclip- 
tic is 

Cq = 23®27'08V26 - 46V84T* - 0V0059T*^ + OV0018lT'^ 

(3-19) 

with T' being the time in Julian Centuries of 3^525 ooordinate-time 
days from the epoch 1900 January 0.5 (J.E.D. 2415020.0). 

The spin matrix (S] has the following form 


[S] 


where 


cos 6 -sin0 0 

sin 0 COS0 0 


]0 0 1 


0 = 0^ + Alp cose (3-20) 

- ^®0 

®o = ® ^ uTi 

0=6" 38"^ 45?836 + 8 , 640 , 184?542T ' + 0?0929T'^ 
d0 

—St ^ 1.002737909265 + 0.589 x 10"*‘-^T' 


where the last two expressions are evaluated at 0.00 UT for 
the day in question and UTl, the (mean solar) rotational 



phase angle of the Earth is related to atomic time through the 
tabulated values of (UTC-UTl) published by the BIH, vthere UTC is a 
defined piecewise linear function of atomic time, A1 (see 
Reference 1, pp. 26-29). Alternatively, the value of (UTC-UTl) 
can be estimated from VLBI data. For historical reasons, 
however, we present in Chapter 4 (Al-UTl) rather than (UTC-UTl) . 

It should be noted that the expressions in Equations (3-17) through 
(3-20) were employed by the BIH in deriving the tabulated values 
of (UTC-UTl) and the values for polar motion, discussed below, 
from photographic zenith-tube observations. Therefore these 
expressions must be employed in order to use the tabulated 
values of (UTC-UTl) and polar motion to reconstruct the angular 
orientation of the Earth as measured by the BIH. 

The polar motion matrix has the following form: 


(W) - 


where: 


== the component of the geocentric angular position, in 
radians, of the pole at time t with respect to the 
mean pole of 1900-05 measured along the meridian 
toward Greenwich 

= the corresponding angular component along the meridian 
90®W of Greenwich. 


1 


X and y are obtained from the tabulated polar-motion 
values published by the BIH. Alternatively we can use VLBI 


data to estimate these values. Both the tabulated values and 
the estimated values for X and Y are corrected for diurnal 
polar motion, as calculated by McClure (Reference 8, p. 83). 

The time derivatives of r^ are calculated by taking suc- 
cessive cross products of the Earth's spin vector, 5, with r^^^, 
where the magnitude of S, is calculated by differentiating 
Equation (3-20) 

ae 




dt 


de 


gO l^(UTC-UTl) + cos c) 


(3-21) 


The derivatives in the second and third terms are calculated 

by numerically differentiating the tabulated values for UTC-UTl 

dUTC 

and Aii^cose. In the first term we have assumed that — = 1. 

This assumption introduces an error in site velocities of the 
-16 

order of 10 sec/sec. Also, and 0 have zero slope, except 

for discontinuities at the beginning of each UTC day, and so 
their derivatives do not appear in Equation (3-21) . 

3.6 ANTENNA-MOTION CORRECTIONS 

Up to this point we have treated the VIBI stations as if each station 
were a single point fixed on a presuned rigid Earth. We must new deal with 
the fact that the receiving antennas are, in general, steerable paraboloids 
^lAiose dimensions are on the order of thirty to sixty meters. These antennas 
are steered or "pointed" by rotating them about either of two mutually ortho- 
gonal axes. Cne of these axes is always motionless with respect to the Earth's 
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crust, and will be referred to as the fixed a>:is. The other axis rotates 
about the fixed axis (Figure 2) . On sane antennas the fixed axis is aligned 
with the local vertical (alt-az mount) ; on others it is aligned 
horizontally (X~Y mount) ; and on still others it is aligned 
with the Earth's spin axis (equatorial mount, as in the figure). 

The analysis of the corrections needed to account for the 
antenna motions will be greatly simplified by the recognition 
that some components of the needed correction are constant, 
independent of the direction in which the antenna is pointing. 
Such components have no effect on the delay rate, and will af- 
fect the delays only by changing the apparent clock offset. If 
the VLBI observations are to be used to synchronize the clocks 
at two stations, then these constant components must be measured 
and removed from the clock-offset estimate. Otherwise, such 
constant components will not concern us. 

In order to begin to analyze ♦-he effects of antenna motion 
on the VLBI observables, it is necessary to carefully define 
what point near the moving antenna will be considered as the 
end of the baseline, or baseline reference point, in the sense 
that the baseline between two antennas will be defined as the 
vector separation of the baseline reference points of the two 
antennas. The correction to the delay for the motion of the 
antennas has a particularly simple form if the baseline refer- 
ence point is defined appropriately. Consider the moving 
axis, which rotates about the fixed axis. Since the axes are 
orthogonal, the moving axis is confined to a plane perpendi- 
cular to the fixed axis. The baseline reference point will b-:’ 
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defined as the point where this plane intersects the fixed axis. 

Now in order to calculate the contribution of the antenna 
motions to the geometric delay, we need to carefully define the 
point on the antenna where the signal will be said to be "re- 
ceived", This point will be taken to be the antenna feed horn 
since the component of the delay arising from the travel time 
of the signal beyond the antenna feed horn has been dealt with 
in Section 2.3.2 (the "cable" delay). Now in Figure 2, it 

can be seen that since the antenna is pointed at the source 
during any observation, the delay caused by the separation of 
the antenna feed from the moving axis will have the same con- 
stant value for all sources. Therefore we need consider only 
the component of the geometric delay caused by the offset, 6, 
between the axes. If the axes intersect, i . e . if |d| = 0, the 
entire effect of the motion of the antenna will have the same 
constant value for all sources. If |d1 ^ 0, the changing correction 
to the delay will be simply the component of D in the direction 
of the radio source 


At 


|'i51oos(^-t(/ 2) = 1^1 sin (J) 


{ 3-22) 


axis offset 

where (J* is the angle between the direction to the source and 
that of the fixed axis of the antenna (see Figure 2) . 

If the antenna mounting is equatorial, that is, xf the 
fixed axis is parallel to the Earth's spin axis, the correction 
takes a particularly simple form, since 

6 = tt/2 - (J) (3~23) 

where 6 is the declination of the source. 
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The correction to the delay rate for the motion of the 
antenna is, of course, the derivative of Equation (3-22), i,e. 


At . ® Id cosd)^^ 

axis offset ' ^dt 


(3-24) 


3.7 PROPAGATION -MEDIUM EFFECTS 

The corrections for the propagation medium are evaluated 
separately for the neutral atmosphere and for the ionosphere, 
the latter being the only appreciable charged-particle contri- 
bution to the values of the observables. For the neutral atmos- 
phere we utilize a modified cosecant law: 


. a a 

“ ■'2 ■ 


where 


a r . . . 0.00143 ,-l 

T. = Pilsine^ + ta H e . V 1 0 ^ 45 ^ 


(3-25) 


(3-26) 


and is the elevation angle of the source as ob served from Sta 
tion _i. The numerical constants in Equations (3-26) represent 
the best-fit values of the results from ray-tracing through a 
standard model atmosphere; the results of the fitting process 
showed that Equation (3-26) agreed with the delays obtained from 
ray-tracing to within 1% for all elevation angles above 1® 
(Reference 2). The coefficient as can be seen from 
Equation (3-26), represents the zenith electrical path length 
(or delay) for the atmosphere over Station Only the water- 
vapor content of the neutral atmosphere presents a serious prob- 
lem for the modcsl because of its variability; its average con- 
tribution to the 7 nsec zenith electrical delay is about 


( 

I 


10-15%. 

The ionosi lere is more difficult to model accurately than 
is the neutral atmosphere, partly because the electron density 
of the ionosphere changes, usually by a factor of about ten from 
day to night. Fortunately, for the X-band frequencies involved 
in these experiments, the total effect of the ionosphere on 
:ielay is never more than about one-tenth of the effect of the 
neutral atmosphere. Since we have not used a model for 
the ionosphere (see Section 2.3.3), seme of the ionospheric correction 
will be absorbed in the neutral atmosphere correction, thereby leading 
to estimates for the neutral atmosphere correction which are 
larger (for group delay data) than the "true" neutral atmosphere 
correction . 

3.8 EARTH TIDES 

To account for the elastic properties of the earth, the 
position of the station at each end of a baseline is modified to 
include the effects of the solid-earth tides. In particular, 
we use the model given by Melchoir (Reference 9 ) . The changes 
in the spherical coordinates of the site position are given by 

Ar = hu 

A* = f (3-27) 

=r ^ 9U 

r cos(P B\ 

where r, (|>, X denote, respectively, the geocentric radial, 
latitudinal, and (east) longitudinal position of the site, and 
where h = 0.584 and i = 0.045 are the relevant Love numbers. 
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The disturbing potential u is given by 
2 

u = A{cos <|) cos(2LHA t + a,2 cos (2LHA-u+w+4»2) 

+ a^cos (2LST+\^2^ "*■ a^cosl2 (GST+X)+4>^] ) + sin2(t) (a^cos (GST 

+X+\p^) tagcos (LHA-y+^pg) +a^cos (LST-o+ij;^) ) - (1.5sin^^-0.5) 
• [ag+a^cos (p-u3)+a^QCOS (2y) +a^j^sin (2o) ] } (3-28) 
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and 

y *= longitude of the moon (mean) = (4.719967 
+ 8399. 709T) radians 

a ** longitude of the sun (mean) == (4.881628 
+ 628.3319T) radians 

(i) ** longitude of the lunar perigee = (5.835152 
+ 71.01803T) radians 
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T = time in centuries from 1900.0 


GST = Greenwich sidereal time 

LHA = lunar hour angle = GST - p + X 

LST = solar hour angle = :iST - a + X 

This treatment of the Earth tides neglects possible inhcnogeneities 
in the response of the earth to the tidal potential, which would 
make h and Z functions of location on the Earth. According to 
Melchior (Reference 9, p. 300) , measurements from different 
parts of the Earth yield values for h that differ froni the ones 
employed here by 0.08, and values for I which differ by 0.03. 




CHAPTER 4 


RESULTS FROM VLBI OBSERVATIONS 


4 . 1 INTRODUCTION 

This chapter consists of a discussion of the results from 
a series of VLBI experiments which began in April of 1972. For 
convenience, these experiments have been divided into two dis- 
tinct sets which I call the Phase I experiments and the Phase II 
experiments. The Phase I experiments were made during the per- 
iod f rc m April 1972 through May 1973, and the Phase II experi- 
ments were made during the period from July 1973 through Jan- 
uary 1975. A brief summary of the pertinent facts about the 
experiments performed is given in Tables 1-3. 

4.2 THE PHASE I EXPERIMENTS 

Sane of the Phase I experiments were undertaken for the purpose of 
measuring transcontinental baselines and others for the purpose of studying 
the structure of the radio sources (see, for example. Reference 21) . Each 
experiment consisted of approximately twenty- four hours of dbservations, 
during vhich time approximtely two hundred observations were made on about 
ten different sources (an observation consists of a single pair of Mark I tape 
recordings: see Reference 21. For experiments involving N stations, there- 

fore, there are N(N-l)/2 dbservations that result from the pairwise ccrobination 
of the N tapes that contain recorded signals from a given 
source at a given time ) . The schedules were carefully designed 
so that, with few exceptions, three or four sources were always 
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kept under observation in fairly rapid succession (<1 hour for 
observations on three sources) . In addition, the observation 
schedules often provided for as wide coverage in hour angle 
and declination of the observed sources as was practical in 
order to optimize the measurement of the delay vs , time sinus- 
oids as discussed in Chapter 2. The schedules were restricted 
in that no sources were observed whose declinations were higher 
than 40®. At the time, no sources were known at higher declina- 
tions which could be reliably used for VLBI observations. Several 
have since been uncovered by our group. Also, the schedules used 
in April and May of 1972 called for a large amount of time for 
observation of the source pairs 3C273-3C279 and 3C345-NRA0512, 
for reasons unrelated to the geodetic purposes of the experi- 
ments. 

4.3 THE PHASE II EXPERIMENTS 

The Phase II experiments began with observations in 
July 1973, and have continued to January 1975. 

There are two characteristics of these experiments which have 
made it convenient to discuss them separately from the discussion 
of the earlier experiments. First, these experiments have gen- 
erally run for a period of four days each rather than one day 
each, and they have generally employed four stations (Hay- 
stack, Goldstone, NRAO and Onsala) rather than two, although 
there is a small number of three station experiments in both 
data sets. Thus, the data from each Phase II experiment com- 
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prised about 2,000 observations (Mark I tape pairs) rather than 
the about 200 observations for each of the Phase I experiments. 
Second, each of the Phase II experiments was scheduled primarily 
for the purpose of studying the structure of the radio sources 
(see Reference 2 3) . The schedules for the Phase II observations 
have therefore tended to concentrate observations on one or two 
sources over periods of time extending for many hours. The 
lack of interleaved observations on many different sources 
tends to make it difficult to separate the effects of clock 
errors from the effects of source coordinate errors, baseline 
errors, and atmospheric effects. Although this problem might 
not have been serious had the clocks at of the stations been 
well behaved, the fact is that the clocks at two of the stations < 
NRAO and Goldstone, were demonstrably very badly behaved during 
most of these experiments. In spite of these problems these 
data proved of some use for geodetic and astrometric measure- 
ments. 

4.4 DATA -PROCESSING PROCEDURES 

Before discussing the specifics of the results of the ana- 
lysis of these data, it will perhaps be useful to discuss the 
general procedures employed in processing each data set. 

The first step in processing the data involves the use of 
two computer programs, VLBIl and VLBI2. These programs are used 
to estimate the "best" values (in the maximum likelihood sense) 


of delay and delay rate for each pair of tape recordings, one 


tape from each of two stations, containing the radio sigiials re- 
ceived from a given source. The development and use of these 
programs is described elsewhere (References 5 and 21). 

4.4.1 Ed iting the D.ita 

The next step in processing any of these data sets in- 
volves the elimination of the group delay ambiguity { see Ref- 
erence 21 ) , In order to detect these ambiguities, VLBI3 is 
used to calculate theoretical values for each observation and 
then to form the residuals or differences between the observed 
values and their corresponding theoretical values. These 
"theoreticals" are based on accurate a priori values for the 
baselinci coordinates, the source coordinates and the neutral 
atmosphere, as well as on preliminary estimates for the co- 
efficients of the offset and rate terms in the clock error poly- 
nomial based on the unedited VLBI data. Ihe delay residuals frcm such a 
corputer run are shewn in Figure 3, and the rate residuals in Figure 4*. 

Most of the observations in Figure 3 cluster alom straight lines, which 
lines are separated by a group delay ambiguity spacing, which in 
this case is 1/3 ysec. The process of eliminating these ambi- 
guities involves nothing more than adding or subtracting N 
times the ambiguity spacing for each point, so as to make all of the 


*If the reader is not already familiar with these plots, the 
explanation in Appendix C will help in understanding the codes 
employed . 

f. 


\ 


I 


7 


SM9 


HS delays, in ns 

9 


Oi 

ca 


UNO 



2.IM 


e 


in 

> 

< 

Q 


2.9M 


S.IN 


e* 


S.9N 


4.100 



-21H. 


•«.0 


Mgure 3. Residuals 



Figure 4. Residuals from Unedited Delay -Rate Data 


parallel lines in Figure 3 coincident. It makes little dif- 
ference which line is chosen ( i.e. which line is not moved at 
all) since this choice will affect only the value of the clock 
offset. There are only two cases where the selection of the 
total ambiguity correction is important. First, if more than 
two stations are involved, the ambiguities must be selected so 
as to preserve closure of the delays around any triangle*. 
Second, if the VLBI observations are to be used to synchronize 
the clocks at the two stations, the ambiguity spacing must be 
large compared to the a priori errors in the more conventional 
methods for synchronizing clocks ( e.g. broadcast time signals 
or portable clocks) . These errors are typically a fraction of 
a microsecond. A second technique for selecting the correct 
£unbiguity could involve the use of an instantaneous bandwidth 
wide enough to resolve the group -delay ambiguity (Reference 21, 
p. 108-110) . The required instantaneous bandwidth would be 
only slightly greater than that obtainable with the present 
recorder systems. A third technique could involve suitably 
changing the separations of the sampled frequency bands 
(Reference 21) , 

It is obvious from Figure 3 that there are some points 
which do not cluster along lines separated by the group-delay 

ambiguity spacing. Th ese observations must be edited out. 

*It should be obvious from the definition of delay (Equa- 
tion (3-4)) that, for Stations 1, 2, and 3; 
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Commonly such "bad” points can be nore easily identified in the delay- 
rate residual plot, where there is no question of ambiguities 
(Figure 4 ) . 

It may well be asked at thi^ point, what causes so many 
"bad" observations, and what criteria are used to eliminate 
them? The answers to the first question are myriad. Bad data 
points can be caused by any failure in the receiv**i -recorder 
hardware at either site, including such things as bad magnetic 
tapes. The only criteria used to eliminate these "bad" obser- 
vations are the magnitude of the residual, and/or its failure to 
lie "close" to a "smooth" curve passing through other residuals 
for the same source. In general, I have tried to be judicious 
in eliminating all suspicious or marginal observations. The 
data sets we have here are so thoroughly over-determined that 
the effect of not including an observation which is, in fact, 
valid will in general be negligible, whereas the effect of 
including a "bad" observation could be significant. Although 
it is possible in principle to investigate all of the myriads 
of combinations of data sets formed by including some data points 
and excluding others, the astronomical number of such possibili- 
ties boggles the mind. .T have therefore sought to eliminate 
as far as possible the "bad" and "marginal" observations, and 
deal only with a single sub-set of each data set which seems to 
contain thoroughly reliable data. In general, about 20 to 25% 
of the observations have been eliminated. Of these observations, 
about 10% were eliminated because the radio sources were too 
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We seek next to deduce the behavior of the station clocks. 
As mentioned in Chapter 2 , if the clocks are not well behaved, 
the determination of the clock behavior can be the most dif- 
ficult problem in the entire process of reducing VLBI observa- 
tions. The best clues to the behavior of the real clocks come 
from the residuals from the type of computer run we have just 
discussed. To preserve these clues, it is desirable to fix 
all of the non-clock parameters because these parameters have 
non-zero correlations with the clock parameters. Therefore if 
all the parameters were adjusted, the mindless least-squares 
estimator would attempt to reduce the effects of the clock on 
the sum of squares of the residuals by adjusting the non-clock 
parameters. This could destroy, or at least confuse, the evi- 
dence of the clock behavior contained in the residuals. For 
the purpose of illustrating this procedure I have included in 
the figures some residual graphs which show some typical effects 
of clock problems. Figure 3 shows the effects of a discontinuity 
in the clock rate at Onsala. (The determination of which one of 
the two station clocks is causing the problem seen in any resi- 
dual set can be difficult. In this case the problem was solved 
by an examination of the station log books, which contained a 
note indicating that the Onsala hydrogen maser clock was re- 
tuned at the time indicated by the break in the residual trend 
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Figure 5. 13elay Residuals Showing the Effect of a 
Discontinuity in the Clock Rate at Onsala 



in Figure 5. In the multi-station case it is often possible to 
determine which stacion clock is causing a particular trend in 
the residuals by noting which station is common to the base- 
lines whose residuals exhibit that trend.) Figure 6 shows the 
effect of estimating the coefficients for a separate first-order 
clock-error polynomial at Onsala for each of the straight lines 
in Figure 5. There is no discernible evidence of any remaining 
clock problems. Figure 7 shows the effect of high-order deri- 
vatives in the clock error at NRAO. This effect would be disas- 
trous for even meter-level geodesy were there no way to correct 
for the effect — its magnitude is nearly 500 ns, and one meter 
is only 3 ns of light-travel-time. Figure 8 shows the same 
data, but this time the coefficient of the second-order term 
in the clock-error polynomial has been estimated as well. The 
parabola shown in Figure 7 is gone, but the remaining drifts are 
still about 40 ns in amplitude. How well these problems can be 
dealt with will be discussed later on in this chapter. 


4.4.3 Testing for Consistenc 




HS 

delays, in ns 

0 10 73 


CODE- 75125- 

9 



® 0 0 « 

’ 9 





e 


d d 0 0 9 

0 


ii 



e 


• 

0 




g 

■!®‘ 


t t ^ to 





1 


CO 

> 

< 

Q 


|o 


2.ID0 


2J00 


S.fOO 

e 

5.|b0 

I 

4.1D0 


i^DO 


d « d 


o<^ 


d d 

0 d ® 


. / V ^ r 

gt fo V 


ad 


J 


0 ^ 


<V 


t i t A ° 

t 'jd*'" . - 
t s 


odJ ^ d 


V. t ♦ I * • 

' fl! 

9 ’Ai j 


Vjo„. ■* •'•'j J ii J j 


d a q 
q a 


t dt : ,v s 

,l V ' . VV , V 


0 4 ^ 


1 0 "l 

0 0 * tt^ 

' { 

d o<\j 


0 0 


## 


J j > 




“ V/4^r 


t t d 


9 9 


*»v 


* f ir 


r r 


V V 


099 


J J 


-!7.42 -15.94 -10.45 -«.970 -5.4* 0 5.4* 9.970 10.45 15.94 17.42 

Figure 6. Delay Residuals from the Same Data as Figure 5, with Separate 
First Order Clock Polynomials Before and After the Rate Discontinuity at 
’ 2; 1 .5 days. 


DAYS 





-2M.5 -taS.T -129.1 -64.97 .0000 64.97 129.1 199.7 2M.9 SZ2.9 



Figure 7. Delay Residuals Showing the Effects of Significant 
High Order Derivatives in the Clock Error at NRAO 
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final step is to check the consistency of the various para- 
meter solutions. It should be recognized, however, that con- 
sistency of results is a necessary condition but not a suf- 
ficient condition for confidence in those results. Convergence 
toward an incorrect result is a well known phenomenon in the 
history of science and it usually results from the fact that a 
series of investigations has been prematurely term.inated when 
those investigations seemed to produce satisfactory results*. 
Nevertheless, the importance of consistency cannot be over- 
stated, for the idea of repeatability of results is deeply in- 
grained into the very concept of measurement, any kind of 
measurement. One expects that use of the concepts of consis- 
tency and repeatability, combined with careful checking and a 
healthy attitude of skepticism will produce reliable results. 

4.4.4 Processing Problems with the Phase I Dat a 

All of the Phase I results which will be discussed in 
Section 4.5, with two exceptions, were obtained from solutions 
in which the clock errors were modeled by estimating the co- 
efficients of the first and second terms of the clock- error 
polynomial at each station except Haystack. (As discussed in 
Chapter 3, the coefficients of the clock -error polynomials at 
Haystack were fixed at zero, and the polynomials at each of 
the other stations were interpreted as the differences between 


This phenomenon is sometimes referred to as the "sheep effect". 
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the clock errors at that station and the clock errors a^ Hay- 
stack.) The two exceptions are the February 4, 1973 data 
(Run 75125-19)* which required the estimation of the coefficient 
of the second-order term in the clock error polynomial, and the 
June 27-28, 1972 data (Run 75125-16)*, which seemed to have some 
inconsistency between the delay and the delay-rate data of the 
type described in Chapter 2. The solution rep>orted in Sec- 
tion 4.5 is a delay-only solution with a clock break at 
1:04 (UTC) on the 28th of June 1972. The residuals are shown 
in Figures 9 (delay) and 10 (delay rate) . The delay 
residuals have no particular bias, whereas the delay-rate resi- 
duals are strongly biased and have a discontinuity at the point 
where the clock break occurs (between 1.0 and 1.125 on the ver- 
tical scale) . It is arguable that these data should simply 
have been eliminated from further discussion, but I have inclu- 
ded them, if for no other reason than to serve as an example 
of the type of problem that can plague these data. The residuals 
on the remaining Phase I data sets exhibit no problrms compa '- 
able to the June 27 data, but neither ire they completely 
trouble-free. Figure 11, for example, she**' the lelay resMualb 
from the August 29-30, 1972 da^'a set. '' e rosiuuals tor 3C 120 
(e) have been connected by c. line so that they may be more 


See Appendix C for an explanation of the run code numbers. 
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easily seen. It is obvious that these residuals are systema- 
tic, yet all attempts to formulate a reasonable explanation 
for this behavior have failed. System problems can be largely, 
^”t not totally , ruled out because the problem does not seem to 
affect any other source even though observations on several 
sources are interleaved in this time period. (Source-dependent 
system problems are not completely impossible; they are only 

extraordinarily unlikely.) Source structure would have to 

-2 

exist on a scale of about 10 arc-second to cause an effect 
of this magnitude. The trouble with a source structure expla- 
nation is that the fringe amplitude data for this day give a 
strong indication that the unresolved component of the source 
is point-like (Reference 16). The remaining obvious possibility 
is that this residual curve represents some sort of atmospheric 
effect. The trouble with an atmospheric effect explanation is 
that none of the other sources, some of which are at comparable 
elevation angles, show any similar effects. While it is not 
possible to rule out pathological situations such as small clouds 
which just happen to cover this source and follow it through the 
day, such explanations have to be regarded as extremely unlikely. 
We are left with reasons for disbelieving all of the likely ex- 
planations for this problem. It is possible, and unfortunately 
even probable, that the delay and delay- rate data simply do not 
contain the information required to unambiguously resolve this 
problem. 

Not all of the peculiar residual effects which are seen 
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defy explanation. Figure 12 shows Haystack -Goldstone delay 
residuals for May 9-10, 1972. The residuals for 3C345 (o) have 
boen connected by lines. It can be seen that the scatter in 
these observations starts at about ±3 ns, and then drops steadi- 
ly until the scatter is a small fraction of a nanosecond. The 
reason for this peculiar behavior is that the structure of the 
source can cause large changes in the correlated flux as the 
projection of the baseline changes (see Reference 23), and, hence, 
in the signal-to-noise ratio for the delay observation. For 
this particular set of observations the correlated flux was 
near zero at the start of the observations, and increased rapid- 
ly throughout the observation period. The decrease in the 
scatter of the residuals is nearly proportional to the increase 
in the signal-to-noise ratio. 

4.4.5 Processing Problems with the Phase II Data 

As was mentioned earlier, the processing of the Phase II 
data was a good deal less straightforward than the processing of 
the Phase I data. The main problem with the Phase II data lay 
in the fact that the clocks at NRAO and Goldstone were very 
badly behaved. A typical effect of the clock error at NRAO can 
be seen in Figures 7-8. Figure 7, as described earlier, re- 
sults from a solution in which only the coefficients of the 
offset and rate terms of the clock-error polynomial were esti- 
mated. Figure 8 shows the same data but with the coefficient 
of the second-order term in the clock- error polynomial also 
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estimated. Figure 13 shows the top part of Figure 8 with the 
vertical scale e.xpanded. In Figure 13 it can be seen that in 
addition to the variation in the clock rate there were two 
discontinuities in the clock offset, each with a magnitude of 
about 10 ns. Not only were there discontinuities in the clock 
offset, but a careful inspection of the lines of residuals will 
show that the two separate lines overlap in the region between 
t = 1.083 and t = 1.16G. The only reasonable explanation for 
this behavior is that the clock offset must have "jumped" back 
and for»_h between two values separated by about 10 ns (Refer- 
ence 12) . 

There are two ways of handling this type of pathological 
clock behavior with the present data-processing program (VLBI3) . 
The first way is to try to fit the curves in Figure 8 with poly- 
nomials, and the second is to form differenced observations, or 
differenced and summed observations, as discussed in Chapter 2. 

Figures 14-20 show some of the stages of the process of 
determining the "best" set of polynomials to use to fit the 
clock errors of a given set of data. Figure 14 shows the Hay- 
stack-NRAO delay residuals front the observations of October 
1973, from a solution in which the coefficients of the first 
three terms of the clock -error polynomial were estimated and 
all the baseline coordinates, source coordinates and atmosphere 
parameters were fixed at a priori values. There are two fairly 
clear discontinuities, one each near t = 3.25 and t = 3.75 on 
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the vertical scale. Since there are three distinct extrema in 


the first section of the curve, it looks as though a fifth-or- 
der clock-error polynomial is called for in this section. The 
second section of the curve is so short that a first-order poly 
nomial (straight-line) will probably be adequate. The last 
section of the curve shows some curvature and will require at 
least a second-order polynomial. Figure 15 shows the residuals 
after fittiiig these polynomials to the data. The magnitude of 
the variation is reduced from 40 ns to about 10 ns, but it is 
obvious that some higher frequency noise is present. For the 
next stage, I tried fitting coefficients of a third-order clock 
error polynomial between each extremum of the first section of 
the curve in Figure 15. The last section of the curve presents 
some special problems, however. It is shown with an expanded 
vertical scale in Figure 16 It is obvious that there is a dis 
continuity of some sort near t = 4.375 on the vertical scale, 
but where exactly is the break? The curves overlap, just as 
they did in Figure 11 The simplest way to deal with this prob- 
. is to eliminate all of the points in the bottom right-hand 
s j • the curve. The remaining points show enough curvature 
, jably require at least a second-order polynomial. The 

j •'. luals from this solution are shov;n in Figure 17. The resi- 
duals are now at a level of about 3 ns and are still quite sys- 
tematic; however, at this level it is very difficult to sort 
out clock effects from the effects of errors in the a priori 
baseline and source-coordinate parameters. These errors 
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as Described In Section 4.4.5 


represent the limiting factor in our ability to determine clock 
effects by this method. 

In Figures 18-20 we see the residuals from a similar set 
of computer runs set up to determine a useful parameterization 
of the clock errors for data obtained in October, 1973, from the 
Haystack-Goldstone baseline. In Figure 18 are shown the resi- 
duals from a solution in which only the coefficients of the off- 
set and rate terms in the clock-error polynomial have been es- 
timated. Since there were two extrema on the residual p. 
the coefficients of a fourth-order polynomial were estimated, and 
the resulting graph of residuals is shown in Figure 19. The am- 
plitude of the extrema decreased by less than a factor of two. 
Next I tried to estimate the coefficients of a separate second- 
order polynomial in each interval between the cusps (at t= 2.5, 
t « 2.75 and t =: 2.85 on the vertical scale) which appear in 
Figure 19. The resulting residuals are shown in Figure 20. Once 
again, I approached the point where clock errors are hard to 
separate from a priori coordinate errors. 

There are many reasons why the above procedure is unsatis- 
factory, not the least of which is that it is very tedious and 
requires large amounts of computer time. More f undcimentally 
this technique is too dependent upon the all-too-fallible judg- 
ment of the operator. The number of possible parameterizations 
of the clock is too large, and thj criteria for selecting one 
over another are too vague. This technique works reasonably well 
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when the clocks are well behaved, that is, when there is at most 
a small number of very distinct problems such as the clock-rate 
discontinuity shown in Figure 5. However, for badly behaved 
clocks the inadequacies of this technique led to the use of a 
scheme involving differenced observations as discussed in Chap- 
ter 2. 

With differenced observables, the sensitivity of the data 
to clock errors drops dramatically. As noted in Chapter 2, if 
the observations to be differenced are exactly simultaneous, the 
clock errors disappear completely from the differenced observa- 
tions. If the observations are not exactly simultaneous then 
the degree to which the clock errors are eliminated depends on 
the spectrum of the fluctuations in the clocks. That is to say, 
if the periods of the main Fourier contributions to the power in 
the fluctuations are short compared to the intervals over which 
differences are formed then there will be little advantage in 
using differer. -ed data. Tf, however, the periods of these fluc- 
tuations are long compared to the differencing intervals, then 
the fluctuations will be largely eliminated by forming dif- 
ferenced observations, or other suitable linear combinations. 

With differenced observations the problem of selecting the 
pairs of observations to be differenced replaces the problem of 
selecting among clock-polynomial possibilities. There is an 

5 

astronomically large number (-2 x 10 ) of ways to combine the 
600-odd observations in Figure 8 by pairs. The general con- 
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straints that we wish to observe in fomiing pairs for differ- 
encing are that the observations be closely spaced in time, and 
that they not be observations on the same source. (Forming dif- 
ferences of closely spaced observations of the same source loses 
almost all the geodetic and astrometric information contained 
in the observations.) To determine a reasonable criterion for 
the maximum spacing in time to be used for the differenced ob- 
servations, it was noted that the clock -rate coefficients were 

-12 

typically around 2 x 10 sec/sec, so that a spacing of 30 min- 
utes would give a maximum clock error for each differenced ob- 
servation of less than five nanoseconds. Couple this result 
with the fact that the clock-rate offset can still be estimated 
with the differenced data (albeit with greatly reduced sensiti- 
vity) and it should be clear that differencing the data can re- 
duce the effect of the clock error to the level of nanoseconds. 
(Differences which span clock discontinuities such as are seen 
in Figure 13 must be weeded out, of course. Fortunately, there 
are only a very small number of such cases.) Since there were 
about ten thousand observations to be differenced in the entire 
Phase II data set, it was not practical to select the lifference 
pairs by hand, and a computer program was therefore written to 
select such pairs according to tlie following algorithm: The 

observations were scanned sequentially in time. If each of the 
first pair of observations was on the same source, the first ob- 
servation was deleted, and the next two observations were con- 
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sidered as a pair. This process was repeated until a consecutive 
pair of observations was found on two different sources. The 
time interval between the observations was then checked. If it 
was found to be greater than 30 min’’.tes, then the first observa- 
tion was again deleted, and the process of searching for a con- 
secutive pair of observations with different sources was resumed. 
If the time interval was found to be less than 30 minutes, then 
a differenced observation was formed from this pair of undif- 
ferenced observations, and the process of searching for a pair 
of observations on different sources was resumed with the obser- 
vations following this pair. This scheme is not optimal; its 
virtue is simplicity. In particular if the observations are 
grouped, for example, by triplets on a given source this scheme 
may result in deleting more observations than is strictly neces- 
sary. That is, by deviating a small amount from a strict time 
sequence it may be possible to construct many more pairs within 
the 30-minutes time constraint if there are a number of clusters 
of observations on the same source within this limit. However, 
the general problem of constructing an optimal scheme for select- 
ing differenced pairs has not been tackled. 

The results quoted in Section 4.5 for the Phase II data up 
through July 1974, are all from solutions in which the data 
from every baseline involving NRTvO or Goldstone were differenced 
according to the above scheme, except for the August 1973 data, 
in which the Goldstone clock was determined to have been well 


enough behaved that only the data involving NRAO were dif- 
ferenced . 

4.5 DISCUSSION OF THE RESULTS 
4.5.1 Introduction 

To obtain the results discussed in this section, the 
formal errors for the individual observations we^e first cal- 
culated from the signal-to-noise ratio as discussed in 
Reference 22. In order to test how closely the formal er- 
rors corresponded to the real scatter in the post fit resi- 
duals, the delay and delay-rate postfit residuals were cal- 
culated from a solution in which only the delay data were 
used to estimate the baseline, source-coordinate, atmosphere, 
and clock-error parameters. The rms value of the weighted 
residuals was then calculated: 



where r is the residual for the i^^ observation, w^ is the 
corresponding formal error, and N is the total number of ob- 
servations. If the formal errors accurately reflect the 
scatter in the residuals, tnen W should be close to unity. 

The values of W fcr the delay residuals were nearly all 
in the range of 2.5 to 3.5, and the values of W for delay- 
rate observations were found to range from a^x^ut 60 to about 
130. The reason that the delay-rate formal-error calculations 
are so at variance with the observed scatter in the dclay- 
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rate residuals is that the formal error calculations do not 
consider instrumental instabilities. In particular, the 
formal errors based on the signai-to-noise ratios have values 
on the order of 10 sec/sec, whereas the frequency standards 
used in the experiments were usually stable to no better than 
10 sec/sec (Reference 12) . 

l*n order to ensure that the relative weighting of the 
delay and delay-rate data is appropriate for the residual 
scatter actually seen in each type of data, the formal er- 
rors for each observation were multiplied by the appropriate 
value of W calculated as discussed above. The resulting 
values for the formal errors were found to produce values for 
W from combined delay and delay-rate solutions which were 
close to unity. Since the clock errors appear to dominate 
the delay-rate errors, the fo„*mal errors for the de3.."y rate 
should have been nearly uniform, rather than being based on 
the signal-to-noise values for those observat ionr . However, 
a qualitative analysis indicates that the necessary charges 
in the delay-rate v'eighting function wou3d not have a signi- 
ficant effect on the final estimates of the parameters. 
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4.5.2 Baseline Results 


The baseline results displayed in Table 4 * give the base- 
line length and the X, Y, Z components in the Earth-fixed co- 
ordinate system described in Chapter 3 ** The reasons for dis- 
playing the length in addition to the three cartesian components 
are two-fold: first, the baseline length is intrinsically in- 

teresting since it is possible to construct the vector separations 
of a net of three or more stations solely from measurements of the 
lengths of the baselines between these stations; second, because 
length is invariant under rotation, it is unaffected by errors 
in the tabulated values for polar motion and UTl. 

It can be seen from Tabic 4 that the baseline length solu- 
tions for each baseline from the different sets of observations 
have an rms scatter about the wsighted mean of less than cne meter. The 
scatter in length ranges from 79 cm for Haystack-Alaska to a mere 
7 cm for Haystack-NRAO. The scatter in the baseline cartesian 
components is considerably larger, ranging up to 324 cm for the Z 
component of the Haystack-Alaska baseline. The reasons for the 
larger scatter on the components are two-fold: first, the com- 

ponents are affected by evrors in the tabulated polar motion and 

* 

See Appendix C for a discussion of the entries in Table 4. 

** 

The metric-system lengths given in Table 4 and elsewhere in 
this thesis were based on the speed of light being 2.997925x 
10 5 km/sec. 
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UTl values; second, the baseline error ellipsoids* tend to be "pan- 
cake” shaped, with the smallest axis roughly parallel to the 
baseline — this is evidenced by the fact that the formal errors 
of some of the baseline components are larger than the formal 
errors of the baseline lengths. 

It will be noticed that the deviation from the mean in some 
cases is very much larger than the quoted formal error (c^. the 
Haystack-Goldstone length from Run Number 75125-11 in Table 4) . 

Are these results to be interpreted as wildly improbable statis- 
tical scatter? No, because the formal errors are based on an 
assumption that the noise in the data is gaussianly distributed 
with zero mean. In actual fact there are almost always systema- 
tic effects clearly visible in the postfit residuals. (These 
systematic effects are probably caused largely by clock varia- 
tions and other instrumental variations and in part by the atmos- 
phere and ionosphere.) When the deviation from the mean is very 
much larger than the formal standard error it is quite likely 
that the systematic errors predominate over the statistical 
errors. For this reason the scatter of the various results is 
likely to be a better measure of the true errors (systematic 
errors plus statistical errors) in the corresponding parameter 
estimates than are the formal errors themselves. It seems fair 
to say that, on the basis of the scatter of the various data 

* 

For a discussion of error ellipsoids, see Reference 15, pp. 51-53. 
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sets, we have determined the length of all of our baselines with 
an uncertainty of less than one meter, and we have determined 
th equatorial components of each baseline with an uncertainty of 
about one meter, and the polar component with an uncertainty of 
several meters. 

At this point, a word needs to be said about the relation- 
ship between these results and the results reported in Referen- 
ces 13 and 18, which results were based on some of these same 
data. The primary difference between these results and those 
reported earlier was brought about by a change in the weighting 
function used in processing these data. In the earlier analysis 
the values for the standard errors in the observations were based 
on the signal strength and on the elevation angles, since it was 
assumed that the atmosphere would degrade the quality of the low- 
elevation data. However, no evidence could be found to support 
the idea that the low-elevation-angle observations were of sig- 
nificantly poorer quality than the high-elevation-angle observa- 
tions. Indeed, a graph of residual magnitude vs , the root-sum- 
square value of the atmosphere corrections for the two sites shows 
no discernible dependence of residual size on atmosphere correc- 
tion. The graphs for delay and delay rate are shown in Figure 21. 
The horizontal scale on these graphs represents the root-sum- 
square value of the atmosphere correction at each site, with the 
zenith electrical path length (p^) of the atn\osphere normalized 
to unity. (To convert the horizontal scale to an approximate 
value for atmospheric correction in seconds of time for delay 
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Figure 21. 


Magnitude of Delay and Delay Rate Residuals 
as a Function of Atmospheric Correction (see 
Section 4.5.2) . 







and seconds per second for delay rjte, simply multiply each scale 

-9 

by 7 X 10 sec.) In order to eliminate the effects of variation 
in the signal-to-noise ratio between observations of different 
sources, each residual has been normalized to a dimensionless 
quantity by dividing the residual by its formal standard error 
based on the signal-to-noise ratio for that particular observa- 
tion (Reference 22). In order to facilitate inter- 

pretation of these graphs, each graph has been divided into ver- 
tical segments with two segments between each fiducial mark on 
the horizontal scale. For each segment the total number of 
points within the segment has been graphed (plain curve) as well as 
the standard deviations of those points (curve marked with boxes 
as: -0 — 0 - Q ) . It can be seen that this latter curve shows no 

pronounced trend toward either larger or smaller residual 
standard deviation as the atmosphere correction increases. This 
graph should not be interpreted to mean that the atmosphere does 
not degrade the data, but rather that the residuals may be domina- 
ted by non-atmospheric effects. Therefore, the results quoted in 
this thesis have been based on relative standard errors for the observa- 
tions derived from the signal-to-noise ratios only. This dif- 
ference caused only small changes in the final estimated para- 
meter values, but sometimes caused large changes in the formal 
errors for these parameter estimates. The other differences between 
these results and the results published earlier are all due to small correc- 
tions in the theoretical models used in VLSI 3. These corrections are five- 
fold. (The maximum magnitudes of these effects are given in parentheses.) 
First, the BIH UTl values have been used instead of the U.S. Naval Observatory 
IfTl values (*•10 msec) . Second, the effects of diurnal polar motion 


have been added (-50 cm). (The diurnal polar motion and 
the change in UTl values will not affect the estimation of 
baseline lengths, since both effects represent rotations of co- 
ordinates.) Third, the Earth-tide model has been improved 
('•S cm) . Fourth, the effects of antenna motion on delay rates 
have been added ("“1.7 ps/sec; this affects only the Alaska 
data) ; and fifth, the effects of the time derivative of 
(Al-UTl) have been included in the delay rates (0.04 ps/sec). 
These effects are described in Chapter 3. 

4.5.3 Source-Coordinate Results 

The source-coordinate solutions are listed in Table 5. 

The reference for right ascension is the ri^t ascension of 3C 273B, which 
has been fixed at 12 2b 33?236 (Reference 13) . The coordinates shewn 
in Table 5 are appropriate for solai'-systcm barycentric coordinates — 
that is, the effects of elliptic aberration are not present. The scatter 
in right ascension for those sources for which there was a reasonable 
number of observations ranges from 0V02 to 0'.'09 (note that the table en- 
tries are in time-seconds rather than arc-seconds) while the 
.latter in declinations ranges fran 0'.'02 to 0V23. The largest scatter in the 
declinations occurs for sources whose decl inations are in the range of 
+5° to -5° (3C120, 3C273, 2134-00) . This increase in scatter for the low 
declination sources is to be ejected, since the sensitivity of delay and 
delay-rate observations to the declination of the sources decreases as the 
declinations go to zero. (The reason for this decrease in sensitivity to 
the source declination is that the derivative of delay rate with respect 
to source declinations goes to zero as the declination ^es to zero. The 
reason for the decrease in sensitivity for delays is slightly more ccmpli- 
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cated-. If we resolve the delay into a sinusoidal ccnponent and a constant 
cxnpcnent [as in Equation (2-3)], then the derivative of the sinusoidal 
oorponent with respect to declination goes to zero as the declination goes 
to zero. Ihe only sensitivity to declination at zero declination fran 
either delay or delay-rate c±>servations thus comes frtm the constant com- 
ponent of delay.) Our previously published results on source coordinates 
(Reference 13) were based on seme of these sane data. Ihe results in Tab- 
le 5 agree with those earlier results for each coordinate to within less 
than 1.5 times the standard errors quo. ^d in the reference. It seems fair 
to say that, except for the declinations of the Icw-declination sources, we 
have determined the positions of these sources with a standard error of less 
than a tenth of an arc-second in all cases. 

4.5.4 Baseline-Closure Tests 

It is all very well to repeat a neasurement N tiirves and then claim 
that the scatter in the results is related to the real i^'=*rtainties in the 
measurement procedure, but actually tt-Ls scatter only places a lower bound 
on the uncertainty. If there are systematic errors which repeat frero run 
to run, the effects of these errors will not show up in the sccitter of the 
results. It is therefore desirable to have additional tests of the con- 
sistency of the data vhich are not based on mere repetition. One inpor- 
tant test involves separating the irulti-station experiment data by baseline 
and solving for each baseline separately, and then checking to see hew well 
the resulting baseline solutions form closed triangles. The results from 
this type of test are given in Table 6. The stations involved in each tri- 
angle are indicated by the same station co-es as are used in Table 4. 
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experiments 


The column labeled "Data" indicates 
the data for each baseline came from 
code; 


which experiments 
/ according to the 


following 


1 = May-June, 1972 

2 = August, 1973 

3 = October, 1973 

4 = January, 1974 

5 = March, 1974 

6 = April-May, 1974 

7 = July, 1974 

8 = January, 1975 

It will be seen that in Table 6 all the data from baselines in 
a given triangle were taken from the same experiment. 

All of the X and Y conponent closure errors for the H-N-S triangle are 
less than one meter, and the Z ccnponent , which is more poorly deter- 
mined, closes to within two meters or better. On the H-N-G tri- 
angle the X components all close to within 30 cm or less, while 
the y and Z components, which are less well determined ( cf . 

Table 4 ) close to within 2.5 meters or less. The remaining 
triangles are too degraded by the poor coverage of data on the 
G-S baseline to match these results, but nevertheless their clo- 
sure results are in fair agreement with their much larger formal 
errors. 

It has been argued that in a three-station experiment the 
delays on the third baseline are not independent of the delays 
on the remaining two baselines, and that therefore the closure 
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test as just examined does not really test the validity of 
the VLSI measurements. To sidestep this problem the baseline 
results used to compile Table 6 were re-combined into triangles 
with the requirement that no more than two of the baselines could 
be obtained from data from a single experiment. A small frac- 
tion of the hundreds of such possible coiTiOinations is shown in 
Table 7. The problem with examining closure in this manner is 
that the results have been degraded by possible errors at the 
meter level in the tabulated polar-motion and UTl information 
used in the solutions. Therefore these results cannot be expect- 
ed to close to much better than one meter or so. Nevertheless 
the results of this closure test place bounds on the total er- 
rors in the baseline parameter estimations, for on the triangles 
which have a reasonable amount of data for all baselines (H-N-S 
and H-N-G) practically all of the closure errors are less than 
two meters. These results at the very least fail to contradict 
our claim to meter level or better accuracy in the baseline 
determinations . 

4 . 5 .5 Antenna-Axis-Offset Solutions 

Another interesting test of the consistency of the VLSI 
measurements concerns the magnitude of the antenna-axis offsets 
discussed in Chapter 3. The axis offset is one of the few para- 
meters which affects the VLSI observables and can be measured 
accurately by conventional techniques to the centimeter level 
or better. If we therefore use the VLBI data to estimate this 
axis offset, and compare the results to the known values we can 
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perhaps set some bounds to the expected errors in the baseline 
solution values. This test is not perfect, as the sensitivity 
of the measurements to the axis offset will not, in general, be 
the same as the sensitivity to the baseline components. Never- 
theless, this test provides an interesting check on the entire 
VLBI procedure. The results from this type of test are shown 
in Table 8 where the data-set codes a::e the same as in Tables 6 
and 7. The measured values for the various antenna axis offsets 
are as follows*: 

Alaska 7.2898 m 

NRAO 140* 14.9273 

Onsala 2.15 

The Alaska and NRAO results seem to confirm our state- 

ments regarding meter- level accuracy of the solutions. A ques- 
tion which requires an answer, of course, is why there is so 
little sensitivity to the Sweden axis offset, especially as com- 
pared to the NRAO sensitivity, since both telescopes are equator- 
ially mounted. A possible reason for the difference is that the 


The value of the Alaska antenna-axis offset was obtained from 
a personal communication from Mr. Sal Rocci of Rohr Industries, 
Soladyne Division, San Diego, California. The value of the 
NRAO antenna-axis offset was obtained from NRAO Drawing #36D00022, 
March 28, 1973. The value of the Onsala antenna-axis offset was 
obtained from Reference 12. 
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Sweden baselines are necessarily poorer in declination cover- 
age because of the high latitude of the Sweden station, and 
poorer in hour-angle coverage because of the longitude of the 
Sweden station. 

4.5.6 Polar Motion and UTl 

It has been recognized ~ince the late 19th century 
that the Earth's rigid crust is net fixed relative to its 
pole of rotation, nor is the Earth's rate of rotation constant. 
The I L S and the B I H have monitored the location of the 
pole of rotation since the early part of the 20th century and 
they publish at 5-day intervals the Cartesian coordinates of 
the pole relative to a defined mean pole. In additxor , by 
the middle of the 20th century, clocks had been developed which 
had sufficient stability to allow determination of the irregu- 
lariti^iS in the Earth's rate of rotation. Thus, in 1956, the 
B I H began measuring and publishing the difference between 
time as kept by the spinning Earth (UTl) and time as kept by 
their atomic clocks (Al or UTC) (Reference 1) . 

We have used both the Phase I and the Phase II data sets 
to estimate both the coordinates of the pole with respect to 
an Earth-fixed coordinate system and the differences between 
Al and UTl. The differences between our estimates and the 
B I H estimates are tabulated in Table 9 . For those experi- 
ments in the Phase I data set that involved only a single 
baseline, only one component of the pole position was esti- 
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mated, along with the Al-UTl offset, since there are only two 
degrees of rotational freedom for each b?seline. In any event, 
the component of the pole position and the (Al-UTl) offset can- 
not be estimated simultaneously with the estimation of the 
orientation of the baseline. To do so would be to attempt to 
determine the orientation of the baseline in addition to its 
3 vector components. To surmount this problem, the coordinates 
of the pole and the (Al-UTl) offset were fixed at the B I H 
values for part of the data in both the Phase I and the Phase II 
data sets. The selection of the data set to be fixed was, by 
nature, somewhat arbitrary, and it should be noted that any 
errors in the B I H polar motion and (Al-UTl) values ^or the 
days selected will bias all cf the polar motion and UTl esti- 
mates. In the Phase I data set, the X component of th 2 polar 
motion and the Al-UTl offset were fixed for the August, 1972 
data, in order to >>e consistent with our previously published 
results (Reference 18 ) . In addition, the Haystack-Goldstone 
data set from January, 1975, because of its similarity to 
most cf the Phase I data sets (one baseline, one day of obser- 
vations, and well-behaved clocks at both stations) was pro- 
cessed with the Phase I data. Three of the Phase I experiments 
involved three stations, and therefore the second (Y) compo- 
nent of the polar motion could be estimated for any two of 
these experiments if it were fixed for the third. The Y com- 
ponent of the polar motion was therefore fixed for the June 6- 
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7, 1972 data. In the Phase II data set (exclusive of the 
Jcinuary, 1975 data) both components of the polar motion and 
the (Al-UTl) offset were fixed for the January, 1974 data. All 
of the Phase I and the Phase II data should have been included 
in a single computer run to estimate all of these polar‘-motion 
and UTl parameters simultaneously. Unfortunately neither of our 
data processing programs is able to handle the number of para- 
meters required for such a solution. We have therefore proces- 
sed the Phase I data in two separate solutions, as described 
below, and the Phase II data ii five separate solutions. It 
would be possible, of course, to reduce our observations to 
"normal points", that is, to condense each data set down to its 
relevant essentials and thereby streamline the computer analysis 
of the entire data set. The computer programs required to 
create such normal points have not yet been written. 

The polar-motion and (Al-UTl) values for the first 8 
experiments given in Table 9 are not identical with the values 
published in Reference 18 , One of them, the X component of 
polar motion for March, 1973, is different from the published 
value by more than twice its formal error. There are several 
differences between the way in which the data were handled 
for the results reported in this thesis and the way in which they 
were handled for results reported in Reference 18, in addition 
to the differences noted in Section 4.4.4. First, all of the 
Phase I data (except the January, 1975 data which had not yet 
been obtained) wore included in a single solution for the results 
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rqxirted in Reference 16. For the results reported in this thesis, the 
data before August 1972 were included in one solution, and the Phase I 
data after August 1972, together with the Haystack-Goldstone data from 
January 1975, were included in a second solution. Second, because of the 
problems noted in Section 4.4.4, the data from June 27-28, 1972, were not 
included. All of these changes affect the way in vhich the systematic er- 
rors present in these data affect the solutions. The fact that three of 
the results changed by about twice their formal errors suggests that the 
"real" uncertainties in many of these results might be at least of the 
order of twice the formal errors (see, also, p. 112) . The mean and stan- 
dard deviation of the differences between our polar motion and UTl esti- 
mates and the corresponding BIH measurements is given in Table 10. The 
corresponding values from Reference 18 are also given. 

The postfit residuals froiri the solutions discussed in this section have 
a semevhat different character than the residuals from the solutions dis- 
cussed in Sections 4 . 5.2 and 4 . 5 . 3. A "typical" set of residuals from a 
polar-motion and (Al-UTl) solution :s shcn^n in Figure 22. The residuals 
shown are the August 29-30, 1972 delay residuals from the second Phase I 
solution discussed above. The residuals are much more systematic than the 
residuals from solutions which include only the August 29-30, 1972 data 
( cf . Figure 11 ) . These systematic trends are apparently caused by varia- 
tions in the instrumental delay since the trends seem to be cenmon to all 
of the data even though the data include i her leaved observations on many 
different sources. Tliese trends are apparently being masked in the 
single-day solutions by adjustments of the baseline - 
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Figure 22. Delay Residuals from a Polar Motion and UTl Solution 


coordinate and source-coordinate parameters. Figure 23 sliows 
residuals from a solution involving the same data set as em- 
ployed in iigure 22 but with the coefficients of a third- 
order clock polynomial being estimated for each day's data. The 
systematic trends shown in Figure 22 are greatly reduced. Two 
of the eight polar motion and (Al-UTl) values estimated in the solution 
using third-order clock polynomials differed from the results 
shown in Table 9 by nearly five times the formal errors for 
those parameters. More work is needed to investigate fully the 
effects of systematic errors on the polar -motion and Al-UTl 
results . 

4.5.7 The Precession Constant 

The action of the gravity fields of the Sun and Moon on 
the Earth's equatorial bulge produces a motion of the Earth's 
spin axis which is referred to as the luni-solar precession 
(Reference 22, pp. 181-184) . In ecliptic coordinates the 
motion of the spin axis can be described by two parameters. In 
effect, the first parameter is t±»e "mean" angle between the spin axis and 
the ecliptic pole, and is referred to as the mean obliquity. (See Reference 
24 for precise definitions.) The second parameter, in. effect the rate at 
vtiich Lhe spin axis rotates about the ecliptic pole, is referred to as the 
precession constant. In Reference 1, the values for these constants at 
the 1950.0 epoch are given as: 

Mean Obliquity = 23® 26' 44 '.'84 
Precession Constant = P = 5026 '.'75/century 
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In other words, the spin ‘axis completes one revolution about 
the ecliptic pole in ~26,000 years. Since VLBI observations 
have virtually no sensitivity to the location of the ecliptic 
plane, we. do not expect to measure the meah obliquity from 
VLBI data alone. However, we can measure that portion of 
precession which results from the motion of the Earth's pole 
in space. 

As mentioned in 

Section 4.5.6, it is not possible with our present processing 
programs to combine all of our data into a single solution. We 
have therefore selected a subset of our data to use in solving 
for the precession constant. We selected these data according 
to two criteria: first, the data should be reasonably free of 

clock problems j second, the data should span a reasonably long 
time interval. We therefore selected the Haystack-Goldstone data 
from August 1972, November 1972, February 1973, March 1973, and 
January 1975, for a single solution to estimate the precession 
constant, P. [The baseline coordinates, source coordLates, clock-error 
coefficients, ana polar-motion and (Al-UTl) parameters were also estimated.] 
Ihe result from this solution was: 

P = 5036+1 arc-sec/century 

To see whether the addition of a -jiird site would improve our sensitivity 
to the precession constant, we used the Haystack-NKAO-Goldstc»ie data from 
August 1973 and January 1975 to estimate P with the follcwinq result; 

P = 5029+2 arc- sec/century 

Ihe difference between this value and the classical value (”2 arc sec, century) 

represents cin angular error of only 0V0 3 over the 1.5 years 

■panned by tho data. For the first solution discussed 
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the angular error is about 072 over the 2.5 years spanned by 
those data. These equivalent differences in angle are coiparable in itagni- 
tude to the scatter in the source-coordinate results reported in Section 4.5.3 . 
More accurate VLBI data or a longer time span of data or both 
will be required to add significantly to our knowledge of the 
precession constant. 

4.5.7 Th e Love Numbers 

The model for the motion of the Earth's crust in response 
to the tidal generating potential of the Sun and the M oon is 
described in Equations (3-27) and (3-28). Since the positions 

of the Sun and Moon are quite well determined, we can calculate 
the disturbing potential, u, to much higher accuracy than is 
required for Earth tide calculations. What is less well deter- 
mined is the response of the Earth to this potential, as expres- 
sed by the Love numbers I and h. We would, therefore, like to 
estimate these numbers from the VLBI observations. 

Now the magnitude of the vertical displacement of the Earth's 
crust [see Equation (3-27) ] is about 30 an, whereas the lateral displacements 
are about an order of magnitude smaller. Since 30 cm is 
about the size of the scatter in most of our baseline results 
it is clear that we will have little chance of obtaining a useful 
estimate of the lateral tidal displacements. However, we did anti- 
cipate obtaining a useful estimate of h. We obtained the following esti- 
mates from different sets of data: 
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h = .54±.02 


h = ,52±.02 
h «= .66±.03 


where the first value is from a data set containing the same 
data as described in the first solution in Section 4.5.7, the 
second value is from a data set containing all of the Haystack, 
Goldstone, and Alaska data from April 1972 through August 1972, 
and the third value is from a data set containing Haysiack, 
NRAO, and Goldstone data from August 1973 and January 1975. 

In all solutions, the baseline coordinates, source coordinates, 
clock-error coefficients, and polar-motion and (Al-UTl) para- 
nieters were also estimated. The total discrepancy between the 
first and third values listed above represents a change in 
station displacement of about 6 cm. It seems clear that with 
markedly improved VLBl data (the prospects for which will 
be described in Chapter 5) VLBI will become a sensitive test 
for measuring Earth tide effects. In fact, it can be shown 
that by exploiting the fact that the phase of the tidal per- 
turbation is often well known, it should be possible to esti- 
mate a separate value for h at each station, and thereby deter- 
mine the inhomogeneities of the Earth's response to the tidal 
potential. The programming required to implement this pos- 
sibility has not yet been undertaken. 
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CHAPTER 5 


CONCLUSIONS AND PROGNOSIS FOR FURTHER PROGPJ:SS 


In Chapter 4 we saw that VLBI observations have been used 
to estimate baseline lengths with uncertainties of less than one 
meter, and to estimate source coordinates with uncertainties leas 
than a tenth of an arc-second. In addition, these data were used 
to estimate polar motion, Al-UTl, the precession constant, and 
the Love number, h. Still, there is a great deal more work that 
needs to be done with these same data. Much of this work involves 
improved methods of dealing with very badly behaved system clocks, 
including such things as optimizing schemes for forming differen- 
ced observations and for forming more general linear combinations 
of observations. In addition, one would want to experiment with 
the use of Kalman filters and stochastic models for 
clock errors. Other things which could be done to improve our 
treatment of these data include improvements in our data proces- 
sing programs which would allov/, for instance, more parameters to 
be estimated simultaneously. (The present programs are restrict- 
ed to a maximum of about 120 parameters.) 

We also anticipate improvements in the equipment used to 
make VLBI observations, improvements which should eliminate many 
of the problems discussed in Chapter 4. The most obvious pro- 
blem is the one caused by badly behaved system clocks. Since we 

have seen that hydrogen maser clocks can be made to operate close 

14 

to their purported stability limit of one part in 1C ( cf . 
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Figure 11), the obvious solution to the clock stability problem 
is to insure that only such clocks are used as meet these 
specifications. Another method of improving clock behavior 
would involve the use of two or more clocks at each station. 

Each clock could be used to monitor the performance of the 
other clock (s). This method is likely to prove to be prohibitively 
expensive, unfortunately. As discussed in Chapter 2, clock prob- 
lems can also be eliminated by using suitable linear combinations 
of observations. 

In the future it should be possible to eliminate the ionos- 
pheric effect on the VLBI observations by using receivers capable 
of making observations simultaneously at two frequencies, say, 
at S-band and at X-band. The neutral atmosphere, being non-dispersive, 
presents a more diffii'nalt problem. One method of '"dealing with the neutral 
atmosphere depends or. resolving cne atmospheric effect into twD ccrnponents, 
one due to the watel^- vapor content of the atmosfhere (the ”wet" component) , 
and the other due to the rest of the atmosphere (the "dry" component) . The 
advantage of this scheme is that the dry component is much more uniform 
spatially and temporally than the wet component. The dry com- 
ponent could therefore be modeled in the same manner as was des- 
cribed in Chapter 3 for the total atmospheric effect, w’hile the 
wet component could be monitored through observations of the 
brightness temperature of the atmosphere near and at the 22 GHz 
water-vapor absorption band. 

Another anticipated improvement in the VLBI equipment in 
the near future is the development of wide band (~50 MHz) tape 
recorder systems. (The tape recorders employed in the experi- 
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merits reported in this thesis are limited to a maximum band- 
width of 360 KHz.) This additional bandwidth could not only 
make possible the elimination of the group-delay ambiguity, as 
mentioned in Chapter 4 , but will also allow the use of very 
much smaller antennas to achieve the S6one signal-to-noise ratio 
as is attainable with larger antennas employing the smaller 
bandwidth. The ability to employ smaller antennas is important, 
since the number of antennas the size of Goldstone or even Hay- 
stack is small, and observing rimo on such antennas is severely 
limited. There is, on the other hand, a substantial number of 
smaller antennas (10m to 25m in diameter) around the world. Also, 
the ability to use smaller antennas opens up the possibility of 
using small transportable antennas, thereby allowing VLBI stations 
to be set up nearly anywhere on the land surface of the Earth, 
as proposed in Reference 17. Recently, some experiments have been 
performed with a relatively large (=10 m diameter) transportable 
antenna (Reference 11); however, use of the wider-band recording 
system will enable at least a threefold reduction in size while 
still yielding higher system sensitivity. 

In addition to increasing the instantaneous, or recorded 
bandwidth, it should be possible to increase the synthesized 
bandwidth by employing wider band amplifiers, thereby reducing 
the uncertainty in the group-delay measurements. 

Finally there is the problem of optimizing our observation 
schedules for the purpose of estimating baseline components 
(or whatever other parameters we wish to estimate from VLBI 
data) . It should be possible to employ linear programming 
analysis to improve, if not optimize, our observation schedules. 
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TABLE 1 


Station Parameters 



HAYSTACK 

GOLDSTONE 


ALASKA 

Location 

Westford, Mass 

Mojave, Calif. 


Gilmore Creek, Alaska 

Organization 

Northeast Radio 
Observatory Corp 

Jet Propulsion 
Laboratory 


Nati onal Oceanic 
Atmospheric Adm. 

Antenna Diameter 

120' 

210' 


85' 

Configuration 

Cassegrain 

Cassegrain 


Prime Focus 

Mount 

AZ/EL with 
intersecting axes 

AZ/EL with 
intersecting axes 

K/Y with Y 
axis offset 287" 

Locations Initially Assumed 
(Geodetic, \ifest Longitudes) 

LON = 71 29 19.201 
LAT = 42 37 23.00 
ELEV = 145 m 

LON = 116 53 
LAT = 35 25 
ELEV = 1030. 

19.150 
33.446 
85 m 

LON = 147 29 41.955 
IjAT = 64 58 42.667 
ELEV = 322 m 

Low Noise Amplifier 

Cooled Paramp 

Maser 


Uncooled Paramp 

Approximate System 
Temperature 

70®K 

30°K 


300®K 

Number of Tape Drives 
Used During Most 
Experiments 

2 

2 


1 

Antenna Slew Rate 

2 ®/sec 

w 0.2®/sec 


2®/sec 
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Location 

Organization 

Antenna Di 2 iineter 
Conf iguration 
Mount 

Locations Initially Assumed 
(Geodetic, West Longitudes) 

Low Noise Amplifier 

Approximate System 
Temperature 

Nut 3r of Tape Drives 
Us a During Most 
Experiments 


TABLE 1 (continued ) 
Station Paraimeters 


SWEDEN 

Onsala , Sweden 

Gialmers University 
of Technology 

84* 

Cassegrain 

Equatorial with DEC 
axis offset 85** 


LON =348 4 47.2 

LAT = 57 23 36.10 

ELEV = 14 m 

Maser 

50°K 


1 


0.25®/sec 


NRAO 

Green Bank, W. Virginia 

National Radio As- 
tronomy Observatory 

140.* 

Prime focus 

Equatorial, with DEC 
axis offset 48*11769 

LON = 79 50 09.47 

LAT = 38 26 15.54 

ELEV = 816.9 m 

Cooled Paramp 

120 ® 

1 


Antenna Slew Rate 


0.25®/sec 


TABLE 2 


Description of Experiments 


DAY 

STATIONS 

SYNTHESIZED 

GROUP- DELAY 


(See Ap- 

BANDWIDTH 

AMBIGUITY 

1972 

pendix C) 

(MHz) 

(ysec) 

Apr U-15 

H-G 

39 

1/3 

May 9-10 

H-G-A 

39-69-39 

1/3 

May 29-30 

H-G-A 

39-69-39 

1/3 

June 3-4 

H-A 

69 

1/3 

June 6-7 

H-G-A 

39-69-39 

1/3 

June 27-28 

H-G 

39 

1/3 

Aug 29- 

H-G 

23 

1 

Nov 7 

H-G 

23 

1 

1973 

Feb 4-5 

H-G 

23 

1 

Mar 30-31 

H-G 

46 

1/2 

Apr 11-13 

H-S 

23 

1 

May 18-23 

H-S 

23 

1 

Aug 10-14 

H-G-N-S 

* 

1/3 

Cct 12-16 

H-G-N-S 

* 

1/6 

1974 

Jan 22-23 

H-G-S 

* 

1/2 

Mar 4-8 

H-G-N-S 

* 

1/2 

Apr 2 9 -May 3 

H-G-N-S 

* 

1/2 

July 8-12 
1975 

H-G-N-S 

* 

1/2 

Jan 14 -17 

H-G-N 

* 

1/2 


* 

For these experiments the synthesized bandwidth was changed 
during the experiment. Full details are available in the 
Haystack archives. 


TABLE 3 


Auxiliary Information on Experiments 


Nominal Center Frequency 

Polarization 

Recorder Bandwidth 

No. Tracks 

Density 

Nominal Record Length 
Nominal Tape Length 
Frequency Switching Period 
Maximum No. Frequencies 


7850 MHz 

Left Hand Circular (IEEE def.) 
720 K Bits/sec or 360 KHz 
6 Data plus 1 parity - 7 
800 characters per inch 
0.2 sec 
3 minutes 
0 . 2 sec 
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RUN 




BASEUIN't: 

VALUE 

DIFF 

sigma 




M 

D 

Y 



CK 

CM 

CM 


75126- 

21 

4 

11 

73 

HS 

U b N G T H 

559938596,0938 

•'164,7344 

87,7759 

* 

753 26- 

18 

5 

18 

73 

HS 

length 

559938733,3906 

-29 ,4375 

59,0161 


75126- 

3 

3 

in 

73 

HS 

Lfc^'GTH 

559938792,7656 

29,9375 

29,8701 

* 

75326- 

6 

10 

12 

73 

HS 

L b G T H 

559936747,1406 

•15,6875 

49,6505 

* 

75126- 

8 

1 

22 

7 4 

HS 

UbNGTH 

55993884'<,7300 

80,9219 

109.3671 

* 

753 30 - 

15 

3 

4 

74 

HS 

L b (>J G T H 

559938760 , 7344 

-2,0938 

27,0566 

* 

75126- 

10 

4 

29 

74 

HS 

LbNGTH 

559938745 , 953l 

- l6 , 875 0 

20 1 9536 

* 

75126- 

l<i 

7 

8 

74 

HS 

L b iv G T H 

559938662,7656 

99,9375 

51,8065 

* 






H3 

LbNGTH 

559938762,8281 


40 .2671 


RUN 

DATt 

BASEL INt 

VALUE 

DIFF 

SIGMA 




M 

D 

Y 



CK 

CM 

CM 


75126- 

21 

4 

11 

7 3 

H S 

X - C 0 M P 

187856081 , 8984 

-18,0900 

136,8134 


75126- 

18 

5 

18 

73 

HS 

X-COHP 

187856176,3672 

76,3794 

124,0570 

* 

75126- 

3 

3 

10 

73 

HS 

X "Gf)Mp 

187356024,5000 

•75 , 4379 

38,8114 


75126- 

6 

m 

12 

?: 

HS 

X-CQMP 

187855884,4609 

-215,5274 

93 . 6595 

* 

75126- 

a 

1 

22 

74 

MS 

X-COHP 

187856177,9453 

77,9597 

70,6074 

* 

75130- 

19 

3 

4 

74 

HS 

X-CQMP 

187856177,6250 

77,6306 

44.8930 

* 

75126- 

10 

4 

29 

74 

HS 

X-CCHP 

107856130,5469 

30,5596 

40,3455 

* 

75126- 

1^ 

7 

8 

74 

HS 

X-COHP 

187856075,5508 

-24 , 4368 

104,3269 

* 


HS X-Cjnr 187856099, 9b8;5 


78,9006 


TABLE 4 


Baseline Results 


RUN 



DATS 


BASEUiNt 

value 

Dirr 

S I gma 



M 

D 

Y 



CM 

CM 

CM 

75126- 

21 

4 

11 

73 

HS 

Y-CQHP 

-516873261,1016 

37,5211 

75,2299 

75126- 

Itf 

5 

la 

73 

HS 

Y-CQHP 

-516073269,9844 

28,6348 

70,7155 

75126- 

3 

8 

10 

73 

HS 

Y-CQMP 

-516873349,1075 

-50.5680 

33,7073 

75126- 

6 

10 

12 

73 

HS 

y-cqhp 

-516873480, 4ii97 

-189.8079 

54,7814 

75126- 

a. 

1 

22 

74 

HS 

Y-COMP 

-516873351,5313 

•52.9045 

69,7287 

75130- 

1’ 

3 

4 

74 

HS 

Y-cqhp 

-536073242,2578 

56.3639 

30 , 3611 

75126- 

lu 

4 

29 

74 

HS 

Y-CQMP 

-516073243,2813 

55,34o2 

24.3166 

75126- 

IJ 

7 

8 

74 

HS 

Y-COhP 

-516873409^3672 

•110i7441 

47 ,8o6o 






HS 

Y-COMP 

-516873298,6250 


76,3961 

RUN 



DATE 


baseline 

VALLE 

DIFF 

sigma 



M 

D 

Y 



CM 

CM 

CM 

75126- 

21 

4 

11 

73 

HS 

Z-CQMp 

105277528,9434 

-681.3496 

271,0302 

75126- 

Iti 

5 

18 

73 

HS 

2-COMP 

105278036,3184 

•173.9744 

175,3865 

75126- 

3 

8 

10 

73 

HS 

Z-COMP 

105270234,3281 

24,0363 

83.5353 

75126- 

6 

10 

12 

73 

HS 

2-COMP 

3 05277557,8750 

-652.4170 

170 ,5898 

75126- 

a 

1 

22 

74 

HS 

Z-COMP 

11)5278220 , 1465 

9,8534 

392,2742 

75130- 

19 

3 

4 

74 

MS 

Z-COMP 

105278315 , 6836 

105 . 39i2 

103,8412 

75126- 

lU 

4 

29 

74 

HS 

Z-C0flp 

105278316,0176 

105,7257 

66,7923 

75126- 

I'J 

7 

8 

74 

HS 

Z-COMP 

11)5278220 , 0408 

9,7561 

162,5737 






HS 

2-CCMP 

105278210,2930 


206,9957 


Table 4 (continued) 


RUN 


DATE 

uaseline 

VALUE 

DIFP 

sigma 



H 

D 

Y 



CM 

CM 

CM 

75125- 

lU 

4 

14 

72 

HG 

length 

389999840,5938 

-42,0078 

12.9144 

75125- 

11 

5 

9 

72 

HG 

LENGTH 

389999815,9063 

-66,6953 

9,7694 

75125- 

2/ 

5 

29 

72 

HG 

length 

389999849,0156 

•33.5659 

11,!?728 

7‘.'125- 

1^ 

6 

3 

72 

HG 

length 

389999893,1094 

10.5078 

31.8757 

75125- 

16 

6 

27 

72 

HG 

length 

389999919,1016 

36 , 50n0 

17,8310 

75125- 

V 

C 

29 

72 

HG 

length 

389999883.3984 

f / v69 

4,9475 

75J 25- 

16 

11 

7 

72 

HG 

L c N G T H 

389999916^7500 

34 ,1^89 

7.9550 

75125- 

IV 

2 

4 

73 

HG 

LENGTH 

389999685,0078 

2,4063 

7,6529 

75125- 

20 

3 

3C 

7 3 

HG 

length 

389999894 , 5625 

ll,96n9 

5.9143 

75125- 

3 

8 

10 

73 

hG 

length 

389999881,6016 

-1,0 0 0 0 

11,9714 

75126- 

6 

10 

12 

73 

hG 

length 

389999874,9766 

-7,6250 

30.i437 

75126- 

0 

1 

22 

74 

HG 

length 

389999833 , 1719 

-49 .4297 

67,2919 

7 513:- 

iv 

3 

4 

7^ 

HG 

L E N C- T H 

389999992 , 3l25 

109,7ln9 

4 3 , 345 0 

75126- 

lu 

4 

29 

74 

HG 

length 

389999888,8047 

6.2031 

12 , 9630 

75126- 

1^ 

7 

8 

74 

HG 

length 

389999917,6953 

35,0938 

38 , 6907 

75130- 

47 

1 

14 

75 

HG 

length 

389999920,4219 

37,8203 

35,2825 






HG 

length 

389999882,6016 


26,3439 


Table 4. (continued) 


RUN 

DATfj 


BASELINE 

VALUE 

DIFF 

sigma 



M 

D 

. Y 



CM 

CM 

CM 

75125- 

10 

4 

14 

72 

HG 

X-COHP 

-384602678,6172 

-33,7871 

16,3282 

75125- 

11 

5 

9 

72 

HG 

X-COMP 

-384602547,1719 

97.6601 

C.9741 

75125- 

27 

5 

29 

72 

HG 

X-COMP 

-364602616,5391 

20,2953 

10.6412 

75125- 

1^ 

6 

7 

W 

72 

HG 

X-COMP 

-384602598, 5l56 

46,3152 

l7.47lb 

75125- 

16 

6 

27 

72 

HG 

X-CCMP 

-384602654,6250 

-9,7935 . 

9.2720 

75125- 

17 

e 

29 

72 

H G 

X-COMP 

-384602649,8672 

-5,0379 

2.8461 

75125- 

1» 

11 

7 

72 

HG 

X-COMP 

-384602662, 5781 

- 1 7 , 7 4 8 3 

3 . 25ic 

75125- 

1^ 

2 

4 

73 

HG 

X-CQMP 

-3846c265i ,7969 

•• 6 , 9 6 2 4 

6,2416 

75125- 

20 

3 

30 

73 

HG 

X-COMP 

-384602624 , 0938 

20.7354 

4.3223 

75126- 

3 

8 

10 

73 

HQ 

X-COMP 

-384602647,2813 

-2.4524 

8.9o9i 

75125- 

6 

10 

12 

73 

HG 

X-CQMP 

-384602668*3203 

-23,4920 

28.5254 

75126- 

b 

1 

22 

74 

HG 

X-COMP 

-384602598,3438 

46,4879 

49 , 7592 

75130- 

l« 

3 

4 

74 

HC 

X-COMP 

-384602750 , C7r»3 

*•105 , 2372 

39,3896 

75126- 

lU 

4 

29 

74 

HG 

X-COMP 

-384602629.7188 

15.1U09 

10.I?59o 

75126- 


7 

6 

74 

HG 

X-COMP 

-3846o2665i67l9 

•20.6417 

31.9294 

75130- 

47 

1 

14 

75 

HG 

X-COMP 

-384602651,6484 

-6,8165 

25,8854 






HQ 

X-COMP 

-384602644,8281 


23,8302 


Table 4 (continued) 


NJ 

00 


RUN 


OATP 

baseline 

VALUE 

diff 

sigma 



M 

1 

. Y 



CM 

CM 

CM 

75125- 

10 

4 

lA 

72 

HG 

y-ccjmp 

1B407543,-4482 

5,1594 

15.9550 

75125- 

U 

5 

9 

72 

HG 

Y-’COrtP 

18407563.7256 

25,4370 

41,4665 

75125- 

27 

5 

29 

72 

HG 

Y-COMP 

18407459,9810 

-78 .3078 

17 ,4608 

75125- 

1^ 

6 

3 

72 

Hr- 

Y-COMp 

184 r755l , ?45i 

12.9566 

?7,l895 

75125- 

1C 

6 

27 

72 

^ G 

V - C Q H P 

I84n7l37.36^'5 

-40 0 .9244 

2 8 , 3 9 2 4 

75125- 

1 7 

6 

29 

7? 

HG 

Y - C G M P 

lB4o763oi6lo8 

92.3224 

9,5536 

75125- 

1^^ 

13 

7 

72 

HG 

Y-CUMP 

l84n755Q,3l2o 

21.0236 

9 ,1859 

75125- 

IV 

2 

4 

73 

HG 

Y - C U H P 

184 n 7692 , 1V24 

1 5 3 , 9 0 3 9 

25,4746 

75125- 

2U 

3 

3C 

73 

HG 

Y " C C H P 

18407319,6226 

-218,6658 

13.0797 

75126- 

3 

8 

10 

73 

HG 

Y-CQMp 

l84n/735,i9o4 

196 , 901 6 

24,9336 

75126- 

4 , 

10 

12 

7 3 

HG 

y-cchp 

ia4r7--^ 3,270U 

4 , 98i5 

57,4610 

75126- 

6 

1 

22 

7 4 

HG 

Y - C C H p 

184q746o , U 65o 

-78.2036 

55,5982 

75lip- 

1^ 

3 

4 

74 

HG 

Y-COMP 

1H4c75«9,q117 

10.7229 

33,9297 

75126- 

10 

4 

29 

74 

HG 

Y-CQMP 

18407547,5723 

9,2839 

22.4395 

75126- 

1^5 

7 

6 

74 

HG 

Y-CCMP 

164c7645 , 0869 

106.7984 

67.8022 

75130- 

47 

1 

14 

75 

HG 

Y-COMP 

18407672,9839 

134,6955 

68,9151 






HG 

Y-COMP 

18407538,2686 


122.2850 


Table 4 


(continued) 


to 

so 


RUN 


DATE 

BASEHNE 

VALUE 

diff 

sigma 




M 

D 

. Y 



CM 

CM 

CM 


75125- 

10 

4 

14 

72 

HG 

Z-COMP 

-61982397,6191 

512 5760 

67,9121 


75125- 

11 

5 

9 

72 

HU 

Z-COMP 

-61983051.8623 

-141.6670 

55,9939 


75125- 

27 

5 

29 

72 

HG 

Z-CQHP 

-6i98286o;5264 

49 , 6t 90 

24,4811 


75125- 

1^ 

6 

3 

72 

HG 

Z-COMP 

-61983222,7451 

-312,5496 

101. ^352 


75125- 

16 

6 

27 

72 

HG 

Z-’COMP 

-61983161, 0684 

-250,8731 

66,3464 


75125- 

17 

8 

29 

72 

HG 

Z-COMP 

-6i9fi28l9 , 4fj92 

90.7657 

21,3948 


75125- 

1« 

11 

7 

7 2 

HG 

Z-COMP 

-61982971,5459 

-61 , 35q6 

45,9399 


751 25- 

19 

2 

4 

73 

HG 

Z-COMP 

-61982799 , 287i 

110.9081 

22,4874 


75125- 

2U 

3 

30 

73 

HG 

Z-COMP 

-61983141,9336 

-231.7387 

22,9573 


75126- 

3 

8 

10 

73 

HG 

Z-CQMP 

-61982793,0986 

117 , 0966 

44 , i472 


75126- 

6 

in 

12 

73 

HG 

Z-COMP 

-61982677,8057 

232,3699 

73,934o 


75126- 

ti 

1 

22 

74 

HG 

Z-COMP 

-61982673.7012 

36.4936 

140,6824 


75130- 

19 

3 

4 

74 

HG 

Z-COMP 

-61982907 , 2q6i 

2,9838 

49,4702 


75126- 

lU 

4 

29 

74 

HG 

Z-COMP 

-6i983on3, 0459 

-92.85n9 

3l.i?022 


75126- 

13 

7 

6 

74 

HU 

Z-CQMP 

-61982932,8662 

"22 , 67 il 

75.2037 


75130- 

47 

1 

14 

75 

HU 

Z-COMP 

-61983028,7422 

-118,5470 

83,9870 







HG 

Z-COMP 

-61982910,1953 


151,4348 


RUN 



DATE 

baseline 

VALUE 

DIFF 

sigma 




M 

D 

Y 



CM 

CM 

CM 


75126- 

3 

8 

10 

73 

HN 

length 

84512988,3516 

-9,5703 

10,0449 

♦ 

75126- 

6 

10 

12 

73 

HN 

length 

64513007,0352 

9,1133 

21,2273 

* 

75130- 

19 

3 

4 

74 

HN 

length 

64513010,8711 

12,9492 

18,4292 

* 

75126- 

10 

4 

29 

74 

HN 

length 

64532999,1016 

1,1797 

7,2016 

* 

75126- 

13 

7 

8 

74 

HN 

length 

64513003,7754 

5,8535 

20,6654 

* 






HN 

length 

64512997,9219 


6,7919 



Table 4 (continued) 


130 


RUN 


DATE 

baseline 

VALUE 

DIPT 

SIGMA 




M 

D 

Y 


CM 

CM 

CM 


75126- 

3 

0 

10 

73 

Hn X"C0MP 

-60952474,6836 

-0,3753 

6,7401 


75126- 

6 

10 

12 

73 

HN X-COMP 

-60 152490,5781 

-16,2700 

16,6035 

* 

75130- 

19 

3 

4 

74 

HN X-COMH 

-6C ;52466v7939 

7,5141 

11,8205 

* 

75126- 

10 

4 

29 

74 

HN X-CCMP 

-60952473,1943 

-3,6859 

4.5913 

* 

75126- 

1^5 

7 

8 

7 4 

HN X-COMP 

-60952448,5830 

25,7255 

11,'^032 

* 






HN X-COMP 

-60952474,3086 


9 ,0324 


RUN 



DATE 


baseline 

value 

DIFP 

sigma 




M 

D 

Y 


CM 

CM 

CM 


75126- 

3 

8 

10 

73 

HN Y-COMP 

46721613,8994 

-18,7764 

13,4023 

* 

75126- 

6 

10 

12 

73 

HN Y-COMP 

46721648,9727 

11,2974 

18,6097 

* 

75l3r- 

19 


4 

74 

HN ^'OOMP 

46721669, 3398 

31 , 6645 

9 , 0557 

* 

75126- 

lU 

4 

29 

7 4 

HN Y-COMP 

46721620 , 8369 

-16 ,8389 

5.9224 

* 

75126- 

I'i 

7 

8 

74 

HN Y-CCMP 

46721662 , 5459 

24 ,87q5 

12,2451 

* 






HN Y-COMP 

46721637,6758 


21,8964 


RUN 


date 


baseline 

VALUE 

diff 

sigma 




M 

D 

Y 


CM 

CM 

CM 


75126- 

3 

8 

10 

73 

HN Z"COMP 

-35275081,2568 

8,8221 

21,9591 

* 

75126- 

6 

10 

12 

73 

HN Z-COMP 

-35275058,7217 

31.3569 

34.4709 

* 

75130- 

19 

3 

4 

74 

HN Z-COMP 

-35275082,0264 

8 , o5i8 

30.7542 

♦ 

75126- 

10 

4 

29 

74 

HN Z-COMP 

-35275093,3760 

-8,2971 

13, 7667 

★ 

75126- 

I'J 

7 

6 

74 

HN Z-COMP 

-35275105,4990 

-l5 , 42o^ 

40.4621 

* 


HN Z-COMP -55275090,0/81 12,6023 


Table 4 (continued) 


131 


RUN 


DATE 

Baseline 

VALLE 

DIPT 

sigma 




M 

D 

Y 



CM 

CM 

CM 


75125- 

11 

5 

9 

72 

HA 

length 

503943202,9609 

95,7578 

124,0440 

* 

75125- 

27 

5 

29 

72 

HA 

length 

5D3948145,6328 

38,4297 

26.2233 

* 

75125- 

1^ 

6 

3 

72 

HA 

length 

503947958, 96q9 

: •148,2422 

48,8597 

* 






HA 

length 

503948107,2031 


76,5198 


KUN 



DATE 


baseline 

VALUE 

DIPT 

sigma 




M 

D 

V 



CM 

CM 

CM 


75125- 

11 

5 

9 

72 

HA 

X-COMP 

-377395121, *7813 

66,1558 

102,8700 

* 

75125- 

27 

5 

29 

72 

HA 

X-CQMP 

-377395183,4844 

4,4600 

20.4095 

* 

75125- 

1^ 

6 

3 

72 

HA 

X-COMP 

-377395251,7266 

•63^7826 

61,3275 

* 






HA 

X-COMP 

-377395187,9375 


23,6399 


RUN 



DATE 


baseline 

VALUE 

DIPP 

sigma 




M 

C 

Y 



CM 

CM 

CM 


75125- 

11 

5 

9 

72 

ha 

y-comp 

-300362066,0781 

66,5933 

81.9875 

* 

75125- 

27 

5 

29 

72 

HA 

Y-COHP 

-300362130,0859 

2.5822 

34,8042 

* 

75125- 


6 

3 

72 

HA 

y^cchp 

-300362166,5313 

•*33.8584 

53 , u33i 

* 






HA 

y-comp 

-300362132,6719 


28,4030 


RUN 


DATE 


baseline 

VALUE 

DIPT 

sigma 




M 

D 

Y 



CM 

CM 

CM 


75125- 

11 

5 

9 

72 

ha 

Z-COMP 

146011446,7539 

576.5602 

275.3502 

* 

75125- 

27 

5 

29 

72 

HA 

Z-COMP 

146010957,7695 

87.5781 

68 . 0187 

* 

75125- 

1^ 

6 

3 

72 

HA 

Z-COMP 

146010062,1445 

-808 , 0465 

174 ,?o99 

* 






HA 

Z-CCMP 

146010870,1914 


323,7698 



Table 4 (continued) 


132 


RUN DATt STAR 

M D Y 


75125- 

11 

5 

9 

72 

7=3C84 

75125- 

27 

5 

29 

72 

7 = 3C04 

75125- 

1^ 

6 

3 

72 

7 = 3CH4 

75126- 

21 

4 

11 

73 

7 = 3 C 6 4 

75126- 

19 

5 

18 

73 

7 = 3Cd4 

75126- 

3 

6 

10 

73 

7 = 3C84 

75126- 

6 

10 

12 

73 

7 = 3C64 

75126- 

6 

1 

22 

74 

7r3C64 

75l30- 

19 

3 

4 

74 

7=3L84 

75l26- 

10 

4 

25 

7 4 

7=3C84 

75126- 

13 

7 

8 

74 

7=3C64 

75130- 

47 

1 

14 

75 

7=3C64 






7=3C84 

RUN 

DATE 

star 



M 

D 

Y 


75125- 

'.1 

5 

9 

72 

7=3C84 

75125- 

27 

3 

29 

72 

7=3C84 

75125- 

14 

6 

3 

72 

7=3Ua4 

75126- 

21 

4 

11 

73 

7=3L84 

75126- 

16 

5 

18 

73 

7=3C84 

75126- 

3 

6 

10 

73 

7=3C84 

75126- 

6 

10 

12 

73 

7 = 3 c a 4 

75126- 

6 

1 

22 

74 

7 = 3 c a 4 

75130- 

19 

3 

4 

74 

7=3UH4 

75126- 

10 

4 

29 

74 

7=3Cd4 

75126- 

13 

7 

6 

74 

7=3Ua4 

75130- 

47 

1 

14 

75 

7 = 3CH4 






7 = 3C84 


value 



M M 

SEC 

RA 

3 .16 

29,54911 

RA 

3 16 

29,54745 

RA 

3 16 

29,53966 

HA 

3 16 

29,55q26 

RA 

3 16 

29,54694 

HA 

3 16 

29 , 55254 

R A 

3 16 

29,55496 

RA 

3 16 

29,54985 

RA 

3 16 

29 ,54646 

RA 

3 16 

29,54032 

RA 

3 16 

29.5513/ 

RA 

3 16 

29,54604 

RA 

3 16 

29,54961 


VALUE 


D M 

SEC 

DEC 

41 19 

51,60653 

DEC 

41 19 

51,71336 

DEC 

41 19 

51,69107 

DEC 

41 15 

5l , 71736 

DEC 

41 19 

5l , 71706 

DEC 

41 19 

51,74206 

DEC 

41 19 

51, 76933 

DEC 

41 19 

5l . 8196U 

DEC 

41 19 

5l , 7704U 

DEC 

41 19 

51, 71965 

DEC 

41 19 

51,75969 

DEC 

41 19 

51,78749 

DEC 

41 19 

51, 69366 


OIFF 

SEC 

-0 , OOOA9 
-0.00216 
-0.00994 
. 0 0066 
-0.00267 
. 00293 
. 00534 
. 00024 
-n . 00314 
-0.00129 
.00176 
-0,00157 


diff 

SEC 

-0.0S713 
.01973 
-0,00259 
,02371 
,02343 
, 0464 0 
,0^568 
,12595 
. 0 7674 
. 0i^600 
.06603 
. 09363 


SIGMA 

SEC 

,00121 
, 00096 
1 0 0973 
I 0 U3u7 
, 0 0267 
, 0 0 0 73 
, 0 0 231 
,00199 
,00117 
, 001U2 
,00255 
,00232 

,00246 

SIGMA 

SEC 

,00532 
,01463 
,10313 
,02271 
,01996 
,00760 
,01973 
, 02712 
,01067 
, 00 794 
,016/7 
,01541 

, 06067 


Table 5 


Source Coordinate Results 


RUN 


DATE 

star 



VALUE 

Dir F 

SIGMA 


M 

D Y 



H 

M 

SEC 

SEC 

SEC 

75125- 11 

5 

9 72 

9SNRA01A0 

KA 

3 

33 

22,38833 

, 00281 

,00129 

75125- lb 

11 

7 72 

93NKA0140 

RA 

3 

33 

22,38333 

-0.00219 

,00114 




9eNRAO140 

RA 

3 

33 

22,33552 


,00248 

RUN 


DATE 

STAR 



VALUE 

DIFF 

SIGMA 


M 

D Y 



0 

M 

SEC 

SEC 

SEC 

75125- 11 

5 

9 72 

9SNHA01A0 

DEC 

32 

8 

36,31305 

-0,06768 

,01509 

75125- lb 

11 

7 72 

9SNHA0140 

DEC 

32 

P, 

36,51637 

.07598 

,01599 




9=NRA01^0 

DEC 

32 

8 

36,44239 


,07171 

RUN 

DATE 

star 



VALUE 

DiFF 

SIGMA 



D Y 



H 

M 

SEC 

SEC 

SEC 

75125- lb 

11 

7 72 

10SCTA26 

RA 

3 

36 

58,93725 

-0,00000 

,00126 




1USCTA26 

RA 

3 

36 

58,93725 


. , 00000 

RUN 

DATE 

STAR 



VALUE 


DIFF 

SlCiMA 


M 

D Y 



D 

M 

SEC 

SEC 

SEC 

75125- lb 

11 

7 72 

1C=CTA26 

DEC 

• 1 

56 

16,78680 

,00000 

,10499 




lOsCT A26 

DEC 

-1 

56 

16,78680 


,00000 


Table 5 (continued) 


134 


RUN 


DATE 

star 



VALUE 




D 

Y 



h 

M 

SEC 

75125- 

11 

5 

9 

72 

11SNKA0150 

RA 

3 

55 

45 , 23q9u 

75126- 

21 

4 

11 

73 

11=NRA0150 

RA 

3 

55 

45 , 23633 

75126- 

Iti 

5 

13 

73 

ll=NHAOi5o 

RA 

3 

55 

45,23^73 

753 26- 

3 

8 

10 

73 

11=NKAC150 

RA 

3 

55 

4 5 , 2 3 V 3 5 

75126- 

6 

10 

12 

73 

11=NRAC150 

RA 

3 

55 

45 , 2 3a 7 u 

75126- 

a 

1 

22 

74 

11=\KAC15U 

RA 

3 

55 

45 , 24l2u 

75130- 

1^ 

3 

4 

74 

ll=NKACl5c 

RA 

3 

55 

4 5 , 2 3 1 9 V 

75i26- 

lU 

4 

29 

74 

11 = NKAGi5[j 

RA 

3 

55 

45 , 232!?6 






li=NRA0150 

RA 

3 

55 

45 , 23526 

RUN 

DATE 

STAR 



VALUE 



M 

C 

Y 



C 

M 

SEC 

75125- 

11 

5 

9 

72 

11=\RAC15D 

DEC 

50 

49 

20 , 0 0672 

75126- 

21 

4 

11 

73 

11=\RAC150 

DEC 

50 

49 

2(J , 1 063a 

75126- 

16 

5 

18 

73 

11=\H AC15U 

DEC 

5o 

49 

20,11485 

"^5126- 

3 

a 

10 

73 

11=NRAC150 

DEC 

50 

49 

20 ,13729 

75126- 

6 

in 

12 

73 

11 = N'K AC15 0 

DEC 

50 

49 

20 , 14a62 

75326- 

a 

1 

22 

74 

11=\'RA0150 

DEC 

50 

49 

?0 , l4o 06 

75130“ 

19 

3 

4 

74 

11=NHAC15u 

DEC 

5o 

49 

2u, 12524 

75126- 

10 

4 

29 

7 4 

ll=\HA0l5u 

DEC 

5n 

49 

20»11602 






11=NRA0150 

DEC 

50 

49 

20,10785 


DlfF 

SEC 

-0 f 00436 

I 00106 

•Of 0 0054 
,00400 
,00344 
.00593 
“0 * 00327 
-0 » 00269 


DIFF 

SEC 


-0 , 09912 
-0,00146 
, 0 0701 
, 02944 
. 04 077 
. 03222 
,01739 
. f)0fll7 


SIGHA 

SEC 

|00156' * 
|0C328 # 
100312 * 
,00108 ♦ 
,00326 * 
,00208 * 
fOUl34 * 
,00128 * 

,003/6 

SIGMA 

SEC 

,01031 * 
,012i>l * 
,01262 * 
,00701 « 

, 03331 
, 01934 ♦ 

, 0 U8V5 ♦ 
t02 Q47 * 

, 04262 


Table 5 (continued) 


RUN 


DATE 

star 



VALUE 

DIFF 

SIGMA 



M 

u 

T 



H 

M 

SEC 

SEC 

SEC 

75125- 

1^ 

6 

3 

72 

14=3C120 

RA 

4 

30 

31,59193 

, 0 0617 

t 00123 

75125- 

17 

8 

29 

72 

14=3C120 

RA 

4 

30 

31,5080/ 

. 0 0311 

t 0 0 026 

75125- 

Iti 

11 

7 

72 

14=3 U 120 

RA 

4 

3o 

3l , 50264 

0 • 0 0 3l3 

t 0 0 024 

75125- 

1^ 

2 

4 

73 

14=3Ci20 

RA 

4 

30 

31,59365 

• 0 0709 

I 0 0 1 0 9 

75125- 

20 • 

3 

30 

73 

14=3C12U 

RA 

4 

30 

31,50508 

*00011 

j 0 0 0 ^ 9 

75126- 

21 

4 

11 

73 

14=3C12U 

RA 

4 

3o 

31,5822u 

"0*00356 

, 0 0 282 

75126- 

Iti 

5 

18 

7 3 

14=3C12o 

RA 

4 

30 

31.5026/ 

-0.00309 

f 0 0 236 

75126- 

3 

8 

1 n 

73 

1 4 3 C 1 2 U 

HA 

4 

30 

31,58659 

.00082 

, 0 0 088 

75126- 

t 

10 

12 

73 

14 = 3C;i20 

KA 

4 

30 

31.50729 

*00153 

, 0 0215 

75126- 

b 

1 

22 

74 

14=3U120 

RA 

4 

30 

31.5047/ 

-0.00099 

,00197 

75130- 

IV 

3 

4 

74 

1453C120 

RA 

4 

30 

31.50095 

-0.00402 

, 00114 

75326" 

lU 

4 

29 

74 

14=3Cl20 

RA 

4 

30 

3l . 56606 

*00030 

1 0 0 0 93 

75126- 

I'J 

7 

a 

74 

14=3Ci20 

RA 

4 

30 

3l , 58300 

-0*00196 

1 0 0251 

75130- 

4 7 

1 

14 

75 

14=3C12Q 

R A 

4 

30 

31,50355 

-0,00221 

, 00231 






14S3C120 

RA 

4 

30 

31,50577 


|00313 


Table 5 (continued) 


136 


run DATfc star 

M D T 


75125- 

14 

6 

3 

72 

14=3C120 

75125. 

17 

8 

29 

72 

14=3^120 

751 :5- 

la 

11 

7 

72 

14=3C’120 

75125- 

19 

2 

4 

73 

14=3C120 

75125- 

2J 

3 

30 

73 

14=3C12U 

75126- 

21 

4 

11 

73 

14=3C12U 

75126- 

18 

5 

18 

73 

14=3 C12U 

75126- 

3 

0 

10 

73 

1 4 = 3 C 1 2 U 

75126- 

6 

in 

12 

/3 

1 4 = 3 C 1 2 U 

75126- 

0 

1 

22 

74 

14-3C120 

75l3n- 

19 

3 

4 

74 

14=3 C 120 

75126- 

1 0 

4 

29 

74 

14=3Ci20 

75126- 

U 

7 

e 

74 

14=3^120 

751 30- 

4/ 

1 

14 

75 

14=3C120 




VALUE 

diff 


D 

« 

SEC 

SEC 

DEC 

5 

14 

59 , 24180 

-0,24371 

DEC 

5 

14 

59 , 25383 

-0,23168 

DEC 

5 

14 

59 , 66012 

,17460 

DEC 

5 

14 

59,41724 

-0 ,06828 

DEC 

5 

14 

59,59603 

,11051 

DEC 

5 

14 

59 ,18290 

-0,30254 

DEC 

5 

14 

59.30991 

-0,09560 

DEC 

5 

14 

59 , 5667U 

f 0O119 

DEC 

5 

14 

59,55203 

. 06651 

U E C 

5 

14 

59 , 6055 u 

, 19999 

DEC 

5 

14 

59 , 02950 

,34398 

DEC 

5 

14 

59,50523 

, n 19 7i 

DEC 

5 

14 

59,84569 

, 36035 

DEC 

5 

14 

59,05191 

,36640 


SIGMA 

SEC 

i04558 * 
,01957 
,01898 
, 046ii4 
, 0if538 
,27259 
, 19976 
, 06969 
,128^3 
, 25396 
I 09412 
,05342 
,14199 
,14071 


14S3C120 


DEC 5 14 59,46552 


, 18627 


Table 5 


(continued) 


RUN 

M 

DATE 
D Y 

star 

75125- 

10 

4 

14 

72 

2/tOJ287 

75123- 

27 

5 

29 

72 

27=0J287 

75125- 

1^ 

6 

3 

72 

27=CJ287 

75125- 

3.6 

. 6 

27 

72 

27=0J287 

75l 25- 

17 • 

8 

29 

72 

27=CJ267 

75125- 

la 

11 

7 

7 2 

27=0J267 

75125- 

l9 

2 

4 

73 

27SOJ287 

75125- 

?o 

3 

30 

73 

27=0J267 

75126- 

3 

8 

10 

73 

27=CJ287 

75126- 

6 

10 

32 

73 

27=CJ287 

75l26- 

8 

1 

22 

74 

27=0 J287 

75130- 

19 

3 

4 

74 

27=0 J287 

75126- 

10 

4 

29 

74 

27=0J287 


27S0J287 




VALUE 

DiFr 


H 

M 

SEC 

SBC 

RA 

8 

51 

57 , 24044 

,00943 

RA 

8 

51 

57.23029 

-0,00072 

RA 

8 

51 

57,2283u 

-0.00271 

R A 

8 

5l 

57 , 22651 

-0 » 00449 

RA 

8 

51 

57,23o99 

-n , C 0 0 02 

RA 

6 

51 

57,2276/ 

— 0 t 00333 

RA 

a 

51 

57 , 23837 

• 0 0 7 3 6 

HA 

8 

51 

57,22984 

-0 . C0117 

RA 

8 

51 

57,23283 

,00162 

RA 

8 

51 

57.23192 

,00092 

RA 

8 

51 

57,2324u 

,00140 

RA 

8 

51 

57,23498 

.00397 

RA 

8 

5l 

57,23206 

.00106 

RA 

8 

51 

57,2310U 



SIGMA 

SEC 

,0015c • 
,00043 ♦ 
fOOO^fl * 
1CUIO2 
, 0 0 028 ♦ 
, 0 0 076 ♦ 
, 0 0 070 » 
, ' 0041 • 
,00089 * 
,00238 * 
, 00237 * 
,00115 # 
,00q73 * 

, 0 0257 


Table 5 (continued) 


138 




RUN 


DATt 

5 1 aR 





M 

D 

Y 



C 

751Z5- 

10 

4 

1^ 

72 

27^0J2J7 

DEC 

2n 

75125- 

27 

5 

29 

72 

27=0J287 

DEC 

2 n 

75125- 

1^ 

6 

3 

72 

27=0J267 

DEC 

2n 

75125- 

16 

6 

27 

72 

27=0J287 

DEC 

2n 

75125- 

1/' 

ft 

29 

72 

27=CJ287 

UEC 

2n 

75l'^5- 

16 

11 

7 

72 

27=0 J287 

DEC 

20 

/ 5 ■ i - 

IV 

2 

4 

73 

27=0J287 

OEC 

2n 

7 5 1 2 5 - 


3 

3n 

7 3 

27=CJ287 

UEC 

2 fi 

75126- 

3 

ft 

10 

73 

27=0J287 

UEC 

20 

75126- 

6 

10 

12 

73 

27=c J?87 

DEC 

20 

75126- 

a 

1 

22 

74 

27=0J287 

DEC 

20 

75130- 

1^ 

3 

4 

74 

27=0J287 

DEC 

20 

75126- 

10 

4 

29 

74 

27=0J287 

DEC 

2n 






27=0J287 

DEC 

20 


VALUE 

diff 

S ; GMA 


M 

SEC 

SEC 

SEC 


17 

58 , 55286 

.10073 

,02745 

.01400 

if 

17 

58,4570V 

» 00497 

if 

17 

58 , 47216 

• 02o^4 

1 02213 

♦ 

17 

58 ,60403 

•15i91 

t 0^419 


17 

58 , 44 q59 

-0 • 01153 

, 0 0877 
,0396ft 
f 02 q 34 
1 0 09Vl 

• 

17 

58 , 64713 

.19501 

■» 

17 

58 ,2466V 

•0.20543 

* 

17 

58 , 4529 (j 

, n 0 0 7 8 

If 

17 

58.50186 

. 04973 

, 03904 

it 

17 

58,62187 

, x6975 

, 0 7334 
,00695 
, 043V3 
1 0l9il 

,06671 

♦ 

17 

58 , 5080 U 

• o5568 

* 

17 

58 ,59744 

.14532 

♦ 

17 

17 

58,48537 
56 , 45212 

. 0^325 

* 



Table 5 (continued) 


RUN 


DATE 
M D Y 


star 


VALUE 

H M SEC 


Ui 

VD 


75125- 

10 

4 

14 

7? 

75125- 

11 

5 

9 

72 

75125- 

27 

5 

29 

72 

75x25- 

14 

6 

3 

7? 

75125- 

16 

6 

27 

72 

75125- 

17 

a 

29 

72 

75125- 

19 

2 

4 

73 

75125- 

20 

3 

30 

73 

75126- 

21 

4 

11 

73 

75126- 

16 

5 

18 

73 

75126- 

3 

8 

10 

73 

75126- 

6 

10 

12 

73 

75126- 

8 

1 

22 

74 

75130- 

19 

3 

4 

74 

75126- 

10 

4 

29 

74 

75126- 

13 

7 

8 

74 


23=4C39,25 

RA 

9 

23 

55,31510 

26=4C39 , 25 

HA 

9 

23 

55 , 29349 

28=4U39 , 25 

RA 

9 

23 

55 , 29561 

28=4C39 , 25 

RA 

9 

23 

55,29427 

28=4C39 , 25 

RA 

9 

23 

55 , 2"'r68 

28=4C39 , 25 

RA 

9 

23 

55,2967/ 

26=4C39,25 

RA 

9 

23 

55,30317 

28=4C39 , 25 

HA 

9 

23 

55,2935/ 

25s4C39,25 

HA 

9 

23 

55,29829 

26 = 4C39 , 25 

HA 

9 

23 

55 ,29644 

23 = 4C39 , 25 

HA 

9 

23 

55,3008/ 

26=4C39,25 

RA 

9 

2 3 

55,30253 

28s4C39,25 

RA 

9 

23 

55,29990 

2BS4C39 , 25 

RA 

9 

23 

55,29831 

23=4C39,25 

RA 

9 

23 

55,29635 

26=4C39,25 

HA 

9 

23 

55 ,30065 

28a4C39,25 

HA 

9 

23 

55.29741 


DIFF 

SEC 

.01769 
-0.00393 
•0 • 0 0179 
•0.00314 
*0.004/3 
* 0 . 0 0 065 
.00576 
'0.00384 
. 000^7 
' 0.0 0097 
,00346 
.00512 
,00249 
.00090 
0.00107 
.00324 


S 1 GMA 
SEC 

,00245 
, 0 04ii7 
I 0 0 0i»4 
, 00872 
I 00131 
. 0 0 049 
,00073 

,00060 

r 00306 
, 00266 
,00069 
,00260 
,00205 
, 001^2 
I 00105 
I 00236 

, 00351 


Table 5 (continued) 


RUN 


DATE 

STAR 



VALUE 

diff 

SIGMA 



M 

L 

T 



D 


SEC 

SEC 

SEC 

75125- 

10 

4 

14 

72 

23=4C39,25 

DEC 

39 

15 

23,70697 

-0,02675 

,02147 

75125- 

11 

5 

9 

72 

28 = 4C39 ,25 

DEC 

39 

15 

23,62621 

-0 ,10751 

, 43316 

75125- 

27 

5 

29 

72 

28=4C39 .25 

DEC 

39 

15 

23 , 7285a 

-0,00514 

,00722 

75125- 

1^ 

6 

3 

7 2 

28=4C39 , 25 

DEC 

39 

15 

23,80443 

, 07071 

,11024 

75125- 

16 

6 

27 

72 

2" =4C39 , 25 

DEC 

39 

15 

23,75362 

, 01990 

,01609 

75125- 

17 

8 

29 

72 

28=4C39 , 25 

DEC 

39 

15 

23,72620 

-0 , 0 0552 

, 0 0 434 

75125- 

ly 

2 

4 

73 

28 = 4 L’ 39, 2 5 

DEC 

39 

15 

23,60366 

- 0 , 1 3 0 0 5 

,01080 

75125- 

20 

3 

30 

73 

28=4039 , 25 

DEC 

39 

15 

23 , 75C12 

f 01641 

, 0 0396 

75126- 

21 

4 

11 

7 3 

28=4C39 , 25 

DEC 

35 

15 

23,73975 

,00604 

, 02231 

75126- 

16 

5 

18 

73 

28=4039 , 25 

DEC 

39 

15 

23,79022 

, 05651 

, 01859 

75126- 

3 

8 

10 

73 

28=4039,25 

DEC 

39 

15 

23,7209/ 

-0 , 01275 

, 0 0869 

75126- 

6 

10 

12 

7 3 

28=4039, 25 

DEC 

39 

15 

23,78625 

, 0^253 

, 0 2896 

75126- 

a 

1 

22 

74 

20=4039 , 25 

DEC 

3$ 

15 

23,74623 

, 01152 

j 0 3 i86 

75130- 

1^ 

3 

4 

74 

28=4039 , 25 

DEC 

39 

15 

23,78x86 

1 0^814 

,01643 

75126- 

10 

4 

29 

74 

28=4039 , 25 

DEC 

39 

15 

23 , 7 3s 0 6 

» 0 0 0 34 

, 0 1 1 u 2 

75 126- 

I'i 

7 

a 

74 

28=4039 , 25 

DEC 

39 

15 

23 , 6li36 

, 0 7765 

i02l23 






2834039,25 

DEC 

39 

15 

23,73371 


,03326 


Table 5 (continued) 


75125 
75125 
75125 
75125 
75125 
75125 
753 25 
75125 

75125 

75126 
75126 
75126 
75126 
75126 

S: 75130- 

^ 753 26- 

753 26- 
75130- 




UATk 

ST AK 



value 

DIFF 

S I GMA 



M 

D 

Y 



U 

M 

SEC 

SEC 

SEC 

- 

10 

4 

14 

72 

33s3C273B 

UEC 

2 

19 

43,57861 

.26853 

,04318 

* 

11 

5 

9 

72 

33=3C273B 

DEC 

2 

19 

43,49699 

,18695 

,05514 

• 

2 ? 

5 

29 

72 

3353C273B 

DEC 

2 

1 ^ 

43 , i465u 

-0 , 16359 

, 03032 

• 


6 

3 

72 

33=3C273d 

UEC 

2 

19 

43,20546 

-0.10463 

, 03528 

• 

16 

6 

27 

72 

33=3C273B 

UEC 

2 

19 

43 , 28166 

-0.02842 

, 1 0 702 

• 

17 

8 

29 

72 

3 3 = 3 2 7 3 B 

UEC 

2 

19 

43,20376 

-0.10632 

, 02535 

• 

lb 

11 

7 

72 

33=3C273B 

UEC 

2 

19 

43,3027(3 

-0 . 00 739 

, 02638 

• 

19 

2 

4 

73 

33=3C273B 

UEC 

2 

19 

43,36660 

.05652 

, 04387 

• 

20 

3 

30 

73 

33=3C273B 

UEC 

2 

19 

43 , 39604 

, 0«596 

,03294 

• 

21 

4 

11 

73 

33=3C273B 

UEC 

2 

19 

42,60881 

-0 , 70128 

, 38408 

— 

lb 

5 

18 

73 

3353C273B 

UEC 

2 

19 

43,14643 

-0,16365 

, 26552 

• 

3 

6 

10 

73 

33=3C273B 

DEC 

2 

19 

43,36296 

, 0!?2e7 

, 0 7432 

• 

6 

in 

12 

73 

33 = 3C:273B 

UEC 

2 

19 

43 , 01935 

-0.29Q73 

, 077b7 


d 

1 

22 

74 

33=3^27iB 

UEC 

2 

19 

43,47625 

f 16616 

, 46834 


19 

3 

4 

74 

33 = 3L'273B 

UEC 

2 

19 

43,3635o 

. 0^341 

,06523 


lU 

4 

29 

74 

33=3C273B 

UEC 

2 

19 

43,36503 

, o5494 

,044b6 


1^ 

7 

8 

74 

33S3C273B 

UEC 

2 

19 

43,61664 

,30655 

,06923 


4 7 

1 

14 

75 

33s3U273b 

DEC 

2 

19 

43,31765 

,00757 

,17131 






33S3C273B 

DEC 

2 

19 

43,3100b 


,13214 


Table 5 (continued) 




run date star 

M D Y 


VALUE 

h M SEC 


75125- 

10 

4 

14 

75125- 

U 

5 

9 

75125- 

1^ 

6 

3 

75:25- 

l7 

8 

29 

75125- 

18 

11 

7 

75125- 


? 

4 

75125- 

ij 

3 

3o 

75126- 

3 

8 

10 

75126- 

6 

10 

12 

75130- 

1^ 

3 

4 

75126- 

18 

4 

29 

75126- 

18 

7 

8 

75130- 

4 / 

1 

14 


72 

35=3C279 

RA 

72 

35=3C279 

RA 

72 

3 5 = 3 C 2 7 9 

RA 

7 

35=3C279 

RA 

72 

35=3C279 

KA 

73 

35=3C279 

RA 

73 

35=3 C 279 

HA 

73 

35=3U27V 

R A 

73 

3 5 = 3 C 2 7 9 

RA 

/4 

35=3U279 

HA 

74 

35=3 C 279 

HA 

74 

35=3C279 

RA 

75 

35=3U279 

HA 


35=3C279 

RA 


12 

53 

35 , 03685 

12 

53 

35 ,83336 

12 

53 

35 , 83586 

12 

53 

^5 , 83 0 16 

12 

5 3 

35 , 83i84 

12 

5 3 

35 , 83491 

12 

5 3 

35 , 33483 

12 

53 

35,8344V 

12 

53 

35,82851 

12 

53 

35,83291 

12 

53 

35,83520 

12 

53 

35 , 83779 

12 

53 

35,83508 

12 

53 

35,83398 


DIFF 

SEC 

.00286 
•0 » 0 0 065 
. 00187 
"0 • 00382 
-0.00215 
.00093 
. 00084 
.00051 
'0,00546 
-0.00107 
,00122 
, 00381 
,00109 


S I GmA 
SEC 

,00082 

,00145 

f 0 0 0i?2 

,00072 
, 0 0 069 
,00129 
, 00105 
,00159 
,00109 
, 0 0062 
, 0 0 0 /8 
j 0 0 0^1 

,00158 

,00249 


Table 5 (continued) 


HUN date star 

M D Y 


75125- 

10 

4 

14 

72 

35=3C279 

DEC 

75125- 

11 

5 

9 

72 

35=3C279 

DEC 

75125- 

1^ 

6 

•» 

72 

35s3C;279 

DEC 

75125- 

17 

8 

29 

72 

35=3C279 

DEC 

75125- 

16 

XI 

7 

72 

35=3C279 

DEC 

75125- 

XV 

2 

4 

73 

35 =3U27 9 

DEC 

75125- 

20 

3 

30 

73 

35=3C27V 

DEC 

75126- 

3 

8 

10 

73 

35 = 3C;279 

DEC 

75126- 

6 

10 

12 

73 

35=3U279 

DEC 

751 50- 

lY 

3 

4 

74 

35=3C279 

DEC 

75126- 

10 

4 

29 

74 

35=3C279 

DEC 

75i26- 

IJ 

7 

3 

74 

35S3U279 

DEC 

75130- 

47 

1 

14 

75 

35=3C27V 

DEC 






3553C27V 

DEC 

HUN 

DATE 

STAR 




M 

0 

Y 



75125- 

16 

6 

27 

72 

36SQU208 

HA 

75125- 

17 

fl 

29 

72 

3630U2q6 

RA 

75125- 

16 

11 

7 

72 

36sOU2o0 

RA 

75125- 

IV 

2 

4 

73 

36soU2o0 

KA 

75125- 

20 

3 

3Q 

73 

36S002Q8 

RA 

75126- 

10 

4 

29 

74 

36=0U2n« 

KA 






36SOU208 

KA 



VALUE 

DiPF 

SIG.MA 

D 

M 

SEC 

SEC 

DEC 

-5 

31 

8,16230 

-C . 17147 

,03766 

•. 1 

31 

7,69309 

.29774 

,08543 

• 5 

31 

8.13761 

-0.14677 

,04828 

- 5 

31 

7,83385 

, 15699 

,04424 

• 5 

31 

7,83526 

.15558 

,03742 

-5 

31 

6, 06927 

-0.07843 

,07356 

-5 

31 

7,94662 

.04421 

,06351 

-5 

31 

8,0525/ 

-0.06173 

,12670 

-5 

31 

8,40817 

-0.41734 

,12300 

-5 

31 

8,00805 

-0.0X721 

,05498 « 

-5 

31 

7,98635 

.00449 

,05958 « 

-5 

31 

8 , 07090 

-0 ’ 0^006 

,0 7699 * 

-5 

31 

7 .97611 

.01472 

,05598 * 

-5 

31 

7,99084 


, 14434 


VALUE 

diff 

SIGMA 

H 

M 

SEC 

SEC 

SEC 

1^ 

4 

45,62320 

-0.00221 

,00121 ♦ 

1^ 

4 

45,62572 

. 00030 

»00056 # 

1^ 

4 

45,62302 

-0*00240 

,00552 ♦ 

1^ 

4 

45,62699 

. 0 U158 

,00076 « 

1^ 

4 

45,62453 

- 0 * 0 0 0 39 

,00064 * 

14 

4 

45,62706 

.00164 

,00308 0 

14 

4 

45,62542 


,00114 


Table 5 (continued) 


WUN DAie STaH 

M D Y 


75125- 


6 

3 

72 

36=OQ208 

PEC 

75125- 

3 6 

6 

27 

72 

36=OU20tt 

DEC 

75125- 

17 

P 

29 

72 

36=oo?na 

DEC 

75125- 

1 a 

11 

7 

72 

36=CU2oa 

DEC 

75125- 

1*^ 

2 


73 

36=CU2ntt 

DEC 

75125- 

2u 

3 

3q 

7 3 

36rog2o« 

DEC 

75126- 

lU 

< 

29 

7- 

3 5 = 0 g 2 0 a 

DEC 






36=OU20tt 

DEC 

Rur 


DATb 

31 AK 




M 

D 

Y 



75125- 

11 

5 

9 

72 

43=3C345 

HA 

75125- 

27 

5 

29 

72 

43=3 C 345 

RA 

75125- 

16 

6 

27 

72 

4j=3C345 

R A 

75125- 

3 7 

6 

29 

72 

43=3u345 

HA 

75125- 

IV 

2 

4 

73 

43=3C345 

HA 

75125- 


3 


73 

43=3C345 

HA 

75126- 

?1 

4 

11 

73 

43=3U345 

RA 

75126- 

la 

5 

le 

73 

43=3C345 

HA 

75126- 

3 

8 

10 

73 

4 3 = 3 g 3 4 5 

RA 

75120- 

6 

10 

12 

73 

43= 3^345 

HA 

75126- 

a 

1 

22 

74 

43=3C345 

HA 

75130- 

39 

3 

4 

74 

43 = 3Cs^45 

RA 

75126- 

3 U 

4 

29 

74 

43r3C345 

HA 

75126- 

1*5 

7 

6 

74 

43=3C345 

HA 

75130- 

A7 

1 

14 

75 

43 = 3C3-«t5 

HA 






43»3C345 

RA 



VALUE 

DIFF 

S I 3 >4 A 

D 

H 

SEC 

SEC 

SEC 

28 

41 

29,41693 

.00578 

0 

28 

43 

29.36272 

-0 , 04843 

,039/5 

28 

4l 

29,41719 

. 006u4 

» 02.U9 

28 

*^1 

29 , 44960 

, p3053 

. 03i57 

28 

4l 

2 9 , 5 3 5 3 n 

. l24i5 

.04289 

28 

4l 

29,38o^U 

-0 . 0*^065 

• 02l9i 

28 

41 

29 ,43541 

, 02426 

,05586 

26 

41 

29 ,41115 


,04546 


value 

difk 

SIGMA 

H 

M 

SEC 

SEC 

SEC 

16 

41 

17,63709 

,00216 

,00123 

16 

41 

17,63402 

-0.00011 

, 0 0098 

16 

4l 

17,62852 

-1) .0 0641 

,00105 

16 

4l 

17,63520 

.00033 

,00 084 

16 

41 

17,63611 

. 00117 

, 0 0 0 9 3 

1 6 

41 

1 7 , 63335 

-n . 00158 

,00082 

16 

4l 

1 7 , 6 3 6 2 ij 

.00127 

, 0 0 3 0 4 

16 

4l 

17,63621 

.00128 

. 0 02/8 

16 

41 

17,63505 

, 00091 

.00083 

16 

41 

37.64251 

,00758 

. 0 0253 

16 

41 

17,63613 

• 00120 

100201 

16 

41 

1 7 » 6 3 0 6 5 

-0.00428 

.00177 

16 

41 

17,63477 

- IJ . 0 0 0 1 6 

,00099 

16 

41 

17,63736 

. 00243 

,00243 

16 

41 

17 , 63l0i 

-0 .00392 

. 1)0 23 2 

16 

41 

17,63493 


,00202 


Table 5 (continued) 


RUN 


DATt 

STAR 



VALUE 

DiFr 

SIGMA 



H 

U 

T 



0 

H 

SEC 

SEC 

SEC 

75125- 

11 

5 

9 

72 

4333 C 345 

DEC 

39 

54 

10,92444 

-0 , 03245 

• 0U708 

75125- 

27 

5 

29 

72 

43=3C345 

DEC 

39 

54 

10.96^23 

,00533 

, n 1 CJ 5 7 

75125- 

16 

6 

27 

72 

43=3U345 

DEC 

39 

54 

10 ,63464 

“0,12226 

, C 1 8 0 6 

75125- 

3.7 

6 

29 

72 

43=3C345 

DEC 

39 

54 

10 , 95385 

-0 , 003C4 

.00646 

75125- 

15 

2 

4 

73 

43=3U345 

DEC 

39 

54 

11,02827 

, 07137 

. 01256 

75125- 

?U 

3 

30 

73 

43=3C345 

DEC 

39 

54 

10,97276 

, 01589 

, 01022 

75l26- 

21 * 

4 

11 

73 

43=3C345 

DEC 

39 

54 

10 . 94991 ; 

"0 . 0C700 

f 023Q9 

75126- 

l6 

5 

18 

73 

43=3C345 

DEC 

39 

54 

10,95641 

- U , 0 0 0 4 H 

, n 2 0 0 8 

75126- 

3 

8 

IG 

73 

43=3U345 

DEC 

39 

54 

1 0 . 9 5 1 U 

. 0 0821 

. 0 07 /0 

75126- 

6 

10 

12 

73 

43=3C345 

DEC 

39 

54 

10,97014 

, 01324 

, 02765 

75126- 

ti 

1 

22 

74 

43=3C345 

DEC 

39 

54 

11 , 022^>4 

.07564 

. 0 3 1 4 1 

75i3o- 

l9 

3 

4 

74 

^3=3C345 

DEC 

39 

54 

lO . H790O 

- 0 • 0 ^ 7 8 3 

• *7 .A* 

. 0 2 n 1 3 

75126- 

10 

4 

29 

74 

43 = 3C;345 

DEC 

v39 

54 

10 . 90636 

, q2949 

, 0 J 94 g 

75126- 

1*1 

7 

8 

74 

43=3U345 

DEC 

39 

54 

11, 

, 08l52 

. 02i47 

75130“ 

47 

1 

14 

75 

43=3C345 

DEC 

39 

54 

10.93332 

“0.02358 

,01768 






43=39345 

DEC 

39 

54 

10.95690 


. 036V9 

RUN 

DATt 

star 



VALUE 


Dirr 

SIGMA 



M 

D 

Y 



H 

M 

SEC 

SEC 

SEC 

75125- 

16 

11 

7 

72 

5es394l8 

KA 

20 

37 

7,49006 

.00024 

,00209 

75125- 

iv 

2 

4 

73 

56=3C4ia 

RA 

20 

37 

7 ,49622 

, 00638 

. n U 133 

75125- 

?U 

3 

30 

73 

56=3C4itt 

RA 

2n 

37 

7,40277 

-1! . 0 0 7u7 

, 0 0 n 9 0 

75126- 

21 

4 

11 

7 3 

56=3C4itt 

RA 

2n 

37 

7 ,49644 

,00859 

,00357 

75126- 

It 

5 

16 

73 

56=3U4i8 

RA 

20 

37 

7,49535 

» 0 055l 

.00321 

75126- 

3 

8 

10 

73 

56s3C4itt 

RA 

2n 

37 

7 , 4«J606 

. 00622 

,00126 






56«3C418 

RA 

20 

37 

7,48984 


,00645 


> 


Table 5 (continued) 


HUN DATE star 

M D Y 


75125- 

le 

11 

7 

72 

56s3C4l8 

DEC 

75125- 

iv 

2 

4 

73 

56=3C4i8 

DEC 

75125- 

20 

3 

30 

73 

56=3C4l8 

DEC 

75126- 

21 

4 

11 

73 

56=3C4lb 

DEC 

75126- 

le 

5 

le 

73 

56=3 U 418 

DEC 

75126- 

3 

8 

10 

73 

5 6 = 3 C 4 1 8 

DEC 






56=3U4lb 

DEC 

RUN 

DATE 

STaH 




M 

D 

V 



75125- 

10 

4 

14 

T r* 

/ C 

59r2134^no 

RA 

75125- 

27 

5 

29 

72 

59=2134+00 

KA 

75125- 


6 

27 

72 

59=2134+00 

R 4 

75:25- 

1 7 

8 

29 

72 

55=2134+00 

HA 

75125- 

10 

11 

7 

72 

59=2134+00 

RA 

75125- 

l9 

? 

4 

73 

59=2134+00 

H 4 

75125- 

2u 

3 

3C 

73 

55=2154+00 

KA 

75126- 

21 

4 

11 

73 

59=2134+00 

HA 

75126- 


5 

18 

73 

55=2134+00 

HA 

75126- 


8 

10 

73 

5 = 2134 + 00 

RA 

75126- 

6 

10 

12 

73 

59=2134+00 

HA 

75126- 

b 

1 

22 

/4 

59 = 213 4 + 0 ' 

HA 

75130- 

19 

3 

4 

74 

59=2134+00 

HA 

75126- 

10 

4 

29 

74 

59=2134+00 

HA 






59e2134^tl0 

RA 



VALUE 

diff 

SIGMA 


C 

M 

SEC 

SEC 

SEC 


51 

0 

35,63900 

-n . 01832 

.09962 

* 

51 

8 

35,67990 

.02258 

,ooa'7o 


51 

8 

35,65238 

-n , 00494 

,01104 


51 

8 

35.62777 

•0,02955 

. (12429 


51 

8 

35,50177 

-0.06555 

, 02125 


51 

8 

35 , 63777 

-0.01955 

.01795 


51 

0 

35.65732 


. 02591 



VALUE 


diff 

S I GMA 


H 

M 

SEC 

SEC 

SEC 


21 

34 

5,23340 

,00759 

,00130 


21 

34 

5.22671 

.00090 

, 0 0 0>8 

* 

21 

34 

5,21457 

-n .01124 

. noii>3 

* 

21 

34 

5,22595 

.00014 

, n 0 0 *5 4 

* 

21 

3 4 

5,22582 

•00001 

,00117 


21 

34 

5 , 2334j 

.00759 

,00095 

• 

21 

3 4 

5,22329 

-0 , 00252 

,00043 

* 

2 l 

34 

5 , ?5^5o 

.02375 

.02774 


21 

34 

5,21539 

•0 . 010^^2 

, 0 0 468 

* 

21 

3 4 

5 , 22735 

.00154 

,00090 


21 

3 4 

5 , 21328 

•U . 01253 

. 00232 

♦ 

2 1 

34 

5 , 22o29 

.00047 

,00229 

it 

21 

34 

5,22382 

• 0 . 0 0 1 9 9 

.03178 

* 

21 

34 

5,2287j 

.00289 

,00102 

* 

21 

34 

5,22501 


,00321 



Table 5 (continued) 


RUN 


DATH 

STAR 




M 

D 

Y 



75125- 

10 

4 

14 

72 

59S21344-00 

DEC 

75125- 

71 

5 

29 

72 

59=2134^00 

DEC 

75125- 

16 

6 

27 

72 

59 = 2134^l)0 

DEC 

75125- 

17 

8 

29 

72 

59=2134^n0 

DEC 

75125- 

Iti 

11 

7 

72 

59 = 2l34-*-0 0 

DEC 

75125- 

IV 

2 

4 

73 

59=2134^n0 

DEC 

75125- 

20 

3 

30 

73 

59=2134+00 

DEC 

75126- 

21 

4 

11 

73 

59 =2l34*nU 

DEC 

75126- 

la 

5 

18 

73 

59=2134+00 

DEC 

75126- 

3 

6 

10 

73 

59=2134+00 

DEC 

75126- 

6 

10 

12 

73 

59=2134^00 

DEC 

75126- 

» 

1 

22 

74 

59=2134+00 

DEC 

75130- 

IV 

3 

4 

74 

59=2134+00 

DEC 

75126- 

10 

4 

29 

74 

59=2134+00 

DEC 


>952X34^00 DEC 



VALUE 

DIFF 

SIGMA 

D 

M 

SEC 

SEC 

SEC 

n 

28 

23,98358 

•1,05o83 

,13717 

,04505 

n 

28 

24,82392 

-0.21049 

n 

28 

25,09430 

. 05996 

.13325 

n 

20 

25,14545 

.11104 

,03761 

n 

28 

24,91358 

-n ,12082 

, 0 9o09 

n 

28 

24,83o8i 

-n , 19759 

, 05 736 

n 

28 

25,25139 

.21698 

.04316 

n 

28 

22,32050 

-2.71390 

2.12687 

n 

28 

25.47512 

.44072 

,47832 

n 

28 

25,04716 

.01275 

.09514 

n 

28 

25,64731 

,61290 

, l4iS6 4 

n 

28 

25,20214 

,16774 

.62522 

n 

28 

25,15790 

.12350 

,15431 

n 

28 

25,02373 

-0,01068 

, 06:501 

n 

28 

25,03440 


,23476 


Table 5 (continued) 


148 


RUN 


DATE 




M 

D 

Y 

75125 - 

10 

4 

14 

72 

75125 - 

27 

5 

29 

7 2 

75125 - 

^4 

6 

3 

7 2 

75125 - 

i- 

6 

27 

7 2 

75125 - 

l 7 

tl 

29 

7 ? 

75125 - 


11 

7 

72 

75125 - 

l 9 

2 

4 

7 3 

75125 - 

20 

.5 

3 u 

73 

75126 - 

21 

4 

11 

73 

75126 - 


5 

le 

7 3 

75126 - 

"3 

8 

IQ 

7 5 

75126 - 

6 

in 

12 

7 3 

75126 - 

b 

1 

22 

74 

75130 - 

1 ^ 

3 

4 

74 

75126 - 

1 3 

7 

B 

7 4 

75 l 3 c- 

4 / 

1 

14 

75 


star 


63sVH422201 
o3=VH422201 NA 
63=VK4222q1 KA 
o3=vK4222o1 
oJ = \yR4222Ql RA 
n3=VR'^222ul HA 
ft 3 = V K 4 'i' 2 2 [] 1 HA 
o 3 - ft’ 4 2 2 2 u 1 ft A 
c3=VK4222u1 ha 
f)T = vH4 222ul HA 
^3=VH4222U1 HA 
03=Vft4222ul HA 
ft3=Vft422201 ha 
ft3=;R4222ul HA 
(53 = V ft 42220 1 HA 
t53 = vH42220l HA 

03=VH4222U1 MA 


VALUE 

H M SEC 


DIFP 

SEC 


22 

0 

39,39675 

, 01131 

22 

0 

39.3S659 

, 00146 

22 

0 

39,30823 

. 00230 

22 

0 

39,37uli> 

-0 • 01^27 

22 

o 

39,38715 

f 001^5 

22 

0 

39,38161 

"0 • 0 0^®^ 

2 2 

0 

39.33998 

» 0Q454 

22 

n 

3 9 , 3 7 / 0 u 

-n , C0B44 

22 

0 

39,33960 

•0042c 

22 

0 

39,38595 

. 0 0 052 

22 

0 

39 , 3Bd32 

. 00288 

22 

n 

39,30364 

. 0 0321 

22 

0 

39 , 30393 

-0 . COi^O 

22 

V 

0 

39 , 33i64 

-0 t 00379 

22 

0 

39.380^4 

, 0 033q 

22 

0 

39,38290 

- n * 0 0 2 5 3 

22 

0 

39,38543 



SIGMA 

SEC 

,00170 * 

, ronvo * 

• 0 0 d 78 • 
.nui^2 • 
•000^3 ♦ 

.P01>5 ♦ 

•0011® 

. PUC^9 
. 0 0 U 9 
.00276 
, 00C73 
. PO?6o 
. 0 0 2 4 7 
. n 0 1 2 

f 002^1 

i 0 02<j3 * 
,0 0515 


Table 5 (continued) 


RUN 


Date 

star 




M 

D 

Y 



75125- 

10 

4 

1^ 

72 

63svH422201 

DEC 

75125- 

27 

5 

29 

72 

63=VR422201 

DEC 

75125- 


6 

3 

72 

63=VR422201 

DEC 

75125- 


6 

27 

72 

63 = VR4222(J1 

DEC 

75125- 

l7 

3 

29 

7? 

63=VK4222o1 

DEC 

75125- 

l5 

11 

7 

72 

63=VK4222o1 

DEC 

75125- 

19 

2 

4 

73 

63=vH422201 

DEC 

75125- 

20 

3 

3C 

73 

63=vR422201 

DEC 

7512^- 

21 

4 

11 

73 

63=vR422201 

DEC 

75l2^>- 

1*1 

5 

18 

73 

63 = \/R4222u1 

DEC 

75126- 

3 

a 

10 

73 

63=vR422201 

DEC 

75126- 

6 

13 

12 

73 

63=VK4222q1 

DEC 

75126- 

ti 

1 

22 

7A 

63=VH422201 

DEC 

75i3o- 

l9 

3 

A 

74 

63=vR422201 

DEC 

75126- 

1 3 

7 

Q 

74 

6'5 = VR4222ji 

DEC 

75130- 

4 / 

1 

lA 

75 

63=vK^222qi 

DEC 






63=VR422201 

DEC 

run 


DATE 


STAR 




M 

0 

Y 



75125- 


6 

3 

72 

66sCTAi02 

RA 






665CTA1U2 

RA 

run 


DA 1 S 


star 




M 

u 

Y 



75125- 

1^ 

6 

3 

72 

6CsCTAiU2 

DEC 






66SCTA102 

DEC 



VALUE 


U IFF 

SIGRA 


D 

M 

SEC 

SEC 

SEC 


42 

2 

8.22872 

• 0 , 1 7 0 76 

,01724 

* 

42 

2 

8,39276 

-0,00672 

,00608 


42 

2 

8,30723 

•0 • 0'^225 

f 0 U 7 / 2 


42 

2 

6 , 4 G 4 3 3 

, 0 Q485 

,01533 


42 

2 

8.411)43 

.01100 

, 0 04 7() 


42 

2 

8.3948/ 

-Of 0 0461 

»0lQi^2 


42 

2 

8,43536 

, n3538 

, 0 075 11 


42 

2 

8.42218 

, o22 7q 

1 006l5 


42 

2 

8 ,41593 

» 0 165q 

,02397 


4 2 

2 

6,37709 

-0 ♦ 0*^239 

1 01^18 


42 

2 

8,4187V 

. 01931 

» 0 U 3 2 1 


42 

2 

8,39823 

-0 • 00120 

,01795 


42 

2 

8,46946 

, 09998 

,0293r 


42 

2 

8 , 4] 73i 

» 0l7S3 

♦ 01369 


42 

2 

8 , 4636[j 

,n54i3 

, n 2 1 3 6 


42 

2 

8,42437 

, n2439 

, o29:>5 


42 

2 

6,39948 


,04109 



VALUE 


DIFF 

SIGMA 


H 

R 

SEC 

SEC 

5FC 


22 

30 

7,82732 

,00000 

. 00327 

* 

22 

30 

7,32732 


,00000 



value 


DIFF 

S I GMA 


C 

M 

SEC 

SEC 

SEC 


11 

28 22,49229 

•0.00000 

,28474 


11 

28 22,49229 


,00000 



Table 


5 (continued) 
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RUN DATR star 

M D V 


75125^ 

10 

4 

14 

72 

67s3C454,3 

75125- 

27 

5 

29 

72 

67=3C454 , 3 
67=3C454 , 3 

75125- 

14 

6 

3 

72 

75125- 

10 

6 

27 

72 

j7r3C454 , 3 

75125- 

l7 

8 

29 

72 

67=3L454,3 

67=3C454,3 

75125- 

10 

11 

7 

72 

75125- 

19 

2 

A 

73 

67=30454,3 

75125- 

2U 

3 

3q 

73 

67 = 3-45 4 j 3 

75126- 

21 

4 

11 

73 

67=3 C 454, 3 
67 = 3C454 , 3 

75126- 

10 

5 

16 

73 

75126- 

6 

0 

10 

73 

67 = 3C454 , 3 

75126- 

6 

10 

12 

73 

67=3 C 454, 3 

75126- 

a 

1 

?2 

7 A 

67=3 C 454 , 3 

75126- 

10 

4 

29 

7 A 

67=3 L 454 , 3 
67=3U454,3 


Table 5 


VALUE DIFF 

H M SEC SEC 


RA 

22 

5l 

29,5347V 

f 0 0169 

RA 

22 

51 

29,53364 

* 0 0 053 

HA 

22 

5l 

29,53376 

t 0 0 565 

RA 

22 

51 

29 . 51 v54 

• 0 * 0 1 356 

HA 

22 

5; 

29 , 5346 u 

.00149 

RA 

22 


29,53101 

•0.00210 

HA 

22 

51 

29.53766 

.00476 

RA 

22 

5 1 

29.52997 

-0 , 0 0313 

RA 

22 

51 

29.52929 

•0.00382 

RA 

22 

5l 

29 , 529U6 

-0 . 0 0 4 c3 

RA 

22 

5l 

29,53440 

.00129 

RA 

22 

51 

29.53204 

-0 . 00106 

RA 

22 

51 

29,53210 

•0,00101 

RA 

22 

5l 

29,53421 

.00110 

RA 

22 

5l 

29,53311 



(continued) 


SIGMA 

SEC 

,00122 ♦ 
.C006A * 
,0061^- * 
,00115 # 

, 000^6 ^ 
, 000^1 * 
, C 0 0 ^ * 
, 0 0 3b6 ♦ 
icC?/7 * 
, 0 0 06 A * 
,001^9 cr 
,C03/A * 

• 0 0 0^1 * 

, 0 0 264 
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RUN 


DATE 

star 



VALUE 

diff 

S I GMA 




M 

D 

Y 



D 

M 

SEC 

SEC 

SEC 


75125- 

lU 

4 

14 

72 

67=3C454 , 3 

DEC 

15 

52 

54 , 0689(J 

-0 , l7057 

,02334 

* 

75125- 

27 

5 

29 

72 

67=3C454,3 

DEC 

15 

52 

54,1086/ 

-0.13079 

,01903 

* 

75125- 

!*♦ 

6 

3 

72 

67 = 3C454 ,3 

DEC 

15 

52 

.54 , 23783 

-0,00164 

,46421 

* 

75125- 

16 

6 

27 

72 

67=30454^3 

DEC 

15 

52 

54 ,40508 

. 16561 

, 02024 

* 

75125- 

1^ 

8 

29 

72 

67=30454,3 

DEC 

15 

52 

54,2400/ 

. 00061 

,00545 

* 

75125- 

16 

11 

7 

7^ 

67 = 30454 ^ 3 

DEC 

15 

5? 

^4,l8l2i 

-0.05826 

, 01143 

* 

75125- 

iv ■ 

2 

4 

73 

67=30454,3 

DEC 

15 

52 

54 ,27305 

, 03358 

,01049 

* 

75125- 

2 0 

3 

30 

73 

67 = 3G454 , 3 

DEC 

15 

52 

54 , 25912 

.01966 

,00770 

f! 

75126- 

21 

4 

11 

7? 

67 = 30454 , 3 

DEC 

15 

52 

54,53/11 

.29765 

,17033 

it 

75126- 

16 

5 

18 

73 

67=30454,3 

DEC 

15 

5? 

54 .32150 

.08204 

,07434 

* 

75126- 


8 

10 

73 

67=30454,3 

DEC 

15 

52 

54 , 29v82 

. 00036 

, 03267 


7 5 1 2 6 • 

6 

10 

12 

73 

67=30454,3 

DEC 

l5 

52 

54 , i5645 

■*0*0^3 01 

1 0^320 

* 

75126- 

b 

1 

22 

7 4 

67 = 30454 j 3 

DEC 

15 

52 

54.43924 

.19977 

,13341 

* 

75126- 

lu 

4 

29 

74 

67 = 30454 , 3 

DEC 

l5 

52 

54 , 2*^534 

. 0 0587 

1 02 4 q7 

it 






67=30454 , 3 

DEC 

15 

52 

54,23946 


1 05 367 


Rl.'N 


DATE 


STAR 



VALUE 

D:Fr 

SIGMA 




M 

D 

Y 



H 

M 

SEC 

SEC 

SEC 


75125- 

11 

5 

9 

72 

75=NKA0512 

RA 

16 

36 

48,20011 

,00121 

, U 0 1 3 1 

* 

75126- 

b 

1 

22 

74 

75=NRA0512 

R A 

16 

38 

48,19806 

-0.00084 

, 0 0 503 

* 

75130- 

19 

3 

4 

74 

75=NKAo512 

RA 

16 

38 

48,1944V 

-0 .00441 

,01438 

* 

75126- 

lU 

4 

29 

74 

75=NRA0512 

RA 

16 

38 

48 , 19796 

-0,00094 

,00120 

it 






75=NHA0512 

RA 

16 

38 

46,19890 


,00110 



Table 5 (continued) 
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K UN- 

DATE 
M D Y 

star 


D 

75125- 11 
75126- e 
75130- 19 
75126- 10 

5 9 72 

1 22 74 

3 ^ 74 

4 29 74 

75=NRA0512 
75=NHA0512 
75=NHA0^12 
7 5 = \- K A 0 5 1 2 

DEC 

DEC 

DEC 

DEC 

39 

39 

39 

39 



75=NKA0t?12 

DEC 

39 

RUN 

DATE 
M C Y 

star 


R 

75130- 19 

3 4 74 

80=OU172 

RA 

14 



60=CU172 

RA 

14 

RUN 

DATE 
M p Y 

star 


D 

75130- 19 

3 4 74 

b0=OUl72 

DEC 

in 



eOsOQ172 

DEC 

in 


VALUE 

niFT 

SIGi^A 

M 

SEC 

SEC 

SEC 

52 

3o , 22601 

• n . 00387 

,00782 

,0583a 

52 

30 , 24900 

,01912 

52 

30,11481 

-0.11507 

, 11390 

52 

30,25221 

.02233 

,01830 

52 

30 , 2298ti 


,01205 

VALUE 

DIFF 

S I GM7 

M 

SEC 

SEC 

SEC 

42 

50 .49610 

0 

. 01018 

42 

50 , 49610 


0 

VALUE 


DiFF 

SIGMA 

M 

SEC 

Sec 

SEC 

11 

12,52911 

.00000 

,67097 

11 

12,52911 


, 0 0 0 00 


Table 5 (continued) 


TRPUET 

DATA 

X 

SICMA 

Y 

sigma 

Z 

SIGMA 



CM 

CM 

CM, 

CM 

CM 

CM 

H * 3 • A 

111 

60 ,86 

67,0 

119,33 

CO 

-186.21 

132,7 

TRPlET 

DATA 

X 

SICMA 

Y 

s I gma 

2 

SIGMA 



CM 

CM 

CM 

Cm 

CM 

CM 

H - - s 

222 

24,31 

104,5 

42,55 

107,8 

-217,20 

158.5 

H^\-S 

333 

-16,25 

43,9 

-59.14 

44,3 

-62 ,59 

69,9 

H ^ -’g - 3 

555 

130,07 

47,6 

-98,53 

30,3 

-45 . 09 

102,0 

M-N-S 

666 

56,70 

21.4 

-18 , 92 

22,5 

-119.20 

49,4 

H-N-S 

777 

57,50 

20,6 

-82,17 

28,0 

130,93 

65,8 

TRPlET 

DATA 

X 

SICMA 

Y 

sigma 

2 

SIGMA 



CM. 

CM 

CM 

CM 

CM 

CM 

j 

22? 

-25,68 

33,8 

71.63 

100,0 

51,62 

65,6 

H-n-6 

333 

-30,23 

120.9 

255,95 

117,9 

174,01 

151,9 

H - N - G 

555 

10,35 

48,1 

80,03 

54,3 

32,27 

103,5 

H-'J-G 

666 

-27,03 

25,7 

134.18 

39,6 

65.53 

52,6 

I 

1 

1 

<7> 

777 

-27,95 

38.2 

64,40 

51,8 

16.81 

65,2 

h^n-G 

688 

3,16 

^,8 

69,43 

19,7 

20.24 

21,0 

TRPLET 

DATA 

X 

SICMA 

Y 

SIGMA 

Z 

SjGMA 



CM 

CM 

CM 

CM 

CM 

CM 

H ^ 3 “ G 

222 

-3,52 

43,8 

71,62 

76,3 

151.79 

93.2 

H-3-G 

333 

21,31 

247,5 

448,64 

174,3 

351,93 

254,2 

H-S-G 

44A 

11,64 

72,0 

-81,82 

78,7 

-1323,44 

1445 , 8 

H-S-G 

555 

-151,06 

86,5 

311,09 

125,2 

104,84 

197,8 

H-3-G 

^66 

326,23 

73,4 

113.10 

73,7 

3o4,42 

176.9 

H-3-G 

m 

-66,81 

44,8 

131,31 

56,5 

30.79 

175,1 

TRPlET 

DATA 

X 

sigma 

Y 

SIGMA 

1 

SIGMA 



CM 

Cm 

CM 

CM 

CM 

CM 


Table 6 


Baseline Closure Errors 
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TRPlET 

data 

X 

sigma 



. CM 

CM 

G • N “ S 

222 

46 , 47 

110.6 

G • - S 

33.3 

35,20 

24^,7 

G-N-S 

555 

-31,34 

96,5 

G-N-S 

666 

411,95 

74 , 6 

G • iW • S 

ni 

16,64 

45.4 


Y 

SIGMA 

Z 

SIGMA 

CM 

CM 

CM 

CM 

42.53 

150,6 

-117,03 

180,9 

133,55 

176.3 

115,33 

256,3 

132,53 

125,. 3 

27.48 

206,8 

-40,00 

76,9 

179.69 

180 , 0 

-14,77 

59,3 

142,90 

177,9 


Table 6 (continued) 


TRPLET 

data 

X 


sigma 

Y 

sigma 

2 

sigma 



cv. 

CM 

CM. 

CM. 

CM 

CM 

h-n-s 

223 

-28 

.ne 

lOii.4 

-31,55 

104,2 

-187.28 


H*N-S 

224 

-UO 


105^,5 

"117,30 

104,3 

-318.51 

177.8 

M'N-S 

225 

70 

.36 

lOii.4 

-126,81 

103,9 

-182.47 

101,2 


226 

-5 

. 79 

1 0 2 , Z 

-87,00 

104,2 . 

-231.95 

354,5 

P-N-S 

227 

•21 

1 7 V 

101.7 

-102, 95 

104,1 

-150,06 

154,7 

- • N - S 

232 

5U 

.95 


-14,73 

51,0 

-106.27 

75.6 

K - N - S 

233 

‘ 1 

.42 

42,8 

-88,83 

42.9 

-76.35 

68,5 

T 

• 

1 

to 

234 

-83 

. 92 

43,^ 

"3 74.58 

43,0 

-207.58 

110,6 

h - r j - S 

235 

9/ 

.5C 

42,9 

-184,09 

42,1 

-71,54 

81,2 

P-M-S 

236 

20 

.85 

41,6 

-144,20 

42,8 

-121.02 

66,8 

P-N-S 

237 

4 

.91 

41,3 

"160,23 

42,5 

-39.13 

67,4 

P-N-S 

252 

40 

.85 

49,3 

-41,89 

40,3 

-140 , 92 

88 , 9 

P - \ - s 

253 

-11 

• 5 £ 

44,6 

-11'’. 98 

29,4 

-111.00 

82,9 

P - N - S 

254 

-94 

.02 

44,9 

-?ni,73 

29,6 

-242.23 

120.0 

H - N; • 5 

255 

87 

.4C 

44,7 

-211 .25 

23,1 

-106.19 

93.7 

P-\-S 

256 

lU 

.75 

43,7 

-371.44 

29,2 

-155,67 

81,5 

P-N-S 

257 

-5 

.20 

43,2 

"187,39 

28,8 

-73.78 

81 , 9 

P**N-S 

262 

63 

.26 

31.1 

53,73 

35,6 

-127,48 

60,1 

P N • S 

2t3 

15 

,86 

23, C 

-20,36 

22,6 

-97.56 

50,9 

P-N-S 

204 

• 6 6 

. oc 

2>3.7 

- 3. n 6 , 1 1 

22,8 

•228.79 

100,6 

P-M-S 

265 

1! 4 

,3C 

23.1 

-115.63 

20.9 

-92.75 

67,0 

P»*N-S 

266 

38 

.U 

21.1 

-75,81 

22,3 

-142,23 

48 , 6 

P * N • S 

207 

22 

.23 

20.1 

-91.77 

21,8 

-60.34 

49 , 3 

P-N-S 

272 

6U 

.77 

31,2 

30,86 

39,6 

48,37 

71,7 

H-N-S 

273 

8 

,40 

2-M 

-43,23 

28,4 

78.29 

64,1 

H-N-S 

274 

-74 

.10 

23.6 

-128,98 

28,5 

-52,94 

107,9 

P - N' - S 

275 

10 7 

.32 

23.3 

-138,50 

27,0 

83,10 

77,5 

P-N-S 

2 76 

30 

.67 

21,3 

-98,69 

28,2 

33.62 

62,3 

P-N-S 

277 

14 

, 73 

20.3 

-114,64 

27.7 

118.51 

62,9 

P - N - S 

322 

9 

,48 

104,9 

72. ?3 

108,4 

-203,44 

159,1 

P-N-S 

323 

-42 

, 89 

10 2.8 

-1 ,86 

104,8 

-173.52 

155,8 

P-N-S 

324 

-125 

. 3 9 

103.0 

-87,61 

104.9 

-304.75 

178,4 

P-N-S 

325 

56 

,03 

10 2,6 

-97,13 

104,5 

-168,71 

161.8 


Table 7. Baseline Closure Errors 
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TRPLET 

DATA 

X 

SIliMA 



CH 

CM 

H • \ • G 

223 

148.15 

91.4 

H-N-G 

224 

-115,15 

45.8 

H-N-G 

225 

121. T7 

42 , « 

H-N-G 

226 

4 , 4 C 

35,3 

u-M-G 

22 7 

-3C ,55 

4 1,7 

h‘ • N - G 

22H 

14,73 

31,4 

h-N-G 

232 

-169.23 

85,4 

H-N-G 

233 

-15, 4C 

120,5 

H - .N - G 

234 

• 2 7 1 1 7 C 

9 0.9 

H-N-G 

235 

--^2 , 

P9.2 

h-N-G 

2 3 6 

-159,15 

86,1 

H - N - G 

237 

-194 , l3 

86 . V 

H-N-G 

238 

- 1 4 6 ,62 

84,6 

M - N - G 

262 

-179 , ll 

36 , 9 

h-N-G 

26 3 

-5,24 

9^.6 

H-N-G 

254 

-266,54 

46,2 

H - N - G 

255 

-32,32 

45.8 

H-N-G 

256 

-14b , 99 

36,4 

H-N-G 

257 

-163,96 

44.3 

H-N-G 

256 

-136 , te 

34,8 

H-N-G 

262 

-V5,t5 

2«J,3 

H-N-G 

263 

9 6 , 1 a 

86 , 1 

H-N-G 

264 

-165,12 

36 , 8 

H-N-G 

265 

71,10 

34,7 

H-N-G 

266 

-45,17 

25,5 

H-N-G 

267 

-60,55 

33.6 

H-N-G 

268 

-35 ,24 

19 , 8 

h-N-G 

272 

-65,62 

29.1 

H-N-G 

273 

1^6.01 

89 , 6 

H-N-G 

274 

-155,29 

42,5 

H-N-G 

275 

6 U , 5 3 

36,6 

H-N-G 

276 

-35,74 

. 9 

H-N-G 

277 

- 7 U , V 3 

.3 6 , .1 


Y 

SIGMA 

1 

SIGMA 


CM 

CM 

CM 

30 t ,37 

125,8 

163.52 

124,0 

104.35 

97 , 6 

239.83 

184,9 

-21,55 

101,2 

63.09 

86,2 

116.14 

95,7 . 

204,35 

69,6 

202.12 

99,8 

206,52 

78,8 

-75,84 

94,0 

21.37 

61,7 

-5.4b 

89.1 

48.35 

108,7 

226.26 

117,4 

160 ,2i' 

151,3 

27.24 

86 , 5 

236,56 

204,2 

-98.66 

90,5 

59,82 

122,2 

39.03 

84 , 4 

201 . C8 

111 , 1 

125,01 

a'^.o 

203,25 

117.1 

-152 , 95 

82,4 

18,10 

106,4 

60,50 

50,7 

-40.30 

77,3 

292,24 

91 . 7 

71.60 

130,6 

93.22 

45,8 

147.91 

189,4 

-32.68 

53,1 

-28.83 

95.4 

1 r. 3 . n 1 

41 , 8 

112.43 

00.7 

19C .99 

50.4 

114,60 

88.0 

-86.97 

37,5 

-70.55 

74 , 0 

32.78 

40,8 

-110,23 

46,3 

264 ,G2 

90,7 

1.67 

115,0 

65.50 

43,7 

77.90 

179,0 

-60.40 

51,3 

-98.76 

72,6 

77,29 

. 39,5 

42.50 

51,8 

163,27 


44,67 

63,7 

-114,69 

34,9 

-140.48 

40,6 

-98.55 

52 , U 

-151.51 

44.3 

133.19 

92,5 

-39.61 

114,2 

-65,03 

47,3 

36. 70 

178.5 

-1.91,73 

54,4 

-140,04 

71 . 3 

- 5 4 ,04 

43,4 

1.22 

50 . 0 

31.94 

51.7 

3.39 

62,2 
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TPPUET 

lata 

X 

S I G A 

Y 

sigma 

Z 

sigma 



cx 

CM 

CM 

CM 

CM 

CM 

M*»S-G 

2c'3 

17U.31 

95.6 

303,36 

100.0 

263,69 

140,6 

H-S-G 

224 

-P2 . 78 

53.6 

104,34 

73,1 

340.00 

196 , 4 

H-S-G 

225 

1^3 .23 

5J,7 

• 21.56 

77,9 

163.26 

108,7 

H-S-G 

226 

26,36 

45.0 

■16,13 

70.7 

304,52 

96 , 1 

H-S-G 

227 

-8,42 

50.2 

202.11 

76,1 

306,69 

102,9 

H-S-G 

22«^ 

36,09 

42,0 

-75, C5 

68,2 

121 .54 

90,5 

H - S - G 

232 

- 2 n 4 , d 9 

233,4 

142.80 

159.0 

269.95 

233 , 6 

H-S'G 

2 33 

-3i , n6 

24d,4 

374 . S4 

176.5 

381,85 

256,2 

H-S-G 

2 54 

-2V4, 36 

235,5 

175.52 

157,6 

458.16 

290 , 6 

H-S-G 

235 

-bd, I4 

234,0 

49.62 

159.8 

281 ,42 

240,2 

H-S-G 

236 

-174,01 

233,7 

187,31 

156,4 

422.68 

234,0 

H-S-G 

237 

-207,30 

234,7 

273.29 

159,0 

424,05 

237, 7 

h-S-G 

238 

-164,40 

233,1 

-4,67 

155 . 4 

239.70 

232,5 

H-S-G 

?42 

-33,77 

66,1 

-274.30 

>36 , 1 

-1612,96 

1433,1 

H-S-G 

243 

1 4 J , 0 6 

106,9 

-42,64 

115,1 

-1501 , 06 

1437,0 

H-S-G 

244 

-123 , 23 

74.0 

-241,66 

8 3,3 

-1424.75 

1443,5 

H-S-G 

245 

112,98 

72,0 

-367 ,h6 

8 7,6 

-1601,49 

1434,2 

H-S-G 

246 

-3,69 

60 , 9 

-229,87 

01,2 

-1460 , 23 

1433,3 

H-S-G 

247 

-3d ,67 

72,4 

-143.89 

85,9 

- 1 4 •> 8 . 0 6 

1433,8 

H-S-G 

248 

6,64 

6 6,9 

-421,85 

79,1 

-1643,21 

1432,9 

H-S-G 

252 

-251,27 

85,2 

: 3 4 . 9 1 

127.5 

128 ,10 

187 , 4 

H-S-G 

233 

-77,44 

120,3 

466.65 

148.7 

24 0 . UO 

214,9 

H-S-G 

234 

-340,73 

90,6 

2 6 7 ,6,1 

125,7 

316,31 

254,9 

H-S-G 

2*55 

-104,32 

80,9 

141.73 

128,5 

139.57 

195.5 

H-S-G 

236 

-221.19 

85,0 

279,42 

124,3 

200.83 

150,8 

H-S-G 

2'>7 

•256,17 

00 , 6 

365.40 

127,4 

253 . 00 

192.4 

H-S-G 

2od 

•21J,fc6 

64,3 

87, 44 

122,9 

97.85 

186.0 

H-S-G 

262 

26d.f;3 

7c , 1 

-60.96 

d3,e 

196,94 

178,9 

H-S-G 

23 3 

4«i ,ee 

114,1 

170,70 

113,4 

3 0 8 , 8 4 

207,6 

H-S-G 

264 

170,38 

82,2 

-28,24 

d0,9 

385.15 

240,0 

H-S-G 

2^5 

414,00 

80.3 

-154,14 

85.3 

208.41 

107,4 

H-S-G 

2 66 

293,15 

76,8 

- 1 , 4 3 

78,7 

349.67 

180 , 4 

H-S-G 

2 a 7 

263,14 

80.0 

69,53 

83,6 

351,84 

184.2 
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TRPLET 

data 

X 

SIGMA 

Y 

sigma 

Z 

SIGMA 



CM 

CM 

CM 

CM 

CM 

CM 

G-K'-S 

2?.:i 

-154,91 

254,6 

113.72 

205,3 

1 , 13 

280 , 4 

G - M - S 

22^ 

16.22 

122.3 

-303,47 

155,8 

-1881 . 70 

1441 . 5 

G-N-S 

225 

-2ni . 28 

132,6 

205,83 

182,0 

-140,72 

243.2 

G - M - S 

226 

3 1 b , n 3 

126,9 

-90,05 

154,5 = 

-71 . 88 

236 , 8 

G-N-S 

22 7 

-57,97 

110,9 

-173,27 

144,4 

-325,80 

232,3 

G-N-S 

232 

73,11 

55,9 

-14,75 

116,8 

•6,10 

115,5 

G “ N - S 

233 

-128,27 

237,0 

56,44 

162,0 

112,06 

243 , 4 

G " N - S 

23 

42,86 

7V.4 

-360 , 75 

3 23.5 

-1770,85 

1434 , 7 

G"f.'-S 

2 35 

-174,64 

94 , 5 

148,55 

155,2 

-29 , 79 

199 , 4 

G “ N - S 

2 3 6 

344,67 

86 , 4 

-147,33 

121,9 

39 , 05 

191 . 5 

G-N-S 

237 

-31,33 

60,4 

-230.55 

108,8 

-214,87 

186 , 0 

G-N-S 

25 2 

63,02 

61,2 

-41,91 

112,6 

-40.75 

124 , 6 

G - N - S 

253 

-136 , 36 

237,3 

2'‘^ , 2. 

179.3 

77.41 

247,8 

G-N-S 

25 4 

32.77 

80 . 4 

-387,91 

119.4 

-1305,50 

1435,5 

G - \ - S 

255 

•1«4,73 

95.3 

12J .39 

152,0 

-64,44 

204,8 

G-N-S 

25h 

334,58 

87,3 

-174,48 

117,6 

4,40 

197,1 

G - N - S 

25 7 

-41 . 42 

61,7 

-257.70 

104,2 

-249.52 

191 , 7 

n - N - s 

262 

y 0 . 4 2 

47,8 

53,72 

111 , 0 

-27,31 

106,0 

G - \ - S 
G-N-S 

26 3 
26-4 

-lie . 95 
6C .17 

234.2 

70,8 

124,91 

-292,28 

178,3 

117,9 

90,85 

-1792,06 

239,0 
1434 , 0 

G-N-S 

G-N-S 

2 6 5' 
266 

-157,33 

361,96 

87,3 

76,5 

217,02 

-73.86 

150,9 

116.3 

-51,00 

17.84 

194 . 0 
165 , 9 

-i - N - S 

267 

-14,02 

48.5 

-15:?, oa 

:02,5 

-236 . 08 

180 , 2 

G - -N - S 

27 2 

82,94 

47.9 

30,64 

312,3 

148,54 

112 , 9 

G-N-S 

2 73 

-118,44 

234,2 

1!)2,03 

179,1 

266 . 7o 

242 , 1 

G-N-S 

274 

52,65 

70 , 8 

-3lf .16 

119,2 

-1616.21 

1434,5 

G-N-S 

21b 

-164 , ai 

87,4 

194,14 

351,8 

124 . 85 

197 , 9 

G-N-S 

276 

354,50 

78,6 

-11)1.73 

317.5 

193 . 69 

189 , 9 

G-N-S 

211 

-21 , 5P 

48.5 

-184,95 

10 3,9 

-60,23 

184 , 4 

G-^. -S 

322 

2 1 U , 0 2 

135.6 

119,64 

3 43,6 

•113 . 76 

200,6 

G - f . “ S 

32 3 

8,o4 

266 , 4 

190.83 

200.2 

4,40 

293 , 5 

G-N-S 

324 

179,77 

143.3 

-22* ,36 

3 49.0 

-1878,51 

1444 , 1 

G-N-S 

325 

-37.73 

154,1 

231'. 94 

176,3 

-137.45 

258,2 


able 7 (continued) 


station 

NRAO 


Onsala 


* 

The measured 

** 


TA BLE 8 

Antenna Axis Offset Results * 

Data** Axis Offset 

Measured- VLB I Estimate 


1 

-0.55 

± 

0.46 m 

2 

-0.03 

± 

0.29 

3 

0.34 

± 

0.67 

5 

0.18 

± 

0.55 

6 

-0.11 

± 

0.29 

7 

-0.39 

± 

0.34 

2 

-2. 

± 

3. 

3 

5. 

± 

5. 

4 

11 

± 

.6. 

5 

11. 

± 

7. 

6 

-10. 

± 

3. 

7 

6. 

±10. 


axis offset values are given in Section 4.5.5 


The "data" codes are given in Section 4.5.4 



TABLE 9 



( 


Estimates of Corrections to BIH 
Values for Polar Motion and Al-UTl 


Date 



Polar Motion 





AX Comp . * 

AY Comp . * 

A (Al-UTl) 



(0?01) 

(OVOl) 


msec 

Phase I 






Apr 14-15 

1972 

- 7.0 ±0.7 

— 

- 

1.1±0.7 

May 9-10 

1972 

3. 8 ±1.6 

- 15.5±10. 


1.1±1.5 

May 29-30 

1972 

2.0 ±0.9 

- 12. 7± 9. 


2. 9 ±1.2 

June 6-7 

1972 

5.0 ±1.8 

— 


1.2±1.4 

Aug 29-30 

1972 

— 

— 


— 

Nov 7 

1972 

- 1.1±0.8 

— 


1.2 ±0.8 

Feb 4-5 

1973 

- 3.2±0.9 

— 

- 

1.2±0.8 

Mar 30-31 

1973 

4.4±0.7 

— 


5.7±0.6 

Jan 16-17 

1975 

0.6±2.2 

— 


5.4±1.3 

Weighted mean and 
rms scatter about 
weighted mean 

- 0.6±4.3 

- 14. 0± 1.4 


2.0 ±2. 8 

Phase II 






Aug 10-14 

1973 

3.3±9.6 

4.1±16.2 

- 

5.1±6.4 

Oct 12-16 

1973 

- 8.2±6.0 

5.6± 8.1 

- 

2.5±1.4 

Jan 22-23 

1974 

— 

— 


— 

Mar 4-*8 

1974 

- 7.7±3.1 

- 3.8± 4.1 

- 

4.2±0.8 

Apr 29- 
May 3 

1974 

- 2.9±4.4 

1.1± 6.4 

— 

2.6±1.7 

July 8-12 

1974 

1.2±3.1 

6.1± 5.5 

- 

2.2±1.1 

Wei ^ ted mean and 
rms scatter about 
weighted mean 

- 3. 4±4.2 

0.8± 4.4 

— 

3.3±0.9 


* 

To convert these angular measurements into meters at the 
surface of the Earth, note that OVOl-0.3 m, and 1 msec ~0.5 m. 
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APPENDIX A 


USE OF DIFFERENCED OBSERVABLES 

In Chapter 2 the use of differenced observations to elim- 
inate systematic clock errors is discussed. Here we wish to 
ar.alyze in detail the problem of the determination of baselines 
and source coordinates from differenced observations. 

Let us suppose that our receiver stations are capable of 
making simultaneous observations on two separate sources and 
that the observations we wish to deal with will be the difference 
between the delay for one source and the delay for the other. 

(We will discuss the use of differences of non-simultaneous ob- 
servations later.) Following the discussion in Chapter 2, we 
notice that the delays from observations of a given source follow 
a diurnal sinusoid, therefore the differenced observations also 
follow a diurnal sinusoid, and can be characterized by three 
independent parameters, e.g. phase, amplitude and offset from 
zero. (The frequency is assumed to be known.) Each source pair 
therefore contributes three independent observations. There are 
six parameters that we wish to determine from these data: 
three components of the baseline, and two coordinates for each 
source, less the right ascension of one source. If we try to 
improve the solution by adding a new pair of sources, we add 
four unknown parameters and only three independent observations 
for each new pair. If, instead, we keep one source the same 



I 


and only add one new source, we add only two unknown parameters, 
but still we get three new independent observations. It is 
easy to show that with a single reference source and four other 
sources paired with it the number of unxnowns is equal to th® 
number of independent observations. 

For comparison purposes, we should note that with undif- 
ferenced observations, three sources were sufficient to deter- 
mine the baseline and source coordinates. The fact that five are 
needed now is a measure of the real information that is lost 
along with the clock effects in the process of forming differ- 
enced observations. 

For real observations it is possible, of course, to observe 
a reference source continuously only if the stations are all 
located on the same side of the equator and the magnitude of the 
declination of the reference source is greater than the co- 
latitude of the station furthest from the pole. Since it is not 
always possible to meet such conditions, it is necessary to in- 
vestigate the conditions under which more than one reference 
source can be employed. Each new reference star adds two un- 
knowns. Each new non-reference star adds two unknowns. The 
baseline has three unknowns. The number of independent ob- 
servations is three times the number of non-reference stars. 

The necessary condition for a unique soluti- i is that the num- 
ber of independent observations be greater than or equal to the 
number of unknowns; in other words: 
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3Ng 2 3 + (2Nj^ - 1) + 2Ng (A-1) 

or 

Ng 2 2 H 2Nj^ (A-2) 

wherv^ 

= number of reference stars 
Ng = number of non-reference stars 

Thus if three reference sources are required, at least eleven 
sources must be observed in all. Eleven sources is only the 
minimum number required to determine the baseline and source co- 
ordinates. Again, in actual practice, it is highly desirable to 
far exceed the minimum number of sources required in order to 
provide redundancy and internal consistency checks. None of the 
observation schedules that we have employed have called for ob- 
servation on a number of sources far in excess of eleven. 

Up to this point we have considered only the case of pairs 
of simultaneous observations. If the observations being dif- 
ferenced are not exactly simultaneous, then we must consider the 
problem of clock errors again. If we write the delay a- the sum 
of the geometric delay plus a polynomial for the difference in 
clock errors for the two sites, 

T . = Tp. + a + 0 t . + ytf + . . . (A-3) 

i. uX X X 

then the differenced observable can be written as 

^2 ‘ *^1 " ^G2 “ *^G1 ^ “ ^1^ + Y(t2 - t^) + ... (A-4 
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If the magnitudes of the high order derivatives of che clock 
error are not negligible, it is a simple matter to add the ap- 
propriate polynomial coefficients to the set of parameters to 
be estimated. In order to be able to estimate these additional 
parameters, it is necessary to add to the above minimu^Ti criteria 
a single observation for each clock term required in the solution 
Again, as always, a large surplus of observations beyond the 
theoretical minimum is desirable. 
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APPENDIX B 

EXPRESS 10^^ OF 7aOMIC TIME ^3 FUNCTION OF COORDINATE TIME* 


B.l INTRODUCTION 

In this Appendij. w' derive a simx"d.e expression for 
the readings of atomic clocks in teruis of our choice of 
coordinate timt; (see Section 1.2). As stated, our context is 
th 2 solar system \»hose dynaiiu.cs cim bs described utilizing the 

weak-field solutiorj to Einstein's field equat.ions. Tne 
equations of motion so derived can be 

cribed in tei-ms of harrr.onic coordinates (sac- Reference ?o i 
p. 162) , with origin at the solar-syste.a bar /center. As 
with any choice of coordinates . these are not "obser- 
vables" (i.e., they are nvot directly measurable in actual 
experiiponi.s) , Since all of our calculations are carried 
out with coordinates, we must derive convenient ex- 
pressions for tli 2 observables in terms of the coordinates. 

VvV. restrict ourselves here to the coordinate expression 
for the readings of Earth-based atomiic clocks such 
as those utilized in our very-long-Laseline -inrer- 
I'erometry (VLBI) experiments. 

♦This work was carried out: in collaboration with R. D. 
Reasenberg . 
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As it will be necessary in the following analysis 
to distinguish among many apparently similar quantities, 
the following conventions will be employed 

t = coordinate time 

AT = time kept by an atomic clock on the Earth 
= displacement of the center of mass of a 

gravitating body, i, from the center of mass 
of the Earth 

-► _ 

X = three-vector displacements not covered by other 
conventions; defined in diagrams as used 

X = derivative of x with respect to t 

— ► ^ 

v = three-velocity of a body with respect to the 
solar-system barycenter 


B.2 ANALYSIS 

B . 2 . 1 General Formulation 

In the post-Newtonian approximation to 
the solution of the Einstein field equations, the 
assumption is made that one is dealing with a system 
of slow-moving particles (v << c) bound together by 
their mutual weak gravitational interactions, a system 
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like our solar system (see Reference 20, p. 212). 

In this approximation, the relation betv/een time kept 
by an arbitrarily moving atomic clock and coordinate 
time, t, is given by 

2 

AT - t ^ fi<t> - ^)dt (B.l) 

where terms of order v^ and higher have been neglected 
[see Reference 20, p. 221, Equaticns (9.2.2 and 9.2.'^)^* 

Note that in the system of units used here, c = 1. In 
other \^rds, the unit of length is the light-second, as in Section 2.3.1. 
We can write 

AT - t = at + e + k (B.2) 

per 

where represents the periodic terms in the solution. 

In this Appendix we will be concerned solely with the 

e portion of this solution. Constant terms in the 
per 

integrand and in the integral will be neglected as 

they arise. Treating the rate in this manner neglects 

2 

a constant factor of order (1 + v ) in the amplitude 
of the pel Iodic terms. Itie deviation from unity of this factor is 
the same order as other terms which will be neglected. 

B 2.2 Separation into Diurnal and Non-diurnal Terms 
In order to integrate Equation (B.l) for 
clocks associated with earth observatories, it has 
been found useful to separate the integrand into two 


mm 
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types of time-varying terms. The first type depends 
on the displacement of the atomic clock from the center 
of the Earth. These terms tend to have diurnal period- 
icity, and will be referred to as diurnal terms. The 
second type depends on the displacement of the center 
of the Earth from the solar system barycenter and 
will be referred to as non-diuinal terms. 

In order to separate the velocity terms, we write 


(Figure B.l) 


v = V + X 
ec r 


(B.3) 


viiere v^^ is the velocity of the center of mass of the Earth, and 
vdiere, in this weak-field aj^oidnation, we neke the further approxi- 
mation that these vectors can be added in a Euclidean manner. Then, 


o .2 

2 V X 

v ec , r , 

2 2 2 ec 


(B.4) 


Of these three terms, the first is a non-diurnal 
term, the second is, for all practical purposes, a 
constw.n\:, and the third is a diurnal term. 

In order to separate the potential terms we make 
a Taylor expansion about the center of mass of the 
Ihrth (Figure B.l) 


'•’i i ^ 2 H. 1/2 

(r^ + x^ - 2r^ • x^) 


where is the contribution of body i to the potential, is the 
ness of body i, and G is tfie gravitational constant. 


1 


r 




Only the second term has a diurnal signature. 







If the la MS of motion are taken to be ap^aroximately Newtonian 



(B.8] 


which is sufficiently accurate for the correction terms, 
then the integral takes the simple form* 


(AT - t) 


diurnal 



(B.9] 


Before analyzing the approximations employed in 

deriving Equation (B.9), we consider the approximation 

involved in Equation (B.l). As can be seen from 

Reference 20, p. 221, Equation (9.2.3), the largest 

terms omitted from Equati<^n (B.l) are smaller than the 
2 “8 

terms kept by v - 10 . Since the expression in 

Equation (B.9) represents a near-sinusoid wit,h ampli- 
tude of 1.5 visec, the dxurnal terms neglected in 
Equation (B.l) should have an effect on the order of 
1.5 X 10 sec. 

Turning to the approximations made in Equation (B.7) 
we see that the largest short-period (semi-diurnal) 
term neglected is contained in 



* 

This result was observed independently by Thcmas (Reference 19 ) . 



If we expand, integrate and evaluate using 


GM - 5 X 10 ^sec, I r . I ~ 5 x 10^ sec, lx I ^ 2 x 10 ^sec, 

0 ' X ' ' r ' 


the eunplitude of the semi-diurnal term is seen to be 
-14 


7 X 10 sec, so Equation (B.9) should be accurate 
to a fraction of a picosecond. 

There is another assumption hidden in Equation (B.9). 
We have considered the gravity field of all the bodies 
in the solar system external to the Earth. To be 
more accurate, the gravitational potential of 

Equation (B.7) should have contained a term of the 
form 

GM^ 





(B.ll) 


For a rigid, spherical Earth, this term can be 
neglected as it is constant. However, the Earth tides 
cause an effect of the order of 


A I y I -- ^ rj ^ Q 


-9 


(H.12) 


and thus a vn;. potential at 


A(|) = 


"'''e 


-j A I j ' 4 X 1 a ^ 


(P.-.I35 


which is of the same order as the error represented 
by Equation (B.IO). 


Another question which needs to be considered is 
the effect of the non-sphericity of the Earth. This 
introduces a position-dependent variation in the Earth’s 
potential, which, because it is constant on a short time 
scale (<< 10^ years) will not be treated further here. 


B.2.4 Integration of the Non-diurnal Terms 

The terms of Equation (B.l) which are neither 
diurnal# semidiurnal, nor constant are 


(AT - t) = /[- 

non-diurnal 


GM. 

- I ]dt 

^ i l?il 


(B.14) 


There is, at present, no known simple expression for 
this integral for the general case of N gravitating 
bodies. There are, however, some interesting and use- 
ful special cases which can be handled. For the two- 
body, Earth-Sun case (Figure B.2) 

AT - t = /[ ^ ]dt (B.15) 

With the diurnal result in mind, we seek a 
solution of the form: 


rr A 


t P,v 

1 ec 


ec 


(D.16) 


where is a constant to be determined. If we 
differentiate Equation (B.IG) 
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(B.17) 


^(AT - t) = P,v^ + P,x • V 
dt 1 ec 1 ec ec 


and again take the Newtonian laws of notion [Equation (B.8)] 


^(AT - t) = P-v^ + P ,x 
dt 1 ec 1 ec 


GM 


Iv 


and since 


r 

M © 

© 


(B.18) 


ec 

we f ind 


r = -Kr 


(B.19) 


^ , GM 

^(AT - t) = P.v^^ - P,K 
dt 1 ec 1 1 ^ 

0 


(B.20) 


Now for the case of a two-body orbit, Newtonian conser- 
vation Df energy requires 


V® 


.2 

r 

-Q. 

2 


G(M + Mj 
0 ^ 


= E. 


(B.21) 


or 


2K 


ec 

T 


G («0 + 


= e; 


(B.22) 


Where and are constant. Combining 
Equations (B.20) and (B.22) 
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Now by equating the coefficients of the and 


(B.23) 


the 


0 


terms with the respective coefficients of 


the integrand of Equation (B.15), we are able to 
determine P^^. 


^ = - 


2Me + 
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(B.24) 


and thus 


(AT - t) 


2«e + «( 


two body 
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0 
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ec ec 


- fP2^t (B.25) 

and since we have defined t in such a way as to eliminate 
linear terms like the second term on the right-hand 
side 
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with 


K 


2Mg + 


e © 


M 


© 


= (1 t ±) 


(B.27) 


Another interesting case which we can try to 
handle is the Earth-Moon-Sun system. The geometry 
of the situation suggests that we resolve the problem 
into two two-body problems and perhaps a cross-term. 
We write Equation (B.l) as (see Figure 3.3) 
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If we break up the first potential term as before, 

and expand the velocity term, we arrive at 
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The third term has exactly the form of the diurnal 

problem [Equation (B.8)], and has a solution 

-X , • X* . The first and second terms have the form 

emb ec 

of the problem in Equation (B.15). The solution to 
Equation (B.15), however, cannot be directly applied 
in this case because the conservation of energy 
equation [Equation (B.21)] doesn't apply to a three-body 
situation. That is to say, the energy interactions 
are too complicated to be expressed in simple closed 
form as in Equation (B.21). However, we can ignore 
this problem and treat the Moon's orbit and the Earth- 
Moon barycenter orbit as though they were separate 
two-body orbits, and then estimate the error introduced 
by neglecting the energy interaction. We then arrive 
at 


(AT - t) 


e-m-s 


= -X 


emb 


"emb 
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X* ‘X* • K__ 

ec ec 0© 

X , • x' 

emb ec 


(B.30) 


We will consider the errors in this expression 

term by term. The error in the third term [See 

-12 

Equation (B.9)) has an amplitude of 2.7 x 10 sec. 
The error in the second term is associated with the 
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assumption constant energy for the Earth-Moon orbiv:. 
However, the amplitude of this term is only about 
3.6 X 10 seconds. An error in this term of less 
than ten percent will be less than the error in the 
third term, and thus we can neglect errors in term 2. 

In order to estimate the errors in the first term, we 
need to estimate the energy interaction between the 
Earth/Moon system and the Earth-Moon-barycenter/Sun system. 
If we caDculate the magnitude of the quadrupole com- 
ponent of the potential of the Earth-Moon system at 
the position of the Sun, and compare it with the mono- 
pole component of the potential, we find that the qua- 

7 

drupole term is smaller by a factor of 10 . This 
would indicate that errors in the first term of Equa- 
tion (B.25) caused by the quadrupole interaction of 
the Earth-Moon system are of the order of 3 x 10 sec. 
However, in order to use Equation (B.25) for a clock 
on the real Earth, some attention must be paid to the 
perturbations of the Earth's orbit by other planets. 
According to I. I. Shapiro (private communication, 1973), 
studies by numerical integration indicate that the short 
term variations in the energy of the Earth's orbit are 

5 

no more than a part in 10 . This would indicate that the 

first term in Equation (B.25) is accurate to only 
-9 

3 X 10 sec. 
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B,3 CONCLUSION 

If v;e collect the terms of Equations (B.9) and 
(B.29) r we arrive at 
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The first term, an 

-9 

annual sinusoid, should be accurate to 3 x 10 sec. 

The second and third terms, monthly sinusoids, should be 

-12 

accurate to 3 x 10 sec. The fourth term, a diurnal 

-13 

sinusoid, should be accurate to 2 x 10 sec. 

It may well be asked why there is any use in having 
error limits on the short-period terms which are orders 
of magnitude smaller than the errors in the longer 
period terms. The answer is that for some types of 
observations, smaller error limits are needed only for 
the shorter period terms. In VLBI work, for example, 
everything but the longitude -dependent diurnal term can- 
cels out of the observations to very high accuracy. 

It is hoped that Equation (B.30) will be of some 
use in processing various types of astronomical observa- 
tions. 
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APPENDIX C 

ORGANIZATION AND PRESENTATION OF RESULTS 


The results quoted in this thesis represent the quintes- 
sence of the results from a very large number of computer runs. 
Ideally, the entire print-out from each run should be included 
in this thesis, so that interested readers would be able to 
investigate such things as what combination of clock and atmos- 
phere parameters were employed and what the parameter correlation 
matrix looked like. Unfortunately the sheer volume of this 
print-out (several thousand pages) precludes its being included 
in this thesis. In order to make this information as available 
as possible, the computer print-outs are being kept on file at 
the Haystack Obser^^atory , and the results reported here are 
catalogued so as to unambiguously identify the computer run 
which produced each result. Each computer run is assigned a 
unique ruji code which identifies it by date and sequence num- 
ber. Thus, for example" ber 75095-32 is the run code 

for. the 32nd run prcv. ~ 9 5th daj of 19 75 

{April 5) . 

The baseline results . .^-v _ :ed in this thesis are summarized 
in Table l. The first entry in each line of these tables is 
the run code, as explained above. The next entry further iden- 
tifies the run by giving the date of the first observations in 
the data set. The next entry identifies the baseline according 
to the following codes: 
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H - Haystack (Tyngstoro- M^ss.) 

G - Goldstone (Mohave, California) 

N - NRAO (Green Bank, W, Virginia) 

S - Onsala (Sweden) 

A - Alaska (Gilmore Creek, Alaska) 

The next entry in each line is the baseline component, in cen- 
timeters. This value is followed by the difference of that baseline 
component value from the weighted mean value (see below) , and 
finally by the formal standard error (sigma) for that component 
from that particular run. The bottom line of each section of 
the table gives the weighted mean and standard deviation of 
the solutions tabulated above, where the weighted mean is 
defined as 

M= (ZP^/w?)/E(l/w?) 
i i 

and the standard deviation* is defined as 
i i ^ 

where is the estimated value of the parameter being displayed, 

Wj^ is the formal error for that estimate, and 

D. = P. - M 
1 1 

This last line, of course, does not list a run code or date. One 
final note — it is possible to display in this table a value 
which one does not wish included in the weighted mean value, 
either because it is not independent of the other results, or 
because something else is wrong with the result. The asterisks 

*'Jhe standard deviation of the mean is smaller by a factor of the square 
rr»at of the numr^xjr of independent measurements, of course. 
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on the far right-hand side of the table indicate which values 
have been included in the weighted mean. 

The source coordinate results reported in this thesis are 
summarized in Table 5. The entries in this table parallel the 
entries in the baseline table. The values for right ascension 
are given in hours, minutes, and seconds of time, and the dif- 
ferences and formal errors are given in seconds of time. The 
values for declination are given in degrees, minutes, and seconds 
of arc, and the differences and formal errors are given in seconds of 
arc. 

In addition to the solution values, the residuals from 
some of these computer runs are plotted in Figures 3 to 23. 

The heading at the top of each figure gives the baseline ac- 
cording to the above codes. The heading then gives the type of 
residual (delay or delay rate) being plotted, and the horizontal 
scale units (nanoseconds for delay, picoseconds per second for 
delay rates). Finally, the heading gives the date and run code. 
The plots show the delays or rates on the horizontal axis. The 

vertical direction is the time of the observation, in days 

starting with 0:00 UTC on the date of the first observation in 

tl data set. Individual observations are plotted with letters, 
the letters indicating the source under observation according 
to the following code; 
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a 

= 

7 * 

= 

3C84 

P 

— 

56 


3C418 
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- 

9 

= 

NRAO140 
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= 
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2134+00 
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60 

= 

0x161 
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= 

66 

= 
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= 

OJ287 
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67 

= 

3C454.3 
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= 

28 

= 

4C39.25 
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= 

OR103 

h 

= 

29 

= 

OK290 

X 

= 

75 

= 

NRA0512 

j 

= 
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3C273B 
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7*’ 

- 

3C390.3 

k 

= 

34 

= 

3C274 

z 

= 

78 

= 

2048+31 

1 

= 

35 

= 

3C275 

a 


79 

= 

0H4 VI 

m 

= 

36 

= 

OQ208 

P 

= 

80 

r= 

OQ172 

n 

= 

41 

= 

CTD9 3 

Y 

= 

81 

= 

4C11.50 

o 

= 

43 

= 

3C345 

6 

= 

82 

= 

CTD135 

P 


53 

= 

OV080 

c 

= 

03 

= 

2344+09A 






C 

= 

84 

s= 

2344+09B 


In addition, data points which have been edited out of the 
solution (see Chapter 4 ) are indicated by a single dot, and 
points which are off scale are indicated by a 0 symbol at 
the edge of the plot. (Note; many of these sources were not 
observed, and some were observed only a very few times.) 


I 

The numbers in this column are sequence numbers in our 
computer zed tables of source coordinates. 
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